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Annealing of Cold Worked Metals 


By Pauu A. Breck 
University of Illinois, Urbana, U.S.A. 


$1. INTRODUCTION 


THE investigations of Taylor and Quinney, Suzuki, Clarebrough, 
Hargreaves, Michell and West), Bever and Ticknor), Gordon and 
others clearly established that of the order of 1% to 10% of the energy 
expended in plastic deformation remains stored in cold deformed metals. 
For instance in heavily cold worked copper the stored energy of cold work 
may amount to approximately 30 calories per mole, or of the order of 1°, 
of the heat of fusion. Even though the stored energy of cold work is 
relatively small, it is an indication of the instability of the deformed state. 
As a result of this instability, the deformed metal has a tendency to return 
to the low energy soft-annealed condition, when the temperature is high 
enough for sufficient atomic mobility. The present paper deals with the 
various processes whereby the well-annealed condition may be approached, 
namely, recovery, subgrain growth* and recrystallization. Earlier work 
in the whole area has been reviewed extensively by Burgers'® and 
Barrett‘; other reviews cover more specifically the field of subgrains‘®? 
or of recrystallization and grain growth”: 1° 11, 12, 13), The present paper 
aims primarily at the presentation of developments resulting from the 
manifold experimental research activities of the last few years. 


* Some investigators prefer to include the phenomenon of subgrain growth 
(or ‘ polygonization ’) under the heading of recovery. However, true recovery 
phenomena, characterized by their kinetics (figs. 2, 6, 7, and 9), do occur 


‘without any interface migration on a scale detectable by present methods. 


Although the kinetics of subgrain growth is not well known as yet, there is 
some indication that it may be different from that of true recovery. Also, 
subgrain growth, as detected by present microscopic or x-ray diffraction 
methods, apparently requires distinctly higher annealing temperatures (or 
longer times) than recovery. In view of the facts mentioned, the separation 
of these phenomena into two distinct classes appears to be well justified. 
Unfortunately, the usage of nomenclature in this field is rather confused ; some 
authors use the term ‘ recovery ’ in such a broad sense that even recrystalliza- 
tion phenomena are included. 
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§ 2. RECOVERY 
2.1. Recovery of Work Hardening 

In single crystals of Zn, when deformed in such a manner that only one 
glide system operates (‘easy glide’), and no bending is involved, the 
work hardening, as measured by the increase of the resolved shear stress on 
the operating slip system, quickly decreases when the deformed crystal is 
stored at room temperature. This phenomenon of recovery (° Kris- 
tallerholung ’) has been first described by Polanyi and Schmid” and later 
studied in more detail by Haase and Schmid), who showed that in 
extended Zn single crystals it gradually proceeded to completion, so that 
the soft-annealed condition was in effect reached within a period of 
approximately one day at room temperature. Even after repeated cycles 
of stressing and recovery, provided that the deformation was glide on a 
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Nine successive shear stress-strain curves from the same pure zinc single 
crystal, showing complete recovery by intermediate annealing treat- 
ments of 1 hour at 260°c. [Ref. (16).] 


single slip system without bending, recovery went to completion at room 
temperature in about one day, so that the stress-strain curve upon every 
new load application was identical with the original one. Similar results 
with Zn single crystals were recently obtained by Li, Washburn and 
Parker“, who repeated the cycle of shear deformation and complete 
recovery in one hour at 260°c as many as nine times. It is characteristic 
of these experiments that after recovery anneal the well annealed single 
crystal condition was again attained, as shown for instance by the return 
to the original flow stress values (fig. 1). No interfaces of any kind were 
introduced into the single crystal either during the deformation process 
or the subsequent annealing. As pointed out particularly by Masing“”), 
the process of recovery is homogeneous on a microscopic scale, and it does 
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not involve any interface migration. Although one may not as yet be 
certain as to the details of the exact nature of the lattice imperfections that 
are responsible for work hardening in pure shear on a single slip system, it 
is clear that these imperfections must be of an easily recoverable type. It 
is significant that under such conditions of deformation Honeycombe(#) 
found no x-ray asterism in single crystals of Cd, even after 100° extension, 
or more. If the deformation does involve bending and x-ray asterism is 
observed, the resulting work hardening is greater, but only a fraction of 
this work hardening is recoverable’), This is particularly true for 
polycrystalline metals, where the constraints of neighbour grains force each 
crystal to deform in a complex manner. The point is well illustrated by 
the results of Tietz, Anderson and Dorn®>) with polycrystalline high 
purity aluminium specimens deformed in tension at three different 
temperatures to the same flow stress of 12 300 p.s.i., as measured at 
78°K. The recovery of the flow stress after various periods of annealing at 
305°K was measured at 78°K. The extent of recovery at 305°K was 
approximately 5-7%, 10% and 14% of the total work hardening for the 
specimens deformed at 273°, 194° and 78°K, respectively. These results 
show that the easily recoverable component represents a variable fraction 
of the constant total work hardening, and that this fraction decreases with 
increasing temperature of deformation. By extrapolation, one may 
expect that plastic deformation of polycrystalline high purity aluminium 
at higher temperatures, suchas 60°c or 100°C, may result in work hardening 
with no recoverable component at all. 

Cherian, Pietrokowsky and Dorn found that polycrystalline tension 
specimens of commercial purity Al after small tension strains at room 
temperature show very marked recovery of the yield stress on annealing 
at temperatures up to 100°c, without any lowering of the flow stress at 
higher strains (‘metarecovery ’). They reported that annealing in the 
temperature range of about 150°-200°c results, not only in a lowering of 
the yield stress, but also in a lowering of the flow stress at high strains 
(‘ orthorecovery ’). The two phenomena appear to have different 
kinetics. It is to be noted, however, that in their study of the recovery 
of zine single crystals deformed in pure shear, Parker and his associates 
did not observe metarecovery. Also, there appears to be no indication of 
the occurrence of metarecovery in the study by Dorn and his associates of 
recovery in pure Al'25), It is very likely that the occurrence of * meta- 
recovery ’, observed with Al of commercial purity, is not a general 
phenomenon, but that it is associated with the impurities in the particular 
metal tested by Cherian, Pietrokowsky and Dorn. 


2.2. Kinetics of the Recovery of Work Hardening 
Drouard, Washburn and Parker‘!®) measured isothermal recovery at a 
series of different temperatures in zinc single crystals deformed at 78°K in 
pure shear along a single glide system, without bending. The amount of 
recovery was measured by the decrease of the flow stress in percents of the 
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difference in flow stress between the deformed and the fully annealed 
condition. Isothermal curves for various temperatures are shown in fig. 
2. The data clearly show the type of kinetics typical of recovery: at all 
temperatures the rate of the process is highest initially, and then it 
decreases with continued recovery. It is to be noted that the kinetics of 
recovery is significantly different from that of a process involving nucle- 
ation and interface migration, where the rate is at first small, then it 
increases to a maximum and decreases again (see fig. 42). This difference 
further supports the conclusion in the previous paragraph, that recovery 
does not involve the migration of interfaces. 

The isothermal kinetics of recovery is shown more fully in fig. 3, where 
the recovery, expressed in percents of the total recoverable work hardening, 
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Residual strain hardening of pure Zn single crystals as a function of annealing 
time at various temperatures, after 8%, shear strain at —50°c. [Ref. (19).] 


is plotted against the logarithm of the temperature-adjusted recovery 
time, combining all of Drouard, Washburn and Parker’s data for various 
temperatures by means of a single activation energy. The central part of 
the curve has an inflection point and it may be considered as a straight line 
over a range corresponding to a time ratio of about 250. At short 
annealing times the experimental recovery values are larger and at long 
annealing times they are smaller than those extrapolated from the 
straight line relationship. Very similar kinetics were found previously by 
Masing and Raffelsieper’?” for pure Al single crystals deformed in tension 
at room temperature. They considered the logarithmic recovery kinetics 
(corresponding to the straight line central portion in fig. 3) as a confirma- 
tion of the theory of recovery by Kuhlmann, Masing and Raffelsieper (2?) 

As shown by Kuhlmann’), the assumption that the activation energy ion 
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recovery decreases linearly with progressive recovery leads directly to a 
linear relationship between recovery and the logarithm of the annealing 
time. The deviations from this relationship at very short and very long 
annealing times are also accounted for by the theory of Kuhlmann, Masin g 
and Raffelsieper. 


Fig. 3 
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Fraction of strain hardening recovered (f) as a function of the logarithm ot the 
temperature-adjusted annealing time, from the data in fig. 2. [Ref. (19).] 


2.3. Activation Energy of the Recovery of Work Hardening 

According to the results of Drouard, Washburn and Parker), the 
temperature dependence of the rate of recovery of Zn single crystals 
sheared at 78°K and annealed in the temperature range of —20°c to 
+10°c can be described satisfactorily by means of a single activation 
energy of 20 kilocalories/mole (fig. 4). This value is in good agreement 
with that for the self-diffusion in Zn crystals along the hexagonal axis”), 
suggesting that the elementary process in recovery may be similar to that 
in diffusion. However, the independence of the activation energy from 
the extent of recovery in these experiments appears to contradict the 
theory of Kuhlmann, Masing and Raffelsieper. The experimental data of 
Raffelsieper with Al, which indicated a linear decrease of the activation 
energy with progressive recovery (fig. 5), were obtained under less 
favourable conditions than those of Drouard, Washburn, and Parker, 
since in the extended Al single crystals a substantial portion of the work 
hardening was not recoverable. 
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2.4, Recovery of the Stored Energy of Cold Work 


It is significant that the kinetics of the low temperature (60° and 100°c) 
isothermal release of the stored energy of cold work from cold rolled 
polycrystalline Cu and Al, according to Borelius, Berglund and Sjoberg’ 
also obeys a logarithmic law, since the rate of evolution of heat varies 
inversely with the annealing time (see fig. 6 and eqn. (1)). 
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Temperature dependence of the annealing time required for various fractions: 
of recovery in sheared zine single crystals (data from fig. 2). The slope 
of the lines corresponds to a constant activation energy of 20 keal per- 
mole. [Ref. (19).] 
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where P is the rate of energy release at the time ¢. This result establishes 
a close analogy between the energy release and the recovery of flow stress, 
as discussed above, and suggests that the reaction associated with the 
energy release studied by Borelius and his associates is also recovery, not 
involving the migration of interfaces—a very probable inference, con- 
sidering the low temperatures and the short annealing periods involved. 
It would be interesting to have flow stress and calorimetric recovery data 
for identical specimens, in order to ascertain whether or not the time. 
constants for the two recovery curves are the same. Further questions 
that arise in this connection concern the effect of the temperature of 
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Fig. 5 
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Decrease of activation energy for recovery of strain hardening in Al single 
erystals with progressive recovery. [Ref. (22).] 
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Rate of release of stored energy of cold work at 100°c from Cu rolled to double 
length, as a function of annealing time, and in inset reciprocal annealing 
time. [Ref. (26).] 
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deformation on the release of the stored energy of cold work. In analogy 
with the effect of temperature on the recovery of yield stress, one may 
expect, for instance, that with decreasing deformation temperature the 
rate of energy release and also the total energy released during recovery 
increase. Also, it is quite possible that plastic deformation at somewhat 
higher temperatures (60°-100°c for Al) may not lead to the storage of any 
recoverable energy, so that the energy stored under such conditions could 
be released only in recrystallization. 


2.5. Recovery of X-ray Line Broadening 

The kinetics of the isothermal sharpening of x-ray diffraction lines 
broadened by cold working has been studied by Van Arkel and Burgers®”? 
in the case of tungsten. They found that in cold drawn thoriated tungsten 
wires considerable line sharpening takes place prior to recrystallization, 
and that the kinetics of this process is qualitatively similar to that illustrated 
above for recovery in the case of flow stress and energy release. It has 
been held by some investigators that the recovery of line broadening prior 
to recrystallization occurs only in metals of high melting point, such as Mo 
and W. In other metals, for instance in cold worked Cu, Ni, alpha brass, 
and steel, Wilson and Thomassen'?*) found that line sharpening on 
annealing occurs simultaneously with softening. Averbach®®) found that 
in brass filings line sharpening was a result of recrystallization. rather 
than of a recovery process prior to and distinct from recrystallization. 
However, a recent study of isothermal line sharpening in cold rolled pure Al 
single crystals by Lutts and Beck) gave definite separation between line 
sharpening and the major portion of softening (fig. 24); it also indicated 
clearly a typical recovery-type kinetics for line sharpening. The kinetics 
itself is, of course, not sufficient to prove that the observed line sharpening 
in effect takes place by a recovery process in the strict sense, without 
sub-boundary migration. 

Undoubtedly, ‘recovery ’ of line broadening prior to recrystallization 
may occur in low-melting metals, as well as in those of high melting point. 
Whether or not it may be actually observed under certain conditions 
apparently depends mainly on the relative rates of recovery and of 
recrystallization. In the particular experiments with Al, referred to 
. above, recrystallization, which presumably depends on the availability of 
‘nuclei ’ in suitable orientations, was delayed by using single crystals of 
the (110) [112] orientation and by careful techniques, eliminating contami- 
nation of the resulting single orientation rolling texture by crystal 
fragments of other orientations. Recovery, which does not involve 
nucleation or interface migration, is less likely to be affected by such 
procedures. The relative rates of recovery and of recrystallization 
depend also on the presence of finely dispersed particles of a second phase. 
This may be seen for instance in the results of Van Arkel and Burgers”, 
who found much more distinct separation of the two phenomena in 
thoriated tungsten than in pure tungsten. The finely dispersed thoria 
particles slow down the migration of grain boundaries necessary for 
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recrystallization, while their effect on the rate of recovery is apparently 
smaller, if any. In view of the results with brass mentioned above, 
solid solution of Zn in Cu apparently does not tend to delay recrystalliza- 
tion more than recovery. 

Just as the recoverable component of work hardening, the line 
broadening in the cold worked state also increases with decreasing 
deformation temperature, according to the investigations of Wilson and 
Thomassen®Y and of Paterson and Orowan®), It is a reasonable 
assumption that all of the line broadening is, at least potentially, recover- 
able (may be annealed out by a recovery mechanism, provided that 
recrystallization does not set in prematurely). The results of Paterson 
and Orowan indicate that the variation of the amount of line broadening 
with the deformation temperature cannot be explained on the basis of 
increasing recovery at higher deformation temperatures, but that it must be 
considered as a true effect of the temperature of deformation. It is 
therefore an interesting question whether or not the much larger line 
broadening observed in W, as compared with Al, may be accounted for by 
the higher rate of recovery at the deformation temperature in the lower 
melting metal, as suggested by Williamson and Hall), 

2.6. Recovery of Electrical Resistivity 

In view of the recent publication of an extensive review by Broom‘®), 
the present remarks will be restricted to a brief summary of newer 
results. Until recently it was thought that the increase in electrical 
resistivity resulting from plastic deformation of polycrystalline copper 
-does not recover to a measurable extent prior to recrystallization. This 
type of behaviour was indicated for instance by the work of Smart, Smith 
and Phillips®*) with high purity copper wire cold drawn nominally at 
room temperature (but due to the heat evolved in the die, actually probably 
at about 60° to 100°c). However, Molenaar and Aarts‘*) discovered 
recently that, if the plastic deformation is carried out at a low temperature, 
the increase in resistivity is larger and that it recovers to an easily 
measurable degree at temperatures at which no recrystallization, and in 
‘fact no measurable change in the mechanical properties is taking place. 
The increase of the effect of plastic deformation on the resistivity with 
decreasing temperature of deformation has been recently demonstrated 
more fully for several metals®®; the effect for Ni is shown in fig. 7. 
Whether the total amount of this increase is recoverable, or the non- 
recoverable part of the resistivity increment also becomes larger with 
decreasing deformation temperature, is not as yet known. R. R. 
Eggleston”) found (fig. 8 (a)) that for copper deformed at 4°xK the electrical 
resistivity undergoes considerable recovery when annealed for 1 min 
periods at temperatures as low as —150°c to 0°c and that when the 
temperature rises to about 100°c this recovery (in view of the kinetics and 
the low temperatures and short annealing periods involved, it is extremely 
unlikely that interface migration processes would play a role here) 
amounts to approximately 50% of the total resistivity increment originally 
‘produced by the deformation (fig. 8 (a)). Presumably, in his experiments 
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the remaining non-recoverable portion of the resistivity increment could 
have been eliminated by recrystallization on annealing at still higher 
temperatures, as in the work of Smart, Smith and Phillips, if such annealing 
treatments had been carried out. In analogy to the effect on the recover- 
able component of work hardening in aluminium, apparently plastic 
deformation at increasingly higher temperatures also has a decreasing 
effect on the recoverable component of the resistivity increment until, as 
mentioned above, none of the resistivity increment is recoverable prior to 
recrystallization in the case of pure Cu deformed at about 60° to 100 BLE 
The results of Molenaar and Aarts‘?>) indicate that the rate of recovery 
of the electrical resistivity is greater than that of the flow stress, suggesting 
that these two properties may be affected by different types of imper- 


Fig. 7 
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Increase in electrical resistivity in Ni as a function of true strain at 100°c (empty 
circles), 0°c (crosses), —78-5°c (filled circles) and —183°c (squares). 
[Ref. (36).] 
fections resulting from cold work. Nevertheless, the low temperature 
drop of electrical resistivity shows kinetics”) (fig. 8 (6)) similar to that of the 
recovery of flow stress, and of line broadening and of the low temperature 
release of the stored energy of cold work. Furthermore, all these property 
changes take place under conditions of annealing time and temperature 
where no interface migration is observable by microscopic or x-ray 
diffraction methods, so that they all may be ascribed to recovery. As has. 
been already observed by Tammann®*), the recovery of different properties 
may take place at quite different rates, so that a cold worked metal may 
display a whole spectrum of recovery phenomena.* If these correspond to. 


ST eIeeE Se MTEL Tees wie eat iy 
* This is apparently true, even though i i igati 
rue, gh in some of the early investigations. 
the observed differences in rate may have been due to dimartince in iachantaty 
since no clear separation was achieved between recovery and recrystallization. ‘ 17) 
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different elementary processes, it may be that further research will reveal 
not only variations in the quantitative details of the kinetics, but also 
differences in the activation energies. On the other hand, if the various 
recovery phenomena correspond to progressive stages of a single process, 
their activation energies should be the same. 
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(6) 
Recovery of electrical resistivity as a function (a) of annealing temperature, 
and (b) of annealing time at various temperatures, after deformation by 
twisting at 4-2°x. [Ref. (37) and unpublished work by R. R. Eggleston. |] 
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The main results on recovery to date may be summarized as follows : 

(a) In general, fractions of certain property changes resulting from 
plastic deformation, for instance of work hardening, may be recovered 
in annealing periods and at temperatures where no interface migration 
can be detected by methods now available. In the case of single crystals 
deformed in pure shear, without bending, recovery may be carried to 
completion and the absence of any interface migration or reorientation in 
the recovered condition is particularly clearly shown by the accuracy with 
which the original mechanical properties can be re-established. 

(b) The recoverable fraction of work hardening, of electrical resistivity 
increment, and presumably also of other property changes, increases with 
decreasing deformation temperature. The amount of x-ray line 
broadening, all of which is apparently fundamentally recoverable, also 
increases with decreasing deformation temperature. 

(c) The rate of isothermal recovery of each property change is highest 
at.the start, and it gradually decreases as recovery progresses. At least 
within a certain range, the amount of recovery of the flow stress and of the 
stored energy of cold work is a linear function of the logarithm of the 
annealing time, so that, within this range, the rate of recovery is inversely 
proportional to the annealing time. 

(d) The kinetics of the recovery of various property changes appears to 
be very similar, as far as the available data dicate. However, the rate 
of recovery may be quite different for each property. Thus, apparently 
the recovery of electrical resistivity is faster than that of x-ray line 
broadening, and this in turn is faster than the recovery of hardness. 


§ 3. SUBGRAIN GROWTH 
3.1. Subgrains and Subboundaries 


Microscopic evidence for the existence of subboundaries or veining in 
ferrite deformed at high temperature has been reported by Rawdon and 
Berglund®® as far back as 1927-28. Greninger‘*) concluded from x-ray 
diffraction evidence (back reflection Laue patterns) that veining 
corresponds to a subdivision of the ferrite grains into individual crystal 
blocks slightly disoriented relative to each other. That the microscopic 
subboundaries in effect correspond to the limits of disoriented domains in 
aluminium was demonstrated more recently by Guinier and Lacombe, 
Collins and Mathewson‘*”) discovered the appearance of slightly disoriented 
crystal blocks upon annealing at 500°c of Al crystals previously extended 
at room temperature; the Laue asterism streaks, which were continuous 
after deformation, split into a multitude of more or less well defined dots 
upon annealing. The information on subgrains up to about 1950 was 
reviewed by Guinier'®), 


3.2. Subboundary Migration 


Crussard’s early results“ already indicated that the dots replacing the 
original asterism streaks in extended Al crystals on annealing at about 
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350° to 400°c become more sharply defined and less numerous on continued 
annealing at higher temperatures, particularly above 600°c, suggesting a 
decrease in the number of crystal blocks, i.e. an increase in the average 
subgrain size. Also he was the first to note that in Al crystals deformed in 
tension the subgrain growth on annealing is associated with softening, 
i.e., a decrease in yield stress. This annealing phenomenon was designated 
by Crussard as ‘recrystallization in situ’, meaning softening without 
reorventation. However, in the absence of reorientation, that is the migra- 
tion of high angle grain boundaries, it is perhaps better not to use the 
term ~recrystallization’*. For this reason, in the present paper the 
term “subgrain growth’ has been adopted to describe this stage of the 
annealing process. That subgrain growth is apparently not always 
associated with appreciable softening is suggested by the observations of 
Honeycombe"®). The movements of subboundaries in aluminium on 
annealing at around 600°c were studied metallographically by Lacombe 
and co-workers’: 4°), and by means of a high resolution focusing Laue 
method by Guinier and Tennevin'*®). The latter investigators found that 
at the annealing temperatures used the subgrains tend to grow, and that 
this process is much faster in high purity aluminium than in commercial 
metal. The work of Dunn and Daniels‘*”) provided convincing metallo- 
graphic and x-ray evidence for subgrain growth in plastically bent silicon— 
iron crystals on annealing in the temperature range of 950° to 1300°c. 
(figs. 9and 10, Plates 25 and 26). These workers showed that the process of 
subgrain growth (increase of the average subgrain size with annealing 
time and temperature) leads to a decrease in the total subboundary surface 
energy per unit volume. They pointed out that, as in ordinary grain 
growth, interface energy may be considered as the driving energy of sub- 
grain growth. 

Whether or not the mobility of all very low angle boundaries is in fact 
particularly high in comparison with grain boundaries of intermediate 
disorientation, as suggested by Dunn and Daniels, could be definitely 
decided only by experiments in which the rate of migration of various 
boundaries is observed in specimens of the same composition, at the same 
temperature and under the influence of the same driving force. There 
certainly appears to be a significant difference between subgrain growth 
and ordinary grain growth at least in regard to the effect of an external 
stress applied during annealing. Wood and Scrutton ® discovered that 
the rate of subgrain growth is greatly accelerated by an externally applied 
stress, sufficient to produce a small creep strain during the annealing. 
This effect was later confirmed by Hino, Shewmon and Beck‘**) by means 
of both x-ray diffraction and a microscopic method. If the tension yield 
stress is accepted as an indicator of subgrain size, as discussed further 
below, the data of Sherby, Goldberg and Dorn show an example of 
acceleration of subgrain growth by a factor of over 100 when a creep 
strain of 0-004 is applied during the annealing of cold worked pure Al at 
530°K. On the other hand, Sperry and Beck'* were unable to detect an 


258 Paul A. Beck on the 


analogous effect of simultaneous stressing on the rate of ordinary grain 
growth in high purity aluminium. This point is significant since it proves 
that the often irregularly shaped subboundaries formed in polycrystalline 
metals on plastic deformation (fig. 11, Plate 27) are in their behaviour more 
closely related to the crystallographically well defined individual kink 
boundaries in zinc single crystals, which Parker and his associates‘? °*) 
have shown to be mobile under external stress (fig. 12, Plate 28) than to 
ordinary grain boundaries, which are not. The difference in the 
behaviour of subboundaries and of ordinary grain boundaries is probably 
not only due to the difference in the extent of disorientation, but more 
significantly also to the type of disorientation‘), as described in the 
following. Ideal subboundaries of the type described by Burgers®®, 


Fig. 13 
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Mobility of dislocation type boundaries in pure Zn single crystals as a function 
of disorientation at the boundary. (Rate of migration at 350°c under 
a constant shear stress of 9-19 p.s.i.) [Ref. (58).] 


may be constructed of a single set of parallel edge dislocations and the 
orientation relationship between the two adjoining lattices corresponds to 
a rotation around an axis parallel to the rational crystallographic direction 
of the edge dislocations. Such a boundary is highly mobile under stress 

even at low temperatures‘: 5”), and it may be expected to show relatively 
high mobility also under the influence of interfacial tension, without 
externally applied stress. Actual measurement of the mobility of 
boundaries of this simple type as a function of disorientation was recently 
carried out successfully by Bainbridge, Li, Edwards, Washburn and — 
Parker(>®), As seen in fig. 13, the mobility gradually decreases with 
increasing disorientation and shows no sign of reversal of this trend as far 
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as the measurements were extended, namely up to about 16°. On the 
other hand, when the orientation relationship of two adjoining lattices is of 
&® more general nature, even though the disorientation may be still small, 
at least two nonparallel sets of edge dislocations are required to form the 
boundary. In this case, movement of the boundary requires atomic 
readjustments that can be achieved only by diffusion at the interface), 
even if impurities that pin down the dislocations are disregarded. Conse- 
quently, the rate of boundary migration is relatively low, particularly 
at small disorientations, where the boundary atoms are tightly held in 
their positions. (For a discussion of the orientation dependence of grain 
boundary mobility see § 4.5 (b).) 


3.3. Mechanism of Subgrain Formation 


As a result of the investigations of Wood and his co-workers(48> 69: 61, 62) 
it appears to be generally accepted that in creep subgrains form directly 


Distribution of dislocations in a bent crystal (a) before polygonization, (b) after 
polygonization, according to Orowan and Cahn. [Ref. (68).] 


during the slow deformation at high temperature, and not as a result of 
annealing subsequent to it. On the other hand, considerable discussion 
has taken place in recent years concerning the origin of subgrains in metals 
deformed at low temperatures and then annealed. In view of the work of 
Cahn‘: &4), who investigated the subgrains formed in cold bent and 
annealed single crystals of Zn and Al, Orowan and Cahn concluded that the 
subboundaries are formed during annealing by means of the thermally 
activated migration of individual dislocations, which were left behind in 
the plastically bent crystal. On this picture, the continuously bent 
crystal lattice undergoes sectionwise straightening during annealing, as the 
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randomly dispersed excess dislocations of one sign rearrange themselves 
into dislocation walls, or subboundaries, where the orientation change 
becomes concentrated (fig. 14). A theoretical treatment of this annealing 
process, named by Orowan ‘ polygonization ’, has been given by Cottrell'®*) 
who showed on the basis of the dislocation theory that the arrangement of 
edge dislocations one underneath the other, as in a dislocation wall, 
corresponds to increased stability. In view of the success of these investi- 
gators in describing the geometry of subboundaries in cylindrically bent 
single crystals in terms of dislocation walls, the polygonization mechanism 
has been quite generally accepted for the formation of subgrains in metals 
deformed at low temperatures and then annealed. From this point of 
view, it seemed reasonable for Greenough and Smith‘*®, for instance, to 
reinterpret even the subgrain formation in creep as a two-stage process, 
where the intermediate (continuously bent) condition is not detectable only 
because of the occurrence of thermally activated ‘ polygonization ’ 
concurrently with the creep process itself. 

However, Wood'®”) pointed out that there was no compelling empirical 
reason for such an interpretation, and he expressed the view *®) that, on 
the contrary, subgrains may form directly during plastic deformation even 
at low deformation temperatures. Both Guinier“ and Cahn‘®®) noted 
that, by means of the focusing Laue method of superior angular resolution 
introduced by Guinier and Tennevin, the latter investigators were able to 
detect subgrains at a lower annealing temperature (from 450°C up in the 
case of pure Al) than with the conventional Laue technique used by Cahn 
(near 600°C). It may be concluded that in all probability the subgrains 
were already present at an earlier stage of the annealing process, where 
Cahn’s Laue asterism streaks still appeared continuous because of 
inadequate resolution. Guinier and Tennevin also pointed out'*) that 
the focusing Laue method’s limit of resolution is about 1/3 mm, and that 
subgrains smaller than that would give continuous Laue streaks rather 
than resolved individual dots in the x-ray patterns. In fact, further 
improvement in resolution, achieved by Lambot‘” through the use of a 
convergent monochromatic x-ray beam focused at a small area of the 
specimen surface, lead to the detection of a recognizable substructure in 
extended Al crystals at a still lower annealing temperature (200°C) and 
~ even wtihout any annealing above room temperature (fig. 15, Plate 28). 
As a matter of fact, even in ordinary back reflection Laue patterns, using 
a fine beam, Honeycombe®) observed the breakdown of the asterisms 
into several intensity maxima for Al crystals immediately after 2—6°% 
elongation at room temperature. The presence of a substructure with 
quite sharp domain boundaries in pure Al deformed at room temperature 
and not annealed was furthermore convincingly demonstrated by 
Heidenreich” using an electron transmission microscopic method (fig. 16, 
Plate 29). Kellar, Hirsch and Thorp‘*”) found that back reflection Debye 
rings taken of cold worked Al with a microbeam x-ray diffraction technique 
can be resolved into individual dots. Both of these investigations, and 
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also the polarized ligth microscopic observations of Beck and Hu‘) gave 
an approximate subgrain size of 2 microns for pure Al severely deformed 
at room temperature. No indication of substantial subgrain growth was 
found? on annealing up to 200°c. Heidenreich found that sub- 
structure was also present in Al deformed at 78°xK and examined at room 
temperature, although it was distinctly finer than that in the metal 
deformed at room temperature, namely about 1 micron. More recent 
microbeam X-ray diffraction work by Gay and Kelly'’® succeeded in 
detecting the presence of a substructure in polycrystalline Cu, Ni and Fe 
deformed at room temperature and unannealed (fig. 17, Plate 29). Chen 
and Maddin‘®) found that even a metal with as high a melting point as 
Mo, shows detectable substructure when a single crystal is deformed at 
room temperature. Conard, Averbach and Cohen‘’® were able to show 
by means of x-ray diffraction with the specimen held at 78°x, that Zn 
single crystals deformed at 78°K and not allowed to heat up above that 
temperature, have a distinct substructure. These observations clearly 
suggest that substructures are formed on plastic deformation under con- 
ditions where thermally activated processes can not play a significant role. 

As mentioned above, the interpretation of subboundaries as ‘ dislocation 
walls ° or arrays of parallel edge dislocations by Orowan, Cahn'®) and 
Cottrell'®®) was originally associated with the view that such boundaries 
are formed an annealing by ‘ polygonization’, i.e., by means of the 
sectionwise straightening of the crystal lattice continuously curved in 
plastic deformation. Polygonization was viewed as the thermally 
activated migration of individual dislocations from their random positions 
in the deformed crystals and the collection of these dislocations along 
planes nearly perpendicular to the active slip direction. In view of the 
results described in the previous paragraph, it now appears quite certain 
that dislocation walls are actually formed during the deformation process 
itself, even at very low deformation temperatures. It is, of course, not 
necessary to assume that al// the excess dislocations of one sign are necessarily 
collected into dislocation walls during the deformation. The lack of sharp 
extinction contour lines’) and the width of the individual x-ray micro 
beam diffraction spots'’”) suggest that some plastic strain may be present 
within the subgrains themselves in the cold worked condition. Neverthe- 
less, the fact, shown by the various experiments mentioned above, that 
metal crystals are in effect ‘fragmented’ into rather distinct lattice 
domains of differing orientations directly on plastic deformation, 
apparently without requiring much or any thermal activation, suggests 
that a mechanism other than polygonization is involved. Whether or not 
some additional subboundaries actually do form on subsequent annealing 
by polygonization can not be decided with certainty at the present time. 
This question is intimately related to the as yet indefinite interpretation 
of the ‘ background’ darkening between the separated spots in X-ray 
diffraction patterns, such as fig. 17. This background apparently. 
originates from diffraction by crystal regions of orientation intermediate 
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between those giving rise to the separated diffraction spots. These 
intermediate orientations may possibly be represented by distinct 
subgrains too small to be resolved with the diffraction method used, or by 
continuously bent crystal regions. In the event that the former inter- 
pretation should prove correct, the apparent continuity of the background 
would be eliminated by subgrain growth. If the latter interpretation be 
true, the disappearance of the continuous background on annealing may 
result from the movement of excess dislocations of one sign from such 
bent regions either to subboundaries previously formed in the vicinity 
during deformation or to subboundaries newly formed on annealing by 
polygonization. 

As to the mechanism of subboundary formation on plastic deformation, 
it is significant that the deformation mechanism known as ‘ kinking ’ 
according to Orowan‘®), Jillson‘®), Hess and Barrett‘*?, Wsahburn and 
Parker?) and Li, Edwards, Washburn and Parker‘®®), results in the 
formation of dislocation walls identical with those one expects to form by 
thermal activation on the basis of the postulated polygonization 
mechanism. In view of the known mobility of edge dislocations under the 
effect of shear stress even at low temperatures, and of the stability of their 
arrangement in the form of dislocation walls according to Cottrell®, it 
seems rather reasonable to assume that a large majority of the excess 
dislocations of one sign will so align themselves during the deformation 
process itself. Recently Bainbridge, Li, Edwards, Washburn and 
Parker©® observed that two parallel kink boundaries in a Zn single crystal, 
migrating on application of an external shear stress, may join to form a 
single kink boundary of greater disorientation, even when the crystal is 
strained at a temperature as low as 78°k, and that a 1° kink boundary may 
trap additional dislocations, which pass through the crystal on straining, 
and thereby increase the boundary angle to as high as 15°. On the basis 
of the various observations described, it appears quite probable that in 
metals the subboundaries form on plastic deformation essentially by the 
kinking mechanism—an athermal process. 

The formation of a kink boundary requires slip more or less evenly 
distributed over a large number of parallel slip planes. For a subboundary 
with a disorientation of 3° the required average distance of the active slip 
planes in Al would be of the order of 0-02 y, that is much smaller than the 
observed distances of ordinary slip bands (order of 1-10 »). However, this 
requirement represents no difficulty ®® for the proposed kink mechanism 
of subboundary formation, since it has been known already from thé early 
work of Masing and Polanyi'*”) that single crystals of Sn can be sometimes 
extended 1000°% without the formation of any visible slip bands. That 
such ‘ fine slip’ in effect may occur to a substantial degree even in the 
presence of slip band formation, was shown recently for Al by Kuhlmann- 
Wilsdorf and Wilsdorf‘*), The above mentioned results of Parker and 
his associates on the movement of kink boundaries under externally 
applied stress in single crystals of Zn provide additional proof for the 


Annealing of Cold Worked Metals 263 


occurrence of fine slip; the movement of a dislocation wall under such 
conditions can take place only by that mechanism. 


3.4. Factors Affecting the Subgrain Size 

Although the formation of subboundaries is apparently essentially 
athermal, as discussed, nevertheless there is a fairly clear indication that a 
thermally activated process is somehow involved at least in determining 
the subgrain size. Wood and Rachinger‘®?) have shown that the subgrain 
size resulting from creep in pure polycrystalline Al increases with increasing 
deformation temperature from 20° to 350°c and with decreasing strain 
rate from 10°, per min to 0-1°% perh. Such an effect of the deformation 
temperature was also found by Heidenreich‘’*) between 78° and 300°K, at 
much higher strain rates than those used in the experiments of Wood. 
The information available concerning this effect is not sufficiently 
quantitative to allow detailed interpretation at the present time. In any 
event, the phenomenon of subgrain growth, discussed above, definitely 
depends on thermal activation ; the results of Dunn and Daniels‘4”) are 
particularly clear cut in showing the increase in subgrain size on annealing 
both with increasing time at constant temperature and with increasing 
annealing temperature. 

The effect of solute atoms on the subgrain size formed on plastic deform- 
ation, and particularly on subgrain growth, are not well known as yet. 
Heidenreich observed’) that 4° of Cu dissolved in Al greatly reduces the 
subgrain size formed on deformation at room temperature. Hazlett and 
Parker‘** noted that the addition of 1 atomic per cent Ti to Ni refined the 
substructure under otherwise similar conditions. One may expect that 
the solute atoms may be also quite effective in preventing the migration of 
dislocation walls as they are known from the work of Cottrell) to im- 
mobilize individual dislocations. On this basis, a very great decrease in 
the rate of subgrain growth would be expected to result from the presence 
of a relatively small number of solute atoms. There appear to be no data. 
available at the present time to test this expectation. 


3.5. Effect of Substructure on Mechanical Strength 


There is now a rather convincing body of evidence'*® indicating that 
the presence of substructure has a great effect on the mechanical properties. 
As shown in fig. 18, Li, Washburn and Parker® found that in a Zn single 
crystal the flow stress in shear is greatly increased if a substructure is 
introduced by a ball indentation at 350°c. Annealing at 400°c does not 
eliminate this large scale substructure (fig. 19, Plate 30), although other 
types of imperfections, which cause work-hardening when Zn crystals are 
sheared without bending, are known to be annealed out readily at much 
lower temperatures(15. 16), Since the high shear flow stress was actually 
observed after a 400°C anneal, it must have been a result of the sub- 
structure and not of the other types of imperfections mentioned. In their 
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introducing substructure by indentation in the [0001] direction at 350°c. 
Specimens were annealed at temperatures indicated and cooled to room 
temperature in 15 minutes, before testing. [Ref. (16).] 
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Flow stress at 0-0025 strain for pure Ni and for a Ni-Ti alloy as a function of 
prestrain, Specimens were prestrained at room temperature and 
annealed for | hour at 800°c before testing. [Ref. (84).1 
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recent work with Ni and with Ni-Ti solid solution alloys, Hazlett and 
Parker‘S found that a few per cents of prestraining in tension followed by 
annealing at 800°C (which is presumably sufficient to recover work 
hardening due to imperfections other than subgrains) greatly increases the 
yield stress in a subsequent tension test. For instance, prestraining to 5%, 
raises the yield stress for pure Ni from about 3500 p.s.i. to 18000 p.s.i. 
(fig. 20), apparently as a result of decreasing subgrain size. Figure 21 
shows the variation of yield stress as a function of the ‘ index of amount of 
substructure ’ for the same specimens. Sherby, Goldberg and Dorn‘5) 
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Same flow stress data as in fig. 20, plotted as a function of the density of 
subboundaries resulting from prestraining to various amounts and 
annealing for 1 hour at 800°c. [Ref. (84).] 


obtained data for pure Al, which may be also interpreted as an illustration 
of the effect of subgrain size on the nonrecoverable component of work 
hardening. Tension specimens were prestrained under varying conditions, 
so that substructures of varying fineness were presumably obtained. The 
recoverable component of work hardening was eliminated either by 
carrying out the prestraining at temperatures above 422°K, or by annealing 
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for 2 hours at 530°K, after prestraining at lower temperatures. It may be 
assumed that the remaining work hardening in these specimens was 
practically all due to the subboundaries. The observed increase of yield 
stress in tension tests at 78°K from about 1300 p.s.i. to 7400 p.s.i., for 
instance with decreasing temperature of prestraining (fig. 22), may be 
ascribed to the progressive decrease in subgrain size. In other experi- 
ments, specimens prestrained at 78°K and presumably possessing a very 
fine substructure were annealed under varying conditions. In view of the 
effect of such conditions on the substructure, as known qualitatively from 
the previous investigations of Wood and his associates, the resulting 
changes in mechanical properties may be again interpreted as consequences 
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Flow stress for pure Al at 0-01 strain at 298°xK after prestraining 15°% at tempera- 
tures shown and annealing for 2 hours at 530°K. The increase in flow 
stress with decreasing temperature of prestrain is presumably due to 
the increasing density of subboundaries, since the annealing probably 
removed all recoverable work hardening. [Replotted from Ref. (50).] 


of variations in subgrain size. Thus, fig. 23 shows the decrease in tension 
yield stress as a function of the period of annealing at 530°K, presumably 
asa result of subgrain growth. Although the above interpretation appears 
to be a very probable one, it clearly lacks definitiveness for two reasons : 
(1) the absence of adequate subgrain size measurements, and (2) lack of 
proof of the absence of lattice imperfections other than subboundaries. 


It has been observed microscopically by Beck and Hu‘®®) that subgrain 


growth occurs to a significant extent in the annealing of polycrystalline 
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high purity Al, and that it serves as a mechanism of annealing by which the 
orientation of the deformed material may be partially retained. Recently, 
some instances became known where complete softening in highly 
deformed polycrystalline metals occurs essentially without reorientation. 
Thus, Semchyshen and Timmons reported‘87) that in rolled Mo sheet, the 
texture of which has substantially a single component of orientation, this 
orientation is retained without appreciable change on soft annealing, 
although the material then consists of newly formed, microscopically well 
resolved grains. Since the retainment of a single component texture on 
annealing excludes ordinary high angle grain boundary migration, 
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Flow stress for pure Al at 0-01 strain at 298°K after prestraining 15% at 78°K 
and annealing at 530° for periods indicated. The decrease in flow 
stress with continued annealing is presumably due to decreasing density 
of subboundaries attained by subgrain growth during annealing. [Re- 
plotted from Ref. (50).] 


Semchyshen and Timmons’ results may be best interpreted as softening 
by means of low angle boundary migration, i.e., subgrain growth. It can 
be expected that, in the course of this process, the disorientation between 
neighbouring subgrains gradually becomes larger, and the character of the 
interfaces more and more approaches that of ordinary grain boundaries in 
a highly oriented structure, as for instance in Cu with cube texture. In 
such cases subgrain growth may change over into ordinary grain growth 
quite gradually and it may be rather difficult to separate these processes. 
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Another instance apparently similar to this was observed by McHargue'*» 
with rolled Ti, where again complete softening occurred without change in 
orientation. (In this case, the identification of the process as subgrain 
growth is somewhat less definite, since the texture consists of two compo- 
nents.) McGeary and Lustman‘*®) obtained data that may be interpreted 
as indicative of 50°% loss of work hardening in cold-rolled polycrystalline 
Zr on annealing at 400°c, essentially without reorientation. Although 
all these results appear highly significant, the definition of orientation in 
these polycrystalline textures is not very sharp and, therefore, the detection 
of reorientation is not very sensitive. An instance of complete softening 
without reorientation in a deformed single crystal of silicon-iron was 
mentioned by Dunn‘). 
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Vickers hardness (D.P.H.N.) and residual x-ray line broadening (B) for pure Al 
single crystal rolled 80% on (110) plane in [112] direction and annealed 
at 350°c for periods indicated. [Ref. (30).] 


In a recent paper, Lutts and Beck presented data for pure Al single 
crystals rolled 80°% on a (110) plane and in a [112] direction, where a very 
sharp deformation texture results with the same orientation as that of the 
original crystal‘*. On isothermal annealing at 350°c it was found that 
the line broadening recovers completely at an early stage of the annealing 
when approximately 75°, of the work hardening is still present (fig. 24). 
This indicates that at least the considerable remaining fraction of the work 
hardening is not associated with other types of lattice imperfections 
detectable by line broadening measurements (such as strain or stacking 
faults) and should be presumably attributed to the presence of sub- 
boundaries. Furthermore, it was found that after annealing to the 
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completely soft condition at 350°c the material has exactly the same very 
sharp single orientation texture as before annealing, so that, in this case, 
softening is definitely not associated with reorientation, that is the 
migration of ordinary high angle grain boundaries (fig. 25, Plate 31). The 
elimination without recrystallization of work hardening presumably due to 
subboundaries should be attributable to subgrain growth, although on this 
crucial point direct evidence is still missing. 

The principal conclusions with regard to subgrains on the basis of results 
to date may be summarized as follows : 

(a) Subgrains, that is crystal domains slightly disoriented with regard 
to each other, are formed in metallic crystals on plastic deformation, 
presumably mainly by kinking, even at temperatures as low as 1/10 of the 
absolute melting point. The only deformation so far known not to be 
connected with subgrain formation is pure shear (without bending) in 
single crystals. 

(6) The subgrain size formed under otherwise similar conditions 
decreases with decreasing deformation temperature, increasing strain rate, 
increasing total strain (within certain limits) and with the addition of 
solute atoms. 

(c) Annealing under suitable conditions may produce subgrain growth, 
leading to a decrease of total subboundary energy per unit volume, and 
thus to increased thermodynamic stability. Subgrain growth is often 
prevented, particularly in deformed polycrystalline metals, by the pre- 
mature onset of recrystallization, which wipes out the substructure. 

(d) The rate of subgrain growth is greatly increased by an external 
stress applied during annealing, sufficient to produce a slight strain. An 
analogous effect in ordinary grain growth was not detected. 

(e) Large work hardening effects are connected with the substructures 
produced in plastic deformation. The fraction of the total work hardening 
due to this effect increases with the deformation temperature for any one 
metal. Subgrain growth on annealing may account for at least a part of 
the relief of work hardening. The extent of relief by this mechanism 
depends mainly on its rate, in relation to that of recrystallization. Some 
instances of complete softening without recrystallization are now known, 
for which subgrain growth is presumably responsible. 


$4, RECRYSTALLIZATION 
4.1. Primary Recrystallization, Gradual Grain Growth and Coarsening 


The classical picture of recrystallization, as the term has been quite 
generally used in the English language literature‘), involves the growth of 
new, strain-free grains at the expense of the cold worked matrix. On this 
picture, the following changes take place simultaneously, as a result of the 
growth of the new grains: (a) release of stored energy of cold work in pro- 
portion to the volume of cold worked matrix replaced by new, unstrained 
grains ; (b) softening, as work-hardened material is replaced by soft annealed 
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grains. In the partially recrystallized material hard and soft areas. 
coexist side by side, as may be ascertained by means of microhardness 
measurements(®) ; (c) local reorientation, as an ordinary high angle grain 
boundary migrates and a deformed grain of a certain orientation is 
replaced by an annealed grain of a different orientation (fig. 26, Plate 30). 
However, recent investigations made it very clear that these changes do 
not necessarily occur simultaneously'** 89. As discussed in § 3.5, instances 
are now known, where complete softening takes place prior to reorient- 
ation, that is essentially without the migration of high angle grain 
boundaries ; recrystallization may then subsequently occur without 
further softening. There have been also reported recently some instances. 
where a large part, or substantially all, of the stored energy of cold work is 
released prior to softening and to reorientation, as will be discussed in 
further detail below. In view of the variable and sometimes negligible, 
amount of residual work hardening and strain energy left in the matrix 
when it is being absorbed by the new grains, the phenomena of recrystal- 
lization in the above defined narrower sense (also called ‘ primary 
recrystallization ’), can not be sharply separated from the types of 
boundary migration (‘secondary recrystallization * or ‘ grain growth ’) 
‘where work hardening or strain energy admittedly play no significant role. 
Because of the lack of sharp distinction between these processes, it is. 
perhaps best to adopt the nomenclature used by Alterthum'®”, van 
Arkel'®§) and Burgers‘, where all boundary migration phenomena 
connected with the movement of ordinary high angle grain boundaries are 
designated by the generic term recrystallization. In the following, the 
term recrystallization will be used in this general sense. It has been 
customary to distinguish, largely on the basis of kinetics, between two. 
types of grain growth phenomena. In‘ gradual grain growth ’ or ‘ normal 
grain growth ’ many grains in a soft annealed matrix grow at the expense 
of others, resulting in an increase of the average grain size, without any 
substantial change in the grain size distribution. On the other hand, in 
‘abnormal grain growth ’, ‘ secondary recrystallization ’, ‘ discontinuous. 
grain growth’ or ‘ coarsening ’ relatively few grains which become very 
large, absorb a soft annealed matrix, which itself remains fine grained until 
absorbed. Characteristic of this process is the marked grain size contrast 
_ between the coarse grains and the matrix (‘ duplex structure’). In the 
following, the terms ‘ primary recrystallization ’, ‘ gradual grain growth ’ 
and “coarsening * will be used to designate the three rather roughly 
outlined types of recrystallization phenomena mentioned above. 


4.2. Driving Energy and Energy Release 
(a) Primary Recrystallization 
As stated in the Introduction, the energy absorbed in a metal specimen 
during cold working is of the order of 1-10% of the energy expended to. 
bring about plastic deformation: 2.3.4.5). Tt has been quite generally 
held’*® that most or all of this stored energy of cold work is released during’ 
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primary recrystallization, and that it represents the driving energy for this 
process. However, the recent investigation of Borelius and his asso- 
ciates®® showed that at least a portion of the stored energy may be 
released at relatively low temperatures and short annealing times, where 
recrystallization does not occur. Suzuki reported) that in 99-96°% pure 
copper deformed in compression to various amounts and heated at the 
rate of 2°c per minute the energy release is practically completed at a 
temperature (270°C), where the softening, as measured by indentation 
hardness or compression yield stress, just begins. It may be assumed 
that in this case the energy release is associated mainly with processes 
occurring prior to grain boundary migration. Preliminary results by 
Stansbury? with a silver bearing copper appear to confirm Suzuki’s 
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Isothermal release of stored energy of cold work at 181°c from pure Cu extended 
15%. (Unpublished data by P. Gordon.) 


findings. On the other hand, Clarebrough, Hargreaves, Michell and 
West) found with phosphorized coppers of various purity that on heating 
up at 2° per min, by far the largest part of the energy release took place 
concurrently with softening and with microscopically observable grain 
boundary migration. In an as yet unpublished investigation, Gordon'?) 
observed that in high purity copper deformed in tension the largest part 
of the energy release was also associated with recrystallization. Gordon's 
isothermal energy release curves at 200°c (figs. 27) are interesting in that 
they show an initial energy release with typically recovery type kinetics, 
similar to that found by Borelius and his associates®®, followed by a 
clearly separated large peak of energy release corresponding to recrystal- 
lization, and with kinetics typical of the latter process. The cause of the 
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discrepancy between the various investigators’ results has not as yet been 
clarified. ‘The conditions of deformation, the impurities in the copper, 
and the method of measurement were different in each case. 

In contrast to their results with copper, Clarebrough, Hargreaves and 
West“) find that deformed nickel releases about 60%, of the stored 
energy of cold work prior to the onset of recrystallization. This suggests 
that the unexplained differences in the results for copper, mentioned 
above, may be due to real differences in the behaviour of the various 
specimens rather than to the methods of measurement. At the present 
stage of this work it appears safe to conclude only that in some cases 
most of the energy release is associated with the absorbtion of the deformed 
matrix by new grains, as in the classical picture of recrystallization, but 
that this is not true in other cases. Under certain conditions, new grains 
may grow in a cold worked matrix, which has already lost the largest 
part of the stored energy of cold work in the early phases of annealing 
without grain boundary migration. The remaining minor fraction of 
the stored energy of cold work, which in such cases serves as the driving 
energy of the recrystallization process, is probably associated with the 
network of subboundaries in the matrix, as suggested by C. S. Smith), 
In the experiments of Lutts and Beck, referred to in the section on Sub- 
erain Growth, ‘ recrystallization ’ i.e. the migration of high angle boundaries 
takes place at a stage of annealing after the completion of softening 
(figs. 24 and 25), where not only recovery had been completed, but 
presumably even some of the subboundaries had been annealed out by 
means of subgrain growth. It is very likely that in this case, too, most 
of the stored energy of cold work was released prior to recrystallization. 


(b) Grain Growth and Coarsening 


The early suggestion of Ewing and Rosenhein™), Tammann(‘104), 
Vogel) and Alterthum'” that the driving energy of grain growth in 
soft annealed metals is the interface energy associated with the grain 
boundaries has received strong support in recent years through the work 
of Harker and Parker ®, C. 8S. Smith?) and J. E. Burke), Burke 
and Schiau” have shown that small residual strains, such as might 
remain in an annealed metal, have little if any effect on gradual grain 
growth. As estimated from the value of the interfacial energy per unit 
interface area and from the change in total interface area per unit volume, 
the energy released when a metal with an average grain diameter of 
10 cm undergoes grain growth to an average grain diameter ten times 
larger is of the order of magnitude of 0-2 cal/mole. This is in the order 
of 100 times smaller than the stored energy of cold work, which is released 
prior to or during primary recrystallization. The driving energy of grain 
growth or coarsening is, of course, even smaller if the initial grain size 
is, as commonly happens, larger than that assumed above. In view of 
the small magnitude of this energy, it is understandable that no successful 
attempts to measure it directly are recorded in the literature. . 
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4.3. Kinetics 
(a) Gradual Grain Growth (Normal Grain Growth) 

Vogel) showed a tendency of grain boundaries to migrate in such a way 
as to decrease their curvature, and ascribed this tendency to the require- 
ment of minimizing the interface energy associated with the boundary. 
Harker and Parker“®® in their more complete analysis later found that in 
addition to the tendency of grain boundaries to straighten themselves, 
approach to equilibrium requires also the migration of grain boundaries so 
as to minimize the resultant surface tension forces exerted at a grain edge 
by the three grain boundaries intersecting at it. They recognized that 
grain boundary migration is brought about by annealing whenever the 
system of grain boundaries in a metal does not satisfy both equilibrium 
conditions mentioned. Smith” showed that true equilibrium arrange- 
ment of grain boundaries in three-dimensional space would be realized 
only if each grain had the shape of a minimum-area tetrakaidecahedron, 
in complete analogy to the stable arrangement of soap films in a froth. 
He, and later in more detail Burke‘!°8), pointed out that a deviation from 
true boundary equilibrium configuration at any point in a polycrystalline 
metal must result in continued boundary migration in the whole specimen 
because of the tendency of the deviation to spread to other areas. Since 
the deviations from equilibrium are not eliminated from the metal despite 
the boundary migration, true equilibrium is in practice not achieved 
until the whole specimen is converted into a single crystal. As described 
by Smith, a curved grain boundary will gradually migrate on annealing 
toward its center of curvature, requiring further continual readjustments 
at the edges. This process takes place relatively slowly, the rate in a 
specific case depending on the local curvature of the boundary and on the 
temperature. At the moment this adjustment process leads to a locally 
highly unstable configuration where more than three surfaces meet at an 
edge, or more than four edges meet at a point, the grain boundaries rapidly 
rearrange themselves, approaching new positions of local equilibrium. 
As confirmed by direct microscopic observation of grain boundary migra- 
tion in pure Zn slightly above room temperature, locally the process 
consists in alternating slow migration and quick readjustment. Small, 
tetrahedral grains have convex boundaries. These necessarily move 
toward the center of such a grain, until it disappears. As the number of 
grains per unit volume decreases the average grain size increases. 

In pure metals or solid solutions with a small initial grain size, annealing 
leads to an increase of the average grain size, but it leaves the type of the 
frequency distribution of the grain sizes more or less unchanged. Hu and 
Beck 125) made a quantitative study of this phenomenon for isothermal 
grain growth in pure aluminium. Planimetric measurement of the area 
of each of approximately 500-1000 grains in the center section of strip 
specimens at three successive stages of grain growth gave grain size 
distribution curves that are similar to each other (fig. 29), and can be 
made to coincide fairly well by proportional changes in the grain size 
scale (fig. 30). As seen above, the increase in the average grain size 1s 
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accomplished by the disappearance of the smallest tetrahedral grains. 
At any point during the annealing some larger grains with more than four 
‘boundary surfaces are being partially absorbed by their neighbours until 
they lose some of their boundary surfaces, eventually become tetrahedral, 
and finally disappear. The largest grains, i.e. those with the largest 
number of boundary surfaces are the ones that tend to grow most. 
During grain growth the number of grains of a given size range con- 
tinually decreases, since some grains leave the group by becoming 
too small or too large, while at the same time this number is increased 
the grains just becoming small enough or large enough to fall within 
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The relative frequency of grains with various numbers of sides in plane section 
remains essentially unchanged during grain growth in pure Al. Rolled 
Al specimens annealed at 500°c for periods shown (see fig. 32 and 


Ref. (113), Mmicrographic tracings by H. Hu, unpublished data on 
number of sides by P. R. Sperry and N. Gottesmann). 


the size range considered. As seen in fig. 29, the net result is that 
the frequency of the largest grains at any time increases, and that 
of the smallest grains decreases. The frequency of grains of certain 
intermediate sizes may first increase and then decrease in the course of 
continuing grain growth. Although the local boundary movements are 
quite discontinuous, the overall change with annealing time of the 
structure as a whole, i.e. of the statistics of the grain size distribution 
is a very gradual or continuous one. The frequency of grains with various 
numbers of sides (in a plane section through the specimen) remains 
practically unchanged during gradual grain growth (fig. 28). 
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The increase of the average grain size with time on isothermal annealing 
follows a parabolic relationship and it can be represented empirically1™ 
by the following formula : 


Din _D"=C . t, te 8) Bee (2) 


where ¢ is the annealing period for grain growth, D, is the average grain 
diameter before and D is that after annealing, C and n are parameters, 
which depend on the temperature, but are independent of D. Figures 
31 and 32 give experimental data for high purity brass and aluminium. 
For a given temperature the rate of grain growth depends only on the 
instantaneous grain size4) (fig. 33). From (2) the instantaneous rate 
of grain growth 

aan. C.D, Site: week Bat ie) 
The value of » for pure brass is approximately 0-5, independently of the 
temperature. For pure aluminium n varies from 0-056 at 350°c to 
0-322 at 600°C, extrapolating to a value of near 0-5 at the melting point. 
Since n is always less than one, (3) gives rates of grain growth decreasing 
with increasing grain size (figs. 34 and 35). 

As seen in fig. 36, if log D is plotted against the log of the annealing 
time for cold rolled strips of brass, the initial points often deviate from a 
straight line relationship. This is to be expected from (2)‘11”, since the 
total isothermal annealing time includes a period for recrystallization 
and this, in general, differs from the hypothetical annealing time that 
would be required to produce by grain growth from an initial grain size 
of zero, that grain size which is in fact produced by primary recrystalliza- 
tion. Depending on the amount of deformation, and thus on the time 
necessary for recrystallization, this initial deviation may be either upward 
or downward (fig. 36). However, if the time scale is shifted so as to 
plot the log of the time for grain growth alone, but include in this also the 
time for the hypothetical grain growth from zero size to the recrystallized 
grain size, then in accordance with (2), all points fall on a straight line 
(fig. 37). The results here reported were obtained by Towers and Beck‘?”), 
who carried out the necessary correction in the time scale for each curve 
by means of independent determinations of the time for complete 
recrystallization and of the grain size as recrystallized. 

A relationship, such as (2) with n=0-5, may be derived") if one 
assumes that the instantaneous rate of grain growth is proportional to 
the instantaneous value of the excess free energy per unit volume associated 
with the grain boundaries. Since this energy may be taken as propor- 
tional to the grain boundary surface area per unit volume, it follows : 


— =A, Se Ce ee 
Th ee ) 
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GRAIN DIAMETER. D, mm. 

Relative frequency of grains of various diameters in plane section, after 
grain growth at 500°c for 1 minute (solid circles), 25 minutes (half-filled 
circles) and 625 minutes (empty circles). Same pure Al specimens as 
in fig. 28. See also fig. 32. (Unpublished data from 500-1000 grains 
in each specimen by H. Hu and P. A. Beck.) 
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ADJUSTED GRAIN DIAMETER, D, mm. 


The relative frequency data of fig. 29, replotted over grain diameter adjusted | 


proportionally to make the average grain diameter in each specimen 
0:96mm. The type of grain size frequency distribution remains 


essentially unchanged during grain growth. ‘(Unpublished data b: 
H. Hu and P. A. Beck.) 8 (Unpublished data by 
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Fig. 31 
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Isothermal grain growth in high purity 70-30 brass. (Unpublished data by 
J. Towers and P. A. Beck.) 
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Rate of isothermal grain growth in high purity aluminium at 450°c depends 
only on instantaneous grain size, regardless of whether specimens were 
previously annealed at 450°c (Nos. 1890 and 1896) or at 400°c (Nos. 
1882 and 1886). Data for 400° and 450°c are replotted from fig. 32 for 
comparison, showing high degree of reproducibility. (Unpublished work 
by M. L. Holzworth and P. A. Beck.) 
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Average instantaneous rate of grain growth for high purity ( —) and for 


commercial (— — — -) 70-30 brass, as a function of grain size. [Un- 
published data by J. Towers and P. A. Beck and from Ref. (114).] 
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Burke noted) that (4) also follows from the assumption that the 
instantaneous rate of grain growth is proportional to the average grain 
boundary curvature at the instant. More recently, Machlin‘!15) showed 
that the assumption underlying the derivation given above may be 
justified on the basis of the thermodynamic theory of irreversible processes. 
Fullman and Fischer!” obtained experimental evidence for the validity 
of (5) for bubble growth in soap froth. As stated above, the value of 
n=0-5 is also approached in pure brass over a wide temperature range, 
and in pure aluminium near the melting point. A discussion of the various 
factors affecting the value of n was reviewed by Fullman‘!!”), 
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Average instantaneous rate of grain growth for high purity aluminium, as a 
function of grain size. [Ref. (114).] 


It was observed‘119) that grain growth on annealing of a rolled poly- 
crystalline metal strip becomes slower than expected, and eventually 
stops altogether, when the grain size approaches the thickness of the 
strip (fig. 38). The deviation from the straight line relationship at large 
grain sizes in fig. 32 results from this effect. This ‘ specimen thickness 
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effect ’ is also observed in wires, where grain growth stops when the grain 
boundaries become approximately perpendicular to the wire axis, and 
the grains occupy the whole cross section of the wire, having an average 
length of the order of the wire diameter. ; 

The grain growth behaviour discussed above is found only in pure 
metals or solid solutions.* When finely dispersed particles of a second 
phase are present in the structure, the grain boundaries tend to cling to 
these particles‘1°7, 118), fig. 39 (Plate 32), so that grain boundary 
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Average grain diameter as a function of total annealing time for high purity 
70-30 brass with large penultimate grain size and final rolling of 17:-4% 
(filled circles) with large penultimate grain size and final rolling of 
33°3% (crosses) and with small penultimate grain size and final rolling 
of 85-4°% (empty circles). (Unpublished data by J. Towers and P. A. 
Beck.) 


migration is inhibited. In such cases, the rate of grain growth becomes 
much smaller than without the inhibiting dispersion, as may be seen for 
instance by comparing the data for commercial brass with those for pure 
brass (fig. 34). Also the exponent , ie. the slope of the logarithmic 


* The solid solutions most extensively studied are alpha brass, as described 
above, and Al-Mg of various compositions up to 2°% Mg‘), Both of these 
exhibit gradual grain growth, as that in pure metals. The view has been 
expressed("1® that solid solution of Ag in Zn inhibits grain growth, but this 
was based on inadequate experimental evidence, and may be considered unlikely. 
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Same grain sizes as in fig. 36, replotted as a function of time for grain growth 
(i.e. total annealing time /ess time for primary recrystallization plus time 
for imaginary grain growth to produce grain size actually attained by 
primary recrystallization). (Unpublished data by J. Towers and 
P. A. Beck.) 


Fig. 38 


MAXIMUM GRAIN SIZE,MM. 


I 2 
SPECIMEN THICKNESS,MM. 


Maximum grain size attainable by gradual grain growth in high purity aluminium, 
as a function of specimen thickness, showing ‘specimen thickness 
effect’. (a) Specimens not etched before annealing. (>) Specimens 
etched before annealing. Empty circle: specimen extremely deeply 
etched before annealing. [Ref. (113).] 
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grain growth lines often decreases ; for commercial brass n=02, as 
compared with n=0-4 for high purity brass. Still another important 
effect of inhibition by a dispersed second phase Is ete grain eee 
stops completely when a certain ‘ultimate grain size" 1s reached! : 
As the grain size increases during gradual grain growth, the driving energy 
for further grain growth becomes smaller until, at a certain point, the 


Fig. 40 
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Ultimate grain size as a function of final annealing temperature in Al alloys 
with 06% (4), 0:-75% (3/4), 1:12% (1) and 1-65% (14) Mn, after 
penultimate anneal at 530°c (A), penultimate anneal at 500°c (B) and 
penultimate anneal at 600°c (D). For the same final annealing tempera- 
ture, the ultimate grain size decreases with increasing number of 
inhibiting particles (i.e. with increasing Mn-content and with decreasing 
penultimate annealing temperature). [Ref. (118).] 


inhibiting effect of the dispersion is just sufficient to prevent further 
growth), The ultimate grain size so reached is stable over long 
periods of annealing. As seen in fig. 40, the ultimate grain size decreases 
with increasing number of dispersed particles. It increases with increasing 
annealing temperature, although this effect, as observed experimen- 
tally“), may be at least partly due to a decrease in the number of 
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dispersed particles through their re-solution and coalescence. It is note- 
worthy that even very few dispersed particles have observable effects on 
gradual grain growth when the grain size becomes large enough, so that 
gram growth may be considered as one of the most sensitive indicators 
of such particles. Inhibition effects have been observed even in high 
purity brass. Figure 41 shows the somewhat gradual change-over in 
high purity brass from grain growth limitation by specimen thickness to 
that by inhibition. 
Fig. 41 
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Ultimate grain size in high purity (Brass A) and in commercial (Burke) 70-30 
brass, as a function of specimen thickness, for the annealing temperatures 
indicated. When the thickness is small it limits gradual grain growth, 
independently of temperature and of purity. For larger thicknesses grain 
growth is limited by inhibition, so that the ultimate grain size depends 
on both temperature and impurities. [Unpublished figure, by J. 
Towers and P. A. Beck, data for commercial brass from Ref. (108).] 


(6) Primary Recrystallization 

Studies by Mehl and his associates'1!®, 12° of the volume percentage 
of material recrystallized on annealing of slightly extended fine grained 
polycrystalline specimens as a function of annealing time indicate that. 
the kinetics of the process is typically represented by a ‘ sigmoidal ’ 
curve, as shown in fig. 42. This kinetics is in sharp contrast to that 
typical for recovery, fig. 6. Johnson and Mehl have shown?) that the 
‘sigmoidal ’ kinetics, where the reaction rate is at first very low, then 
increases to a maximum and again decreases, is characteristic of trans- 
formations that may be formally described in terms of two simultaneous 
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processes, namely the nucleation of new grains and their growth. Both 
primary recrystallization and coarsening (secondary recrystallization) 
may be described in that manner. Goeler and Sachs‘») first carried out 
a quantitative analysis for primary recrystallization, assuming rates of 
nucleation and of growth independent of time, and making no allowance 
for the decrease of transformation rate due to impingement of the 
recrystallized grains. In the more recent and more complete studies 
by Stanley and Mehl!) and by Anderson and Mehl‘? the rate of 
growth @ was still assumed to be constant and the same for all new 
grains, but no assumption was made regarding the rate of nucleation, J, 
which was determined experimentally. In these studies, G has been 
determined by measuring the diameter of the largest grains in a series 
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Isothermal recrystallization curve at-350°c for pure Al extended 5-1 %. 
[Ref. (120).] 
of specimens annealed for various periods at constant temperature. The 
results fit well the interpretation, according to which the first new grain 
in each of these specimens formed after the same ‘ period of incubation ’ 
and then grew at the same rate, so that this rate could be determined 
from the slope of lines like those in fig. 43, without additional assumptions. 
After plotting the total number of recrystallized grains as a function of | 
annealing time and graphical differentiation, the rate of increase of the 
number of recrystallized grains per unit of untransformed volume, i.e. 
the rate of nucleation N, could be calculated as a function of time, since 
the volume fraction transformed was measured as a function of annealing 
time. ‘Typical results are shown in fig. 44. Anderson and Mehl deter- 
mined both the rate of growth and the rate of nucleation in pure Al as 
a function of deformation and of annealing temperature, figs. 45 and 46. 
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Radius of largest grains vs. annealing time at 350°c of pure Al specimens, after 
extensions indicated. [Ref. (120).] 


Fig. 44 
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Rate of nucleation as a function of annealing time at 350°c, after 5% 
extension. |[Ref. (120).] 
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From the temperature dependence of these rates, activation energies 
were calculated for both N and @ in the range of deformations studied, 
from about 1% to 15% extension (fig. 47). A proof of the correctness 
of the interpretation of the process in terms of nucleation and growth 
and of the assumptions used was obtained by Anderson and Mehl by 
calculating the amount of material transformed as a function of time 


Fig. 45 
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Rate of growth at temperatures indicated vs. per cent extension of pure Al 
specimens. [Ref. (120).] 


from the independently determined values of G and N, with the help of 
eqn. (6) of Johnson and Mehl”2) for the three dimensional case 
and of eqn. (7) derived by Stanley and Mehl?!) for the two dimensional 


case 
Pcie 8G8a/,, G8 bi 
fHj)=1 exp. bf (c ee “§. “See a9) —bi— 1) err Ss (6). 
2G2%a fet bt? 1 
ha) exn= Aare 
f()=1—exp— (5 : -!-;) (7). 


where f(t) is the volume fraction recrystallized at time ¢ and a and 6 are- 
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the empirical constants giving the rate of nucleation as a function of 
time according to eqn. (8) : 

VOX) Oh oat eee i Se aeats d+ | (8) 
Equation (8) approximates the ascending portion of the curve showing 
the experimentally determined data in fig. 44. Equations (6) and (7) 
were derived on the basis of the nucleation and growth picture and of a 
suitable assumption to take into account the effect on the transformed 
volume of the impingement of the recrystallized grains on each other. 


Fig. 46 
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Rate of nucleation at temperatures indicated at 20% recrystallization vs. per 
cent extension of pure Al specimens. [Ref. (120).] 


The curves so calculated were in good agreement with the values of the 
volume fraction transformed, directly measured as a function of annealing 
time. 

In arecent paper, H. Suzuki") indicated that in high purity aluminium 
the rate of growth G of recrystallized grains after smal] extensions decreases 
with gradual grain growth in the deformed matrix material. This obser- 
vation was interpreted by Suzuki in accordance with the view expressed 
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earlier") that the grain boundary interface energy may contribute a 
substantial part of the total driving energy for primary recrystallization 
in a polycrystalline metal, when the matrix grain size and the deformation 
are very small. In addition, the residual stored energy of cold work 
itself may be gradually eliminated by continuous grain growth "** Oe 
the portions of the matrix swept by moving boundaries become essentially 
free of strain and of substructure. For both of these reasons, the 
driving energy for the growth of large recrystallized grains may be 
expected to decrease as a result of the concurrent gradual grain growth 
in the matrix. Suzuki’s results for low deformations, as shown in fig. 48, 
differ radically from those of Anderson and Mehl, fig. 43. The cause of 
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Activation energies for nucleation and growth in pure Al vs. per cent extension. 
[Ref. (120).] 


the discrepancy may be the higher purity of the aluminium used by 
Suzuki. The possibility that the metal employed in the work of Anderson 
and Mehl may have contained dispersed second phase particles in sufficient 
quantity to inhibit continuous grain growth in the matrix is suggested 
also by the small matrix grain sizes reported by these investigators. 

With increasing deformation, the number of recrystallized grains 
becomes larger, that is the grain size smaller‘) (fig. 49), until the new 
grains can no longer be separated from the original deformed grains by 
means of simple criteria, such as size. Also, a considerable fraction of 
the total growth of each new grain may take place within a single deformed 
matrix grain, and the rates of nucleation and growth vary greatly with 
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Fig. 48 
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z ‘Radius of grains artificially nucleated vs. annealing time in minutes at 506°C, 
a: after 1-25°% extension. [Ref. (123).] 
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the orientation of the individual matrix grains with respect to the stress 
axes, Because of this variation and of that resulting from the nhomo- 
geneity of the plastic deformation due to grain boundaries and kink 
bands, the local differences in the rates of nucleation and growth are very 
great. For the various reasons mentioned a detailed analysis of nuclea- 
tion and growth rates, such as that accomplished by Mehl and_ his 
associates for primary recrystallization after small deformations (where 
the described local variations are averaged out due to the large size of the 
recrystallized grains), would be difficult or impossible to carry out for 
metals more severely deformed. However, the average overall kinetics, 
that is the volume fraction of material recrystallized as a function of 
annealing time, can be determined accurately even after severe defor- 
mation in cases where recrystallization is associated with a pronounced 


Fig. 50 
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Isothermal recrystallization curves for pure Cu at temperatures indicated after 
98% cold rolling. [Replotted from Ref. (126).] 


reorientation of the texture. Decker and Harker‘2® were thus able to 
measure quantitatively the kinetics of reorientation in 98°% rolled copper 
where the recrystallized material has a very sharp ‘cube texture 
Figure 50 shows that the kinetics is typically ‘ sigmoidal’ like that found 
after small deformations (fig. 42). Similar results were obtained by 
Seymour and Harker?” for a highly rolled Ni-Fe alloy, and by McGeary 
and Lustman'*®) for zirconium. These results suggest that primary 
recrystallization even after severe cold working may be described as a 
nucleation and growth process, a conclusion which is, of course, in 
excellent agreement with numerous qualitative microscopic observations: 
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(c) Coarsening (Secondary Recrystallization) 


As mentioned above, grain growth in soft annealed metals may take 
place by a gradual or ‘ continuous ’ mechanism, which leaves the type 
of grain size frequency distribution relatively unchanged. However, if 
conditions are such that most grain boundaries in the structure are not 
able to migrate, although there are a few which are free to do so, the 
kinetics of the process of grain growth become quite different from those 
discussed in (a). The relatively few growing grains become very much 
larger than the surrounding matrix grains (fig. 51, Plate 33), so that the 
type of grain size frequency distribution is also radically altered. The 
splitting of the structure into two clearly separated portions, namely the 
new, coarse grains and the remaining relatively unaltered matrix grains, 
gave rise to the designation ‘ duplex structure ’. There are now known(‘!29) 
two fundamentally different causes that may be responsible for the 
immobility of grain boundaries and, hereby, for the type of grain growth 
discussed.* One of these is inhibition by dispersed particles of a second 
phase, as described above. The other one is connected with the orienta- 
tion dependence of boundary mobility, to be discussed in §4.4; if the 
orientation of almost all grains in the structure is nearly the same, most 
grain boundaries have very low mobility, and only the few grains with 
deviating orientations can grow. These two causes of boundary im- 
mobility may be operative independently from each other, but they are 
in many technological situations both active at the same time. For this 
reason, it is not always easy to tell whether a given case of coarsening is 
of the inhibition-dependent or of the texture-dependent type. Also in 
some cases the dispersed particles are so small that they can not be 
detected microscopically. For instance, the well known coarsening 
phenomena in Al-killed low carbon steels are typical of inhibition- 
dependent coarsening, yet for many years it has not been possible to 
identify the dispersed particles responsible for them. Recent investiga- 
tions indicate that this dispersed phase is connected with the nitrogen 
content, so that it is most likely aluminium nitride“?°). 

A detailed investigation®) of grain growth in Al-Mn alloys containing 
a dispersed second phase showed that coarsening occurs when a small 
initial grain size is stabilized by precipitating a large amount of inhibiting 
second phase and the structure is then annealed for prolonged periods 
of time at a temperature just below the solvus point (fig. 52). The latter 
treatment causes the precipitate to redissolve partially and to coalesce, 
so that the number of inhibiting particles is slowly reduced to a lower 
value. In the course of annealing, the inhibition initially sufficient to 
immobilize the grain boundaries at a small grain size level is gradually 
decreased, until at certain localized areas, where conditions happen to be 
particularly favourable, grain growth may start, while it is still inhibited 
OE ee, ee eee ee Se 


* The third known cause of grain growth stoppage, namely the ‘ specimen 


thickness effect ’, is not known to give rise to coarsening. 


P.M. SUPPL.—JULY 1954 x 


292, Paul A. Beck on the 


elsewhere. Once a grain becomes considerably larger than its neighbours 
it can grow at their expense at a level of inhibition still sufficient to 
prevent the migration of boundaries between the small grains. This 
condition obtains, since at the edges of the large grains the surface tension 
forces of the intersecting boundaries are further removed from equilibrium, 
and also the boundaries of the larger grains have greater curvatures 
than those between the small matrix grains‘. It is significant that, if 
the inhibition is removed too fast at too many points, for instance by 
annealing at a temperature slightly above the solvus point, the grain 
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Isothermal grain growth at final annealing temperatures indicates, in Al+1-1% 
Mn alloy vs. annealing time. Vertical dash-dot line indicates approxi- 
mate start of coarsening. [Ref. (118).] 
growth then taking place is of the gradual rather than of the discontinuous 
type. Since gradual grain growth is limited by the specimen thickness, 
while the grain size attained in coarsening is not, it is possible under 
suitable conditions to arrive at much larger grain sizes by annealing at 
a somewhat lower temperature than those developed at a higher tempera- 
ture (see fig. 52), 
Texture-dependent coarsening occurs only when the matrix has a 
well developed single orientation texture. If the matrix texture has more 
than one component, no matter how sharply each is developed, in general 
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the structure contains many high angle boundaries and these are free to 
migrate, so that the resulting grain growth is of the gradual type. The 
only examples of pure texture-dependent coarsening now known are 
those occurring in a cube texture matrix, or in a matrix obtained by 
primary recrystallization from a homogeneously deformed single crystal. 

The kinetics of coarsening has not yet been investigated in as much 
detail as that of gradual growth and of primary recrystallization. The 
early investigations of W. Feitknecht(3, Karnop and Sachs) and 
others established the general features of the phenomenon, and showed 
that its kinetics may be similar in many respects to a nucleation and 
growth type process. This was later demonstrated particularly clearly 

Fig. 53 
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Diameter of largest grain formed by coarsening in silver vs. annealing time at 
temperatures indicated, after final rolling reduction of 50%. [Ref. (133).] 


by the investigations of Rosi, Alexander and Dube"). The analogy 
between the kinetics of primary and secondary recrystallization was. 
emphasized by Burke and Turnbull’. Although there are so far no 
quantitative data available on the volume transformed as a function of 
isothermal annealing time, in analogy to the relationships in primary 
recrystallization a ‘ sigmoidal * type transformation curve is certainly to 
be expected, at least when the specimen size is sufficiently large in com- | 
parison with the resulting grain size. 

The best published data on the rate of growth in coarsening are those 
of Rosi, Alexander and Dube for silver, shown in figs. 53, 54, and 55, 
As seen in fig. 55 these data define a thermal activation type temperature 
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dependence, with an activation energy of 28 kcal/mole. According to the 
calculations of Turnbull'!8" the free energy of activation for this boundary 
migration is 21-5 kcal/mole as compared with 18-6 kcal/mole for the free 
energy of activation for self diffusion along high energy grain boundaries. 
The reasonably good agreement here supports the interpretation?) 
according to which the elementary process determining boundary mobility 
is similar to that determining the rate of self diffusion along the boundary. 
Certain theoretical questions involved are discussed fully in the review 
by Burke and Turnbull™, which also gives a tabulation of available 
data for free energies and entropies of activation, calculated on the basis 
of the absolute rate theory. 

Although statistical rates of ‘nucleation’ have not been determined 
for coarsening, Rosi, Alexander and Dube measured incubation periods 
for individual secondary grains in silver as a function of temperature. 
They found a temperature dependence of the reciprocal incubation period 
very similar to that of the rate of growth in the same material, with an 
activation energy of 30 kcal/mole. From this, they concluded that the 
incubation period corresponds to ‘ invisible growth ’ at a lower rate when 
the secondary grain is not sufficiently larger than the neighbouring 
matrix grains and it has not yet emerged as a secondary grain. They 
found that in order to be able to grow even very slowly, a secondary grain 
had to be at least twice as large as its neighbours, and that its rate of 
growth increased with its size up to a grain diameter approximately 
13 times larger than the average matrix grain size. Beyond that size, 
the rate of growth became independent of the size of the secondary grain. 
Even though the material used in this investigation was high purity 
silver, the occurrence of distinct coarsening in a matrix containing 
several orientations suggests that second phase inhibition may have played 
arole. (In 70-30 brass, under analogous conditions of preferred orienta- 
tion, gradual grain growth instead of coarsening takes place.) For 
inhibition-dependent coarsening, the measured incubation periods are 
often influenced by such factors as the rate of re-solution of the dispersed. 
phase, and the extent of gradual grain growth in the matrix prior to 
coarsening. A rather clear indication that even small differences in 
matrix grain size may cause large changes in the incubation period in 
coarsening was obtained in the investigation with Al-Mn alloys‘?®. 
Complex relationships sometimes encountered are shown in fig. 56, from a. 
study by Sperry and Beck'1?® of coarsening in a high purity Al-alloy 
with 0-35°% Fe. When a high penultimate annealing temperature was. 
used, the incubation period showed a clear-cut minimum at a final 
annealing temperature of 600°c ; at higher final annealing temperatures 
the incubation period actually increased, so that the apparent ‘ activation 
energy ’ between 600°c and 650°C was negative. ‘This was not the case 
when the penultimate annealing was carried out below 606 °C, or at 
600°c for 5 minutes only. Apparently, the amount of Fe dissolved 
during the penultimate anneal and re-precipitated during cold working 
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had an important part in the overall inhibition effect. Presumably, the 
fast re-solution of this very fine precipitate on annealing in the tempera- 
ture range of 600°-650°c allowed some gradual grain growth to take place 
in the matrix and thus delayed the ‘ nucleation ’ of coarse grains. 


Fig. 56 
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Approximate annealing periods for beginning of coarsening in pure Al+0:34% 
Fe alloy at final annealing temperatures from 500° to 650°c, after the 
various penultimate annealing conditions stated and a final rolling 
reduction of 80%. (Unpublished research by P. R. Sperry and P. A. 
Beck.) 


4.4. Problems of Nucleation in Recrystallization 

As mentioned in the previous section, it is commonly recognized that 
primary recrystallization, as well as secondary recrystallization (coarsen- 
ing) may be at least formally described in terms of nucleation and growth, 
although this is not true in the case of gradual grain growth. Even 
though the growth of new grains after a period of incubation is readily 
observed by metallographic means in both primary recrystallization and 
in coarsening, the origin of the ‘nuclei’ of these new grains has been 
the subject of much speculation and it is still not definitely known. 


(a) Primary Recrystallization 


All theories recognize the condition, first formulated by R. Becker(37), 
that a strain-free nucleus is stable and capable of growth in a strained 
matrix if it is not smaller than a certain critical size. At the critical 
size, for a small amount of growth the free energy gained from the 
additional volume transformed is just sufficient to compensate for the 
free energy required to increase the interface area. The various theories 
differ mainly in regard to the mechanism by which nuclei of critical or 
larger size are formed and to the interpretation of the incubation period. 
In their recent review, Burke and Turnbull” gave a rather complete 
description of the various theories that have been proposed to explain 
nucleation in primary recrystallization. 
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In the * conventional nucleation theory ° strain-free nuclei are formed 
by means of thermal activation, at a rate increasing sharply with the 
local lattice strain. The incubation period in this theory is the time 
required at a given temperature for a sequence of thermal fluctuations 
until a sufficiently powerful activation event is encountered leading to 
the formation of a strain-free ‘nucleus’ of not smaller than critical size 
at a certain site in the deformed matrix. Burke and Turnbull estimated 
that a local elastic strain of the order of 0-2 or greater would be required 
in order that a certain site should effectively contribute to the formation 
of recrystallization nuclei by this mechanism. It is necessary to assume 
that such large strains are present in at least a small fraction (for example 
in the order of 10~®) of the total volume. If this volume fraction is small 
enough, the locally very high energy of these sites and the relatively 
high energy regions which it is necessary to assume in their immediate 
vicinity in order that the nuclei can grow, may constitute only a very 
small part of the residual stored energy of cold work in the specimen as 
a whole. Consequently, it is possible to postulate a strain distribution 
in the cold worked matrix that would allow nucleation by this mechanism, 
without requiring more average stored energy than the measured values. 
However, it is uncertain what fraction of the measured stored energy of 
cold work is actually retained at the point in the annealing process where 
recrystallization starts. As stated in §4.2, there are at least some 
instances where apparently most of the stored energy of cold work is 
released prior to the main part of recrystallization. It is, therefore, 
probable that this theory is not generally applicable. 

In the so-called * block theories ’ the nucleus is pictured as a pre-existing 
block of a size sufficient to grow, that has either remained during the 
deformation process relatively strain-free with respect to its surroundings 
(‘low energy block theory ’), or has become relatively strain-free during 
annealing by means of some process that decreases the strain uniformly 
in the whole block (’ high energy block theory’). In contrast to the 
‘ conventional nucleation theory ’, in the ‘ high energy block theory ’ the 
nucleus-block may be quite large, since the underlying assumption is 
that the nucleus is not formed by a single activation event, but rather 
by a gradual process requiring a large number of thermal activation events 
within the block volume, whose effect is cumulative. 

In the low energy block theory, the incubation period is the time for 
initial slow growth before the higher, final growth rate is attained. In 
analogy to the findings of Rosi, Alexander and Dube) for coarsening, 
it has been recognized that in primary recrystallization, too, the rate of 
growth of a nucleus of the critical size must be practically zero, and that 
this rate may be expected to increase gradually with increasing size, as 
the gain of free energy for a small transformed volume becomes consi- 
derably larger than the loss due to the corresponding increase of interface 
area. Recently, Cottrell18%) proposed an additional reason for low 
initial rates of growth, and for increased rates as the new grains become 
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larger. As discussed in the following section, there is now good evidence 
indicating that the mobility of a grain boundary depends on the orientation 
difference between the two grains on either side of it, having very low 
values for small disorientations. Cottrell suggests that originally the 
nucleus may be only slightly disoriented with regard to its environment, 
but that as it grows, the boundaries reach nearby matrix regions (e.g. 
kink bands) where the orientation is quite different. From there on the 
rate of growth may become much larger. In addition to giving a plausible 
interpretation of the incubation period, this picture also has the advantage 
of permitting a considerably smaller size for a stable nucleus, in terms of 
the Becker condition. than the assumption of a nucleus ‘of strongly 
deviating orientation, the radius of a stable spherical nucleus is propor- 
tional to the interface energy”), which in turn increases greatly with the 
disorientation''®®), 

The theory of nucleation proposed by Cahn‘14°), assumes that essentially 
strain-free blocks of very small size form by polygonization in the early 
phase of annealing. It is implied, however, that the incubation period 
is not determined by this process, but rather by the rate at which some of 
the resulting subgrains grow sufficiently large at the expense of neigh- 
bouring subgrains of nearly the same orientation to serve as nuclei 
capable of growing at the expense of matrix material of strongly deviating 
orientation. In this theory, which in part leans on the earlier suggestion 
by Burgers“4), the incubation period corresponds to the annealing 
time necessary for the initial subgrain growth. 

Since thus polygonization itself is not an essential part of the theory, 
the usefulness of Cahn’s idea is not affected by the recent experimental 
confirmation of the picture originally proposed by Bragg'1*”) and supported 
by Burgers“), according to which the subgrains are essentially formed 
during the deformation process itself, rather than during subsequent 
annealing.* In view of microscopic observations of subgrain behaviour'!45» 
and of the occurrence of recrystallization in practically strain-free material, 
it appears that Cahn’s theory, together with the mechanism proposed by 
Cottrell, represents the best basis available at present to describe the 
formation of recrystallization nuclei in deformed single crystals and in 
severely deformed polycrystalline material. 

Microscopic observations by Sperry and Beck(4® showed that in 
polycrystalline high purity aluminium, after deformation of less than 
approximately 40°, by rolling, the formation of new nuclei on annealing 
is relatively rare and that volume elements relatively free of substructure 
form instead by the ‘strain-induced migration’ of already existing 
boundaries. Figure 57 (Plate 33) shows such an area, that became 
largely free of substructure when the boundary migrated from its original 
position 1 to a new position 2. The microscopic technique used (oxide 
film and sensitive tint illumination) indicated that the area between 1 
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* On the other hand, the theory of nucleation proposed in 1949 by Beck(18® 
was based on polygonization, and it should be, therefore, abandoned. 
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and 2 has the same orientation as that of the original grain A, from which 
the boundary moves away. This conclusion was confirmed by means of 
back reflection Laue patterns. Measurements of the rate of this kind 
of boundary migration were carried out“ by an annealing treatment 
interrupted at suitable intervals and by developing surface markings 
showing the successive boundary positions (fig. 58, Plate 34)- Within 
the accuracy of the measurements, no incubation period could be detected ; 
growth started immediately upon the immersion of the specimens in the- 
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Isothermal strain-induced migration at 350°c of several boundaries in the pure 
Al specimen shown in fig. 58. [Ref. (11).] 


molten salt bath at the annealing temperature, and in most cases the 
initial rate of boundary migration was the highest (fig. 59). The volumes 
essentially free of substructure so formed eventually attain dimensions 
commensurate and even larger than the matrix grain size, and on further 
annealing, they begin to absorb other surrounding grains. ie 

It seems reasonable to assume that under conditions similar to those- 
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used in Anderson and Mehl’s2® experiments, that is when both the 
initial grain size and the strain is small, this process represents the pre- 
dominant ‘nucleation’ mechanism. This conclusion is confirmed by 
the findings of Tiedema‘!4*), who has shown by means of x-ray diffraction 
techniques that the central area (‘ nuclear spot ’) of large grains produced 
in pure Al strips by the strain-anneal method often consists of many 
subgrains of an orientation similar to the main portion of the grain. 
This is just what may be expected if the large grain is nucleated by the 
mechanism described, since the recrystallized grain can not absorb the 
deformed grain (A in fig. 57) of nearly the same orientation from which 
it originates“%5), Tiedema also observed that the main portion of the 
large recrystallized grains often consists of two or more slightly disoriented 
domains. As noted by Crussard”), this, too, is in agreement with the 
proposed mechanism, since it has been observed microscopically''*® that 
the relatively substructure-free area swept by the migrating boundary 
(fig. 57) in effect usually does contain a few low angle boundaries. As 
seen in fig. 58, the origin of these boundaries may be traced back to 
independent nucleation along the original grain boundary of several 
‘buds’ of nearly the same orientation, which later meet and become 
joined by low angle boundaries. Apparently, each ‘ bud ° originates ina 
single subgrain, S, of deformed grain A (fig. 60 (a)), directly adjacent to 
the boundary of deformed grains A and B. If this boundary is mobile 
by virtue of the large disorientation between A and B and, if opposite to 
.subgrain 8 of grain A the subgrain size in grain B is much smaller, so that 
there is an unbalanced ‘ pull’ in the direction of B, the conditions are 
favourable for the process to start immediately on annealing. In the 
initial stage of the boundary migration, shown in fig. 60 (6), the driving 
energy is the residual portion of the stored energy of cold work (which 
may be merely the energy of the substructure) in grain B less the energy 
required to extend the migrating high angle grain boundary. After the 
linear dimensions of the substructure-free volume become approximately. 
twice as large as the average matrix grain size, the driving energy is 
increased by the ‘ pull’ of the high angle grain boundaries between 
adjacent matrix grains, which in itself now compensates for the energy 
required to expand the migrating boundary. Consequently, one may 
assume a period of relatively slow growth during the initial stage, when the 
net driving energy is low, and subsequent fast growth after the ‘ re- 
crystallized ’ grain becomes several times larger than the matrix grains. 
On this picture, the ‘incubation period’ would correspond to the 
initial stage of slow growth. The peculiar variation of the rate of 
nucleation with annealing time, as reported by Anderson and Mehl‘29, 
fig. 44, could then be interpreted as the frequency distribution for the 
periods of slow growth for the various ‘ nuclei ’. 

Both the mechanism of ‘ nucleation’ just described, and the theory 
-of Cahn discussed previously, have in common that they do not presuppose 
any lattice stress other than that implied by the presence of subbound- 
-aries. In both cases, the site of ‘ nucleation ’ is neither a high stress point 
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nor a low stress point, but rather it is characterized by relative freedom 
from substructure. According to the Cahn—Cottrell picture, the sub- 
Structure-free volume of sufficient size develops on annealing by subgrain 
growth, and the high energy boundary with the matrix develops gradually 
as the disorientation increases during the slow initial growth. In the 
other mechanism, both the subgrain of sufficiently large size with respect 
to its neighbours and the high angle boundary are present already in the 
initial state, and the migration of the high angle boundary may start 
immediately. The apparent incubation period observed here macro- 
scopically, for instance by Anderson and Mehl2, has a rather different 
physical significance from the incubation period in accordance with 


Fig. 60 


‘Schematic drawing of grains A and B with their substructure (a) before, and 
(6) after the beginning of strain-induced boundary migration. This 
boundary migration is ‘nucleated’ by subgrain S of grain A, which is 
particularly large and is directly adjacent to the grain boundary. 


the Cahn—Cottrell theory. The ‘ strain-induced boundary migration ° 
mechanism is apparently prevalent at low to medium deformations, and 
it certainly requires a polycrystalline matrix. Presumably, its prevalence 
may increase with decreasing matrix grain size, since with a larger 
number of grain boundaries favourable combinations of subgrain size 
contrast and disorientation should become more frequent. On the other ~ 
hand, the Cahn—Cottrell mechanism may predominate in the case of 
single crystals, or of polycrystalline metal severely deformed, where most 
of the recrystallized grains originate within the deformed grains, rather 
than directly at the boundaries. The frequency of coincidence of the 
required local gradients of subgrain size and of orientation may be expected 
to increase with increasing deformation. 

Microscopic observations*) with the oxide film-polarized light 
technique revealed the fact that, at least in severely rolled high purity 
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Al of large grain size, subgrain growth may lead directly to structures: 
that under the microscope appear very much like recrystallized material. 
Under such conditions, the subgrains were observed to grow to sizes of 
the order of 0-02 to 0-03 mm; they can not then be distinguished from 
ordinary grains by their size. Since the disorientation among neighbours 
tends to increase with progressive subgrain growth, this criterion also- 
becomes tncertain. It is quite possible that in some cases subgrain 
growth may in this manner change over imperceptibly into what might 
be designated as gradual grain growth, leading to even larger grain sizes, 
without ever passing through a discernible stage of nucleation and growth, 
characteristic of primary recrystallization. 


(6) Coarsening (Secondary Recrystallization) 

Since, as it is now quite generally recognized, coarsening does not depend. 
on the presence of lattice imperfections in the matrix”, the Becker 
condition must be interpreted entirely in terms of interface energies. 
Accordingly, a stable, growing nucleus must be large enough so that the 
surface energy gained by shortening the interfaces between the matrix 
grains adjoining the nucleus must be at least equal to the surface energy 
lost by the corresponding extension of the nucleus-matrix interface. 

In the case of coarsening in rolled and recrystallized silver, Rosi,. 
Alexander and Dube) found that the nucleus must be at least twice as 
large as the average matrix grain size. When the matrix has a single 
orientation texture, so that most of the matrix grain boundaries have low 
energies, while the interface between the coarse grain and the matrix 
consists of high energy boundaries, as in the case of coarsening in cube- 
texture copper, the nucleus must be presumably much larger. Such 
considerations practically exclude all possibilities of nucleation by thermal 
activation, even if it is assumed that nuclei form at the intersections of 
twin-boundaries, grain boundaries and specimen surface, as proposed by 
Kronberg and Wilson‘!4®).* The only reasonable picture then is that the 
“nucleus ’ for coarsening is a grain already present in the matrix 
(‘preformed nucleus’), that is capable of growth by virtue of its size 
and its orientation®*®.+ The grain structure and orientation, as set up 
in primary recrystallization and in subsequent grain growth, must be 
responsible for the coarsening behaviour on further annealing. 

The interpretation of the incubation period in coarsening depends on 
the type of coarsening considered. In orientation-dependent coarsening 
(e.g. in ‘ cube texture ’ copper), where general grain growth is prevented 
by the low mobility of most of the matrix grain boundaries, the incubation 
period must be due to the initially low rate of growth of the ‘ nuclei ” 
when they are still quite small. However, as shown by Rosi, Alexander, 
and Dube“), an increase of the rate of growth with the size of the 
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# Following a calculation by E. Orowan (private communication), the 
necessary activation energy would be several orders of magnitude larger- 
than that available at ordinary annealing temperatures. 
} This orientation effect will be discussed in the following sections. 
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growing grain is encountered also in cases where general grain growth is 
prevented by second phase inhibition (‘ mhibition-dependent coarsening ’). 

In addition to the initially low rate of growth when the ‘ nuclei’ are 
still small, in the case of inhibition-dependent coarsening the incubation 
period may be also affected by the rate of re-solution and coalescence of 
the dispersed inhibiting particles. When the inhibition is initially very 
Strong, as in the experiments with Al-—Mn alloys !!8), the decrease of 
inhibition during annealing is undoubtedly the predominant factor in 
determining the incubation period. The local conditions of decrease of 
inhibition. play here an important role in the selection of sites for ‘nuclea- 
tion ° by permitting certain large matrix grains to grow in preference to 
others that may be equally large. The picture of ‘ nuclei’ for coarsening 
given above requires, therefore, some correction. In general, the 
“nucleus ’ for coarsening is a grain already present in the matrix, selected 
for growth by virtue of its large size, its orientation, and the relatively 
weak inhibition at its boundaries. Of these three selective factors the 
large relative size is always required, the importance of orientation 
increases with increasing sharpness of the matrix texture and that of 
local conditions of inhibition with the general level of inhibition in the 
material. 


4.5. Problems of Growth Rate 

(a) Theory of Boundary Migration 

A rather detailed discussion of the problems of boundary migration 
was given in the recent review by Burke and Turnbull?). Consequently, 
the following comments are limited to a very brief summary. Anderson 
and Mehl’s quantitative data‘? for the rate of growth of recrystallized 
grains in a fine grained, slightly deformed matrix of high purity Al, as 
a function of annealing temperature, lend themselves very well for 
theoretical analysis. It was pointed out by Mott‘! that the measured 
rates of growth are much too high for a theory in which the thermally 
activated elementary process is the transfer of a single atom from a 
erystal to the boundary layer. He found, however, that the difficulty 
can be eliminated if one assumes that the thermally activated elementary 
process is the transfer of a group of about twenty atoms into the boundary 
zone (‘melting’). This theory accounts also for the surprisingly high 
activation energy for growth observed by Anderson and Mehl. 

Turnbull@34) and Burke and Turnbull treated the problem from the 
point of view of the absolute rate theory and they proposed an alternative 
interpretation, where the observed large variation of the rate of growth 
with temperature is supposed to be the result of progressive re-solution 
with increasing annealing temperature of growth-inhibiting dispersed 
particles. As discussed in § 4.3 of this chapter, the results of A, Suzuki) 
suggest that the metal used in Anderson and Mehl’s investigations may 
in fact have contained some dispersion of a second phase. The apparent 
dependence of the measured activation energy for growth on the impurity 
content was noted by Suzuki. No doubt this correlation could be further 
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clarified if the impurities forming dispersed particles would be taken into 
account separately from those in solid solution. That the presence of 
second phase forming impurities does in fact increase the measured 
activation energy for the rate of boundary migration is seen when compar- 
ing the activation energy for grain growth in high purity brass (approxi- 
mately 40 kcal/mole 8: 112)) with that in brass of commercial purity 
(approximately 60 kcal/mole! 15»), ‘These observations tend to support 
Burke and Turnbull’s point of view. On the other hand, since the effect 
of inhibition on the rate of boundary migration in recrystallization may 
be expected to be much larger if the preceding deformation is small 
than if it is large, the ‘inclusion hypothesis’ of Burke and Turnbull 
would require a great decrease in the measured apparent activation 
energy with increasing deformation. However, as noted by Burke and 
Turnbull themselves, the decrease of activation energy for the rate of 
growth actually observed by Anderson and Mehl is relatively small 
(fig. 47), in disagreement with the theory. Whether or not the unex- 
pectedly high entropies of activation for growth calculated from Anderson. 
and Mehl’s data‘134: 10 are in agreement with the theory can not be 
determined on the basis of available data. It now appears that a fully 
satisfactory solution of these problems will not be possible without. 
additional data. What is needed is accurate growth rate measurements. 
in a matrix, which is itself stable against structural changes during the 
growth anneal, and also free from dispersed phases. 


(b) Orientation Dependence of Boundary Mobility 

Dunn) measured the rate of growth of silicon—iron crystals as a fune- 
tion of their orientation relative to a fine grained matrix with essentially 
single orientation texture. He found a broad maximum of the rate of 
growth for a relative orientation corresponding to a [100] rotation of 45°. 
Independent observations by Lacombe and Berghezan'®®) with high 
purity aluminium and by Tiedema, May and Burgers'!®® with aluminium 
of commercial purity indicate that in recrystallization after small strains 
and in coarsening the large new grains do not absorb all matrix grains. 
at the same rate, but that matrix grains with orientations near that of the 
large growing grain, or nearly twin-related to it, may persist over long 
periods of annealing, even when completely enveloped by the new grain. 
Since the grain boundaries corresponding to these two relative orientations. 
have low interface energies, in the case of coarsening in commercial 
aluminium the above observation might be explained by assuming that 
the driving energy for the movement of these low energy boundaries is not: 
sufficient to overcome the inhibition. However, Lacombe’s results with 
recrystallization in high purity aluminium, where the inhibition, if any, is 
very much smaller, and at least a considerable fraction of the driving 
energy is due to the plastic strain, may be explained only on the assump- 
tion that the mobility of grain boundaries depends on the disorientation — 
across them. Other observations* also support the view that, in 
general, low energy boundaries have low mobility. However, on the 
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basis of Guinier’s results'4), it has been suggested) that low energy 
boundaries of a special type (essentially dislocation wall type subbound- 
aries) may have high mobility in spite of their low energy. That large 
differences in mobility may exist between low energy boundaries of this 
special type and those of a more general kind is quite plausible, since the 
atomic mechanism of interface migration is different in the two cases. 

It has been shown by Rathenau and Custers'5®, Kronberg and 
Wilson"? and Beck and Hu‘!>”) that in face centred cubic metals the 
fastest growing grains on coarsening after random artificial nucleation 
in a polycrystalline matrix with single orientation texture are those with 
an orientation relative to the matrix corresponding to a 30° rotation 
around a [111] axis for Ni—Fe alloy and Cu, and to a 40° rotation around 
the same axis for Al (fig. 61, Plate 35). Beck, Sperry and Hu®) later: 
found the same orientation relationship for maximum boundary mobility 
in the recrystallization of slightly deformed high purity aluminium crystals 
artificially nucleated by scratching (fig. 62, Plate 35). Although additional 
work would be certainly very desirable to determine more quantitatively 
the orientation dependence of boundary mobility, the results of the above 
investigations are rather definite in establishing the orientation relation- 
ship for maximum boundary mobility in various f.c.c. metals. 

This orientation (corresponding to 30° or 40° rotation around [111], 
which is approximately half way between the two relative orientations. 
with very low interface energy mentioned above (corresponding to 0° and 
60° rotations around [111], respectively), may be expected to correspond 
to a maximum of the interface energy. That there may be, at least 
qualitatively, a correlation between boundary mobility and interface 
energy, has been suggested'1#) on the assumption that the rate determin- 
ing process in boundary migration is akin to self diffusion along the 
boundary, and that the rate of boundary self diffusion increases with the 
boundary energy. Experimental confirmation has since been obtained‘1> 
of the proposition that the rate of self diffusion is higher in high energy 
boundaries than in those of low energy. Turnbull’s calculation”, 
showing that in silver the free energy of activation for boundary migration 
is in fair agreement with that for grain boundary self diffusion, suggests. 
that the elementary process in boundary migration may indeed be the 
same as in grain boundary self diffusion. More quantitative data would 
be desirable in order to test the above correlation. 

Another attempt to rationalize the orientation dependence of boundary 
mobility might be based on the findings of Kronberg and Wilson‘*®), who. 
observed that for a series of discrete rotation angles around a [111] axis. 
and also around a [100] axis, in f.c.c. structures the atomic movements 
necessary for transition between the two disoriented lattices are particu- 
larly short and, therefore, presumably easy to accomplish, and 
approximately 14% of the atomic sites actually coincide, requiring no 
movement at all. Consequently, these discrete orientations might be 
assumed to correspond to maximum boundary mobility. Although 
two of the predicted [111] rotations are fairly close to the observed 
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orientations of greatest boundary mobility for Cu, respectively Al, 
a serious difficulty with this theory lies in the fact that most of the 
predicted orientations of maximum boundary mobility have not’ been 
confirmed experimentally. Furthermore, there is apparently no reason 
to assume that the relative orientations given by Kronberg and Wilson 
are the only ones possessing the described geometric features. As a 
matter of fact, in the case of twin orientation as many as one third of 
all atoms in the two f.c.c. lattices fully coincide and all the other 
corresponding sites are also less than the distance of closest packing 
apart when projected on the (111) plane. In spite of this, as mentioned 
above, boundaries between twin-related lattices are notoriously immobile. 
It may be, therefore, concluded that the suggestion just discussed is 
not likely to prove satisfactory. 


4.6. Annealing Textures 


It was noted above that the migration of high angle grain boundaries 
(recrystallization) involves local reorientation, as a grain of a certain 
orientation is replaced by a new grain of a different orientation. The 
observation that metals annealed after plastic deformation generally 
show a pronounced preferred orientation and that this is usually different 
from the deformation texture, indicates both that * nuclei’ for the new 
grains are available in orientations corresponding to the annealing texture, 
and that these ‘ nuclei’ have the capacity to grow at the expense of the 
matrix during annealing. A theory of annealing textures must account 
for the scarcity of new grains in other than the preferred orientations. 
This scarcity may be attributed either to a shortage of nuclei in the 
corresponding orientations (oriented nucleation theory), or to their low 
rate of growth%* 155, 157) (oriented growth theory). 

The assumption that sharp annealing textures result from the availabi- 
lity of nuclei only in specific orientations is, in general, extremely difficult 
to verify directly, since nuclei are by nature small, and can not be 
identified as such until after they have grown and were, therefore, subject 
to growth selectivity. An oriented nucleation theory of annealing 
textures might, however, still have merit, if it would allow the prediction 
of annealing textures on the basis of some general principle regulating 
nucleus orientations, and applicable to a large number of individual 
cases. So far, no such general principle has emerged. Burgers and 
Louwerse™%) attempted to formulate one by assuming that nuclei in 
plastically strained metals are formed in regions of ‘ local curvatures ’, 
which are disoriented with regard to the main part of the lattice by 
rotation around a [112] axis parallel to the active slip plane and per- 
pendicular to the active slip direction. They were able to analyse 
recrystallization textures in compressed and rolled aluminium single 
crystals in terms of such rotations, at least approximately. However, 
this particular oriented nucleation principle certainly did not prove 
generally applicable. Burgers himself, in a later publication™®») found 
it necessary to introduce a new oriented nucleation principle in order 
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to explain the formation of the ‘ cube texture ’ (orientation : (100) [001}) 
in f.c.c. metals. According to this new principle, those elements of the 
rolled structure which have the cube orientation are the first to 
‘ polygonize * (i.e. they sustain subgrain growth at the highest rate) 
and therefore they are the first to become capable of acting as nuclei. 
However, this principle is not generally applicable either, since it does 
not account for annealing textures other than the cube texture. Actually, 
rolled 70-30 brass‘'®), cross-rolled copper''®), rolled single crystals of 
aluminium'?%, 157 and of copper‘t®) etc. are known to recrystallize in 
other orientations. Even in straight-rolled copper, the amount of cube 
texture formed on annealing becomes appreciable only after very high 
deformation’®», A third oriented nucleation principle was introduced 
by Kronberg and Wilson? to account for the texture obtained by 
coarsening in cube texture copper. It is apparent that, in the absence 
of a generally applicable oriented nucleation principle, in various specific 
instances new ad hoc assumptions are required regarding the orientation 
of the nuclei in order to explain the experimentally found annealing 
textures. So far, the oriented nucleation theory apparently did not 
allow the prediction of an unknown annealing texture. 

The three specific oriented nucleation principles mentioned are subject 
to grave objections. The [112] rotations predicted by Burgers and 
Louwerse’s principle were not confirmed by other investigators either 
for compressed%°) or rolled aluminium single crystals“®” or for rolled 
copper single crystals‘). Furthermore, it was shown” that the 
results of Burgers and Louwerse themselves may be described more 
accurately by [111] rotations. The principle of Burgers, Liu and Tiedema 
predicts the formation of the cube texture on annealing, of deformed f.c.c. 
metals, certainly at least when a measurable fraction of the deformed 
material is in the cube orientation*. The investigation of Merlini and 
Beck‘1®) showed that in cross-rolled copper there is present a measurable 
fraction of material in the cube orientation, but the cube texture 
nevertheless does not form on annealing. This condition obtains in 
spite of the fact that the small fraction of cube oriented material in 
straight-rolled copper, where the cube texture does form on annealing, 
is only approximately half of that in cross-rolled copper. The principle 
of oriented nucleation advanced by Kronberg and Wilson*®) to explain 
the orientation of grains resulting from coarsening in cube texture copper 
predicts a large number of orientation relationships, of which actually 
only three are néar orientations experimentally found. No suggestion 
was so far advanced to account for the absence of coarse grains near the 
other predicted orientations. 

On the other hand, the oriented growth ony is based on the ex- 
perimentally observable orientation relationship that corresponds to 
maximum boundary mobility. As stated in the previous section, for 
f.c.c. metals this orientation relationship may be described as a rotation 
of 30° or 40° around a [111]axis. By analogy, for h.c.p. metals maximum 
boundary mobility may be expected for a disorientation of 30° around 
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the [0001] axis1¢®), When the deformed metal has a fairly sharp single 
orientation texture, the oriented growth theory predicts the formation 
of an annealing texture consisting of all the crystallographically equivalent 
orientation corresponding to maximum boundary mobility. In the 
case of copper, there are 8 equivalent orientations related to the ideal 
orientation of the deformation texture according to 30° rotations around 
the four [111] axes. Liu and Hibbard®) found that, when a copper 
single crystal is heavily rolled on a (110) plane, in a [112] direction, 
‘ and a sharp single orientation deformation texture is obtained, the 
resulting annealing texture in fact consists of the eight orientations 
mentioned above. The same eight relative orientations predominate on 
coarsening in f.c.c. Ni-Fe@5® and in Cu®®) with a single orientation cube 
texture matrix. In the case of hexagonal close packed metals, the expected 
orientation relationship corresponding to maximum boundary mobility 
is satisfied by one orientation only. Recent work of McGeary and 
Lustman'8®) showed that a single orientation deformation texture may 
be obtained with polycrystalline zirconium by 60°-cross rolling. They 
also found that the annealing texure developing from this deformation 
texture corresponds to a rotation of 30° around the [0001] axis, as 
predicted on the basis of the expected maximum boundary mobility®®. 
McHargue et al.“ recently obtained an essentially single orientation 
deformation texture by straight-rolling a hexagonal close packed alloy 
of Ti+3-8% Al, and they found that the annealing texture, here too, 
was the one just described(17)). 

When the deformation texture is more complex, consisting of more 
than one orientation and of significant orientation spreads, the annealing 
texture corresponding to maximum growth rate in such a complex 
matrix is, in general, difficult to predict uniquely, without a detailed 
knowledge of the orientation dependence of boundary mobility. On 
the basis of the oriented growth theory, one may expect that the crystal 
orientations most favourable for growth in a complex matrix, and there- 
fore actually developmg in the corresponding annealing texture, will 
not, in general, necessarily have exactly the ideal orientation relationship 
with any one of the matrix texture components, but rather they may 
deviate to some extent because of the presence of other matrix texture 
components. Furthermore, several of the orientations ideally related 
to certain individual matrix texture components may be so unfavourably 
oriented for growth at the expense of other matrix components or of 
important matrix orientation spreads that they do not grow to a signi- 
ficant component in the annealing texture(!6®), 

Despite the conditions described, in most of the known annealing 
textures of f.c.c. metals the various components of the annealing texture 
come surprisingly near to having the ideal orientation relationship with 
respect to at least one of the matrix texture components. In cross-rolled 
Cu) each of the eight main components of the annealing texture 
(type: (519)[4, 11, 1]) is exactly so related to one of the two main 
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deformation texture components (type: (110) [223]). The other eight 
potential annealing texture components, which actually do not form, 
are distinctly less favourably oriented* in relation to the second deformation 
texture component than the eight components that do form (table 1). 
In straight rolled 70-30 brass, the formation of the (112) [110] annealing 
texture, which would be ideally oriented with respect to both main 
deformation texture components (type: (110) [112]), is apparently 
prevented because it is very unfavourably oriented with respect to the 
principal orientation spread of the deformation texture (table 2). 
Instead, four main components develop on annealing, each one of which 


Table 1. Cross-rolled Copper 


Orientation relation- | Angular distance Angular distance 
ship with second of rot. axis from of rot. axis from Angle of 
main def. text. nearest [111] second nearest [111] | rotation 


comp’t. (110) [557] 


For actual recryst. 

text. components, re- 

lated to first main 

def. texture comp. by 16° 54:5° 49° 
30° rotation around 
[111] near _ centre, 
(519) (4. 11, 1] 


For hypothetical re- 

eryst. text. comp.’s, 

related to first main 

def. texture comp’t. 33°5° 39° 35° 
by 30° rotation around 
[111] at periphery 


is fairly close to having the ideal orientation relationship with one of 
the matrix components) and it is not very unfavourably oriented in 
relation to the other one and to the principal orientation spread (table 2). 
Also, in straight-rolled h.c.p. zirconium, which has a two-component 
deformation texture, the recrystallization texture‘®®, 167) corresponds. 
quite well to a rotation around the [0001] axis of each of the deformation 
texture components and a deviation from ideal reorientation appears 
only as a slight difference in the angle of rotation (20° instead of 30°). 
It may be concluded, therefore, that the principal effect, if any, of the: 


* The angular distance of the axis of the rotation, which relates the two 
orientations, from the nearest [111] axis may be taken in a first, rough approxi- 
mation as an index of the extent to which a certain orientation approximates 
that most favourable for growth (except near the identical and the twin 
orientations). 
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presence of additional texture components or major orientation spreads 
in the matrix is to prevent the formation of certain annealing texture 
components, which are unfavourably oriented for growth at their expense. 
The deviations from the ideal orientation relationship with one of the 
matrix texture components (as a result of the presence of other matrix 


Table 2. Rolled 70-30 Brass 


Orientation relationship with 


first main def. comp’t. 
(111) [112] 


Angular distance of 
rotation axis from 


Angle of 
anes seis: Rotation 
{111] fq] 
| For actual recryst. texture Be ee Qg0 
(225) [734] 11-5 63-5 32° 
| For hypothetical recryst. 0° 705° 30° 


| texture (112) [110] 


second main def. comp’t. 
(110) [112] 


| For actual recryst. texture é nese 
| (295) (734) 155 56-5 46-7° 


| For hypothetical recryst. % fee 
texture (112) [110] 0 70:5 30° 


centre of principal def. text. spread 
(110) [001] 


For actual recryst. texture 3 
(225) [734] a = 


For hypothetical recryst. é 
texture (112) [110] 43-5 46° 
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orientations) are usually fairly small. This being so, at least approximate 
predictions usually can be made. The prediction made in the case of 
cross rolled copper prior to the quantitative determination of the annealin 
texture proved to be precise. : 
The same conclusions apply also for the formation of the cube texture 
((100), [001]) in rolled f.c.c. metals. This annealing texture deserves 
particular attention, because of its unique sharpness, and because of the 
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rather detailed knowledge now available concerning its formation. It is 
known”) that in copper the percentage cube texture formed increases 
with decreasing penultimate grain size and with increasing rolling reduction, 
in such a way that various combinations of these two variables leading 
to the same rolled grain thickness give rise to approximately the same 
fraction of cube oriented material, if the annealing treatment remains 
the same. It is also known that the fraction of the cube texture com- 
ponent increases with increasing annealing temperature!73) or time(!7®, 
Finally, it has been observed(7® that the formation of the cube texture 
can be suppressed to a great extent if the rolled Ni-Fe alloy strip is etched 
to a thickness of approximately 1 micron before annealing. It was 
pointed outt*® 16) that not only the orientation of the cube texture 
and its sharpness but also the effect of all the above mentioned factors 


Table 3. Cube Texture 


Orientation relationship with 


Cu, main def. Cu, centre | brass, main | brass centre 
text. _ of spread | def. text. | of spread | 
(5, 8, 14) [11, 5, 7]* | (112) [111] | (110) [112] | (110) [001] } 


Angular distance 


of rot. axis from BiG 23° ie 55° 

nearest [111] 

Angular distance 

of rot. axis from g9° 50° 50° 55° 

second nearest 

[111] 

Angle of Ea ° RRO 45° 
Rotation ab “ ph 


* This ‘ideal orientation’ is more accurate than the (3, 5, 8) [3, 5, 2] 
orientation originally given by Hu, Sperry and Beck (see also ref. 186). 


on its formation may be readily explained on the basis of the oriented. 
growth theory. The cube orientation is not far removed from the ideal 
orientation relationship for maximum boundary mobility with the main 
deformation texture components in the rolled f.c.c. metals in which it does 
form on annealing (table 3), and it is in fact unique in having such an 
orientation relationship with all four deformation texture components'*#°). 
For this reason, the growth of recrystallized grains in the cube orientation. 
is favoured over others after the new grains have grown large enough 
to be in effective contact with matrix elements (e.g. grains, kinkbands) 
in more than one and, preferably, in all four matrix orientations. This. 
critical size is reached relatively earlier and. therefore, selectively favoured. 
growth of the cube grains may take place over a larger fraction of the 
total annealing period, if the elements of the rolled matrix are smaller 
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(hence the effect of penultimate grain size and of the rolling reduction) 
and if the annealing temperature is higher, or the annealing period longer. 
On the other hand, factors inhibiting the growth of the new grains (certain 
impurities) or preventing them from coming into effective contact with 
more matrix texture components (etching to small thickness) will tend 
to suppress the formation of the cube texture. 

An objection raised(!7”) against the oriented growth theory was based 
on the observation that the maximum of the grain boundary energy, 
and presumably also of boundary mobility, as a function of disorientation 
is usually quite broad. Grains deviating from the orientation for 
maximum mobility by 15°, or so, still grow at relatively high rates. 
In view of such observations, it has been doubted that the orientation 
dependence of boundary mobility is sufficient to explain the sharp 
annealing textures known to occur. It has been noted’), however, 
that the maximum becomes sharper if the increase of volume of recrys- 
tallized material rather than the linear rate of boundary migration is 
considered as a function of orientation; on the basis of volume, the 
pertinent variable is the third power of the mobility. If it is assumed 
that the ‘ incubation period ’ for nucleation is the period of slow * invisible 
growth’ and that its rate depends on the orientation relative to the 
matrix in a manner similar to that of visible growth, then the apparent 
rate of nucleation will be also similarly orientation dependent. Assuming 
that both G and WN are proportional to the boundary mobility, the latter 
enters approximately at the fourth power into the exponential term of 
the formula giving the recrystallized volume. This is seen particularly 
well in the simplified form of the Johnson—Mehl equation (9), 


F()=1—exp(— FN.) ME, oS 


where G and WN are assumed to be independent of time. It is clear that 
the orientation dependence of the recrystallized volume may be relatively 
sharp, even if that of the mobility is not. An other effect, which further 
exaggerates the consequences of small differences in mobility, is due to 
the impingement on each other of the recrystallized grains78), As seen 
in fig. 62, fast growing grains ‘ cut in front’ of grains with boundaries 
of lower mobility, so that the latter are ‘ pinched off’. As a result of 
this “pinching off’ due to impingement, grains with boundaries of 
maximum mobility can grow to disproportionately larger sizes than other 
grains, even if the differences in mobility are slight, provided that the 
number of competing new grains per unit volume is large enough to make 
the competition keen.* An additional reason, on the basis of oriented 
—_—_—ee aa ee eee 
Be Burgers and Tiedema‘!?®) proposed an explanation of the varying grain 
sizes obtained in this experiment on the basis of an assumed variation in the 
distance of the ‘ nuclei’ from the edge of the specimen. It is clear, however, 


that this explanation does not account for the most significant effect, namely 
the observed correlation between size and orientation of the new grains. 
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growth, for the formation of sharp annealing textures comes into play 
if the matrix texture has more than one component. As already 
mentioned, the favoured orientations in this case are those for which the 
effective overall mobility in the complex matrix as a whole is largest. 
Such a “compromise texture ’ (e.g. the cube texture) may correspond to 
high mobility effective only in a very narrow range of orientations, 
because of the partially overlapping restrictions on mobility imposed by 
the various individual matrix texture components(!6®. In view of these 
considerations, the objection raised appears to be ineffective. 

An other objection raised against the oriented growth theory(!77: 129) 
was that it is incompatible with the observed effect of annealing tempera- 
ture on the resulting texture. Actually, the information now available 
on this effect is rather limited. It is known that in rolled Cu the fraction 
of cube texture increases with the annealing temperature73) and that 
in cold rolled 70-30 brass annealing at 650°C results in a (110) [112] 
type texture, in contrast to the (225) [734] texture obtained at around 
340°c. In both of these cases the textural changes may be interpreted as 
the results of grain growth after recrystallization, the reorientation being 
in good agreement with the principles of oriented growth'!4>: 157), A tex- 
tural change; similar to that in brass, occurs also in silver, where it is 
associated with discontinuous grain growth‘), Jn such cases where the 
annealing texture changes with the annealing temperature, it is possible 
to test whether or not the texture produced at the lower annealing 
temperature in fact changes on further annealing at the higher temperature 
into the same texture that is characteristic of the high annealing tempera- 
ture. Ifsuch experiments would give positive results, they would confirm 
that the texture change is a result of additional growth. If the results 
turned out to be negative, explanation might be attempted for instance 
along one of the following lines: (1) Precipitation of dispersed 
particles of a second phase takes place during the low temperature 
anneal, inhibiting further growth, at least initially, on subsequent anneal- 
ing at higher temperature, (2) the orientation dependence of boundary 
mobility changes with tenperature, or (3) the orientation dependence of 
nucleation rates changes with annealing temperature. If the incubation 
period is considered as the period of slow ‘invisible growth’, then the 
effect on the rates of nucleation also may be a consequence of the change 
in the orientation dependence of boundary mobility. New data will be 
certainly necessary in order to allow a choice between the various 
possibilities. 

A specific limitation of the oriented growth theory of annealing textures 
should be clearly noted. As seen above, well defined annealing textures 
are predicted only under circumstances where a reasonably large number of 
well disoriented nuclei are available. It was shown by various investi- 
gators‘179 180, 181) that, in the case of slightly deformed single crystals, 
where nuclei, and particularly well disoriented ones, are rather scarce 
the orientation of the recrystallized grains bears no simple relationship 
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with the matrix orientation. Apparently, the various nuclei are so far 
apart that competition between them is ineffective ; almost any nucleus 
can grow, with the exception of those in nearly identical or nearly twin 
orientations. This condition may be so exaggerated that active nuclei 
in favourable orientations are practically absent. For instance, carefully 
prepared specimens of cube texture copper may be annealed for long 
periods at 1000°c without coarsening%®, since the texture is so perfect 
that grains disoriented more than a few degrees are practically absent. 
If the specimen is annealed long enough, coarsening eventually does 
occur and the large grains have orientations in the order of 10° removed 
from that of the matrix. Kronberg and Wilson4® and Dunn and 
Sharp'82) found some coarse grains rotated about 20° approximately 
around (100) axes with respect to the cube texture matrix in Cu, even 
when <111)-rotated coarse grains also formed. It may be assumed that 
in borderline cases, where nuclei in orientations most favourable for growth 
are still fairly scarce, coarse grains also are able to develop in somewhat 
less favourable orientations. From the observed orientations one may con- 
clude that, among the orientations differing from the matrix by up to 
about 20°, those corresponding approximately to <100)-rotations are 
relatively most favourable for-growth. Various investigators'1®’: 149. 189) 
noted that in cases of scarce nucleation the artificial introduction of 
sufficiently random nuclei gives rise to the growth of grains in the usual 
orientations corresponding to maximum boundary mobility. This 
furnishes additional proof that the special conditions discussed above 
are indeed characterized by a scarcity of nuclei. <A different situation 
where nuclei are also absent, even though the matrix is capable of support- 
ing the growth of ‘ nuclei’ artificially introduced, may arise in cases of 
inhibition by a dispersion. As described in § 4.3 (c), a larger grain may 
be able to grow in such a matrix at levels of inhibition and at annealing 
temperatures where the boundaries between the small matrix grains are 
effectively immobilized. Ifa single large grain is produced by local over- 
heating, under properly chosen conditions this may be made to grow 
and to absorb the whole specimen into a single crystal. The method was 
used by Tiedema‘*®) and by Dunn and Lionetti‘’’ for producing 
crystals of predetermined orientations. Clearly, then, under carefully 
controlled conditions it is possible to produce in various ways orientations 
that are not determined by selective growth. 

In a study of coarsening in rolled silicon—iron single crystals, Dunn(18®) 
traced the orientation of the coarse grains to four larger than average 
size grains present in the structure at a stage of annealing before visible 
coarsening. He stated, although published no supporting evidence, that 
these particular grains grew to their large size at the expense of the 
deformed single crystal (primary recrystallization) rather than of the 
surrounding recrystallized grains. Dunn also stated that the orientation 
of these four grains was representative of that of all large grains in the 
structure deviating from the main orientation of the matrix. Thus, in 
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his view, the grains in question served as oriented nuclei for coarsening, 
without regard to the dependence of boundary mobility on the orientation 
of the growing grains relative to the recrystallized matrix. If his view 
is accepted, this case may be considered as somewhat similar to the special 
situations mentioned in the previous paragraph. The cause of the large 
size of recrystallized grains of the particular orientation in question was 
not explained. More recently, Dunn''8® attempted to generalize this 
mechanism to other cases of coarsening. However, there appears to be 
no experimental evidence, nor any evident reason for the assumption 
that in all these cases the structure obtained by primary recrystallization 
contains sufficiently large grains only in the orientations predominating 
in coarsening. For example, unless an additional selective principle is 
introduced in the case of cube texture Cu, one might expect that the 
minor texture components assumed by Dunn to furnish nuclei for 
coarsening would occur with equal frequency an all the 28 orientations 
which, together with the cube texture itself, may be derived from the 
four rolling texture components. (If one prefers to use the 22° and 38° 
[111]-rotations proposed by Kronberg and Wilson, as does Dunn, the 
number of differently oriented minor textures expected is much larger 
yet.) Since the coarsening texture actually consists of only four or 
eight'1®*) components, it appears doubtful that Dunn’s theory is capable 
of predicting, or even explaining, actual coarsening textures in face 
centred cubic metals, without additional ad hoc hypothesis. 

As a summary, it may be stated that the only situations so far definitely 
established, where the orientation of grains formed on annealing is not 
determined by selective growth, but presumably by the available individual 
nucleus orientations, are the special ones characterized by scarcity of 
nuclei in orientations favourable for growth. In these cases the oriented 
growth theory does not apply, since the prerequisites for selective growth 
are not fulfilled, there being not enough nuclei to select from. It is 
interesting that even in slightly deformed single crystals the orientations 
of the naturally available nuclei are apparently quite random since, in 
the absence of selective growth, no well defined annealing textures are 
formed. On the other hand, after severe deformation, where one might 
expect much more random nucleation, since the local disorientations are 
much greater, in fact sharp annealing textures are formed. This may 
be ascribed most likely to selective growth, made possible by the abund- 
ance of nuclei. The observed predominance of the same kind of 
reorientation, corresponding to [111]-rotations of 30° or 40°, in both 
primary recrystallization and coarsening in all f.c.c. metals naturally 
follows from the oriented growth theory and from the observation. that 
this particular orientation relationship also corresponds to maximum 
boundary mobility. In contrast to this, no oriented nucleation principle 
so far published proved to be suitable for generalization, i.e. for explaining 
more than the one case for which it was originally invented. The close 
similarity or identity of the predominant orientation relationship in a 
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number of carefully documented cases of both coarsening and primary 
recrystallization in f.c.c. metals with that for maximum boundary 
mobility would seem to require of a more general oriented nucleation 
theory that it assume for the nuclei in both primary recrystallization and 
coarsening the same orientations that are known to correspond to 
maximum boundary mobility. An ‘ oriented nucleation ’ theory of this 
kind, as proposed recently by Liu and Hibbard‘18”, is somewhat difficult 
to distinguish from the oriented growth theory, except insofar as it is based 
on hypothetical, and not independently verifiable nucleus orientations, 
instead of on the experimentally observed orientation dependence of 
boundary mobility. In the case of a matrix with a sharp single orientation 
texture it leads to the same predictions as the oriented growth scheme, 
but it can not account without additional ad hoc hypotheses for the 
absence of some of the predicted reorientations in the case of complex 
matrix textures. Also, such an oriented nucleation theory would have 
great difficulty in explaining the extreme sharpness of the cube texture’?°) 
formed on annealing from a rolling texture that includes very considerable 
orientation spreads. If it is assumed that the nuclei have a definite 
erystallographic relation to the matrix, one would expect each nucleus 
to be so related with respect-to that microscopic portion of the matrix 
directly surrounding it ; consequently, the spread of nucleus orientations 
should correspond approximately to the spread of matrix orientations. 
In the theory of Liu and Hibbard ‘’” it is assumed that on recrystallization 
of rolled Cu nuclei are contributed mainly by matrix elements having the 
‘modified (358) [352] ’ ideal orientations of the rolling texture, and to a 
very minor extent by matrix elements with the (110) [112] orientation. 
The continuous spread of matrix orientations, for instance that between 
the two orientations mentioned, is arbitrarily assumed to contribute no 
nuclei. 
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81. INTRODUCTION 


AN electric point charge e brought into a dielectric medium polarizes its 
surroundings. If this point charge is static then the polarization is 
determined by the static dielectric constant. For a moving point charge 
the polarization can only be determined if the dynamic properties of the 
medium are known. If furthermore the point charge is replaced by a 
particle with a certain mass and charge, then the polarization produced 
by the particle reacts on it. A situation arises which can be described 
in terms of a field, the polarization field, which interacts with the particle ; 
the latter then possesses a self energy in the field. A simple and in- 
teresting application of this idea is presented in the motion of a free 
electron in the conduction band of an ionic crystal with which the present 
article is mainly concerned. 3 

The interest in discussing this question is twofold. Firstly, of course. 
a calculation of the properties of free electrons in ionic crystals has an 
intrinsic interest. Secondly, however, this case provides a very simple 
example for a non relativistic field theory, and in view of its simplicity 
it might be expected to lead to the discovery of a number of new features 
of such fields and to the development of new methods. In fact after the 
first application of the methods of field theory to electrons in ionic crystals 
(Frohlich, Pelzer and Zienau 1950), the use of these methods in metals 
led to an important step in the theory of superconductivity and to the 
prediction of the isotope effect (Frohlich 1950). Subsequent discussion 
has shown, however (Fréhlich 1953) that new methods are required to 
deal with all aspects of superconductivity, and at present it seems that 
a method which would combine the results of the two methods described 
in § 4 and §5 might fit the requirements (cf. end of § 6). 

The methods of § 4 and §5 can best be described as dynamic and static 
respectively. The static method is essentially an application of Hartree’s 
self consistent field method. Thus an electron with mass m in an ionic 
lattice may be described by an electronic wave function giving rise to an 
average charge distribution p. Treating this p as a static charge, it 
establishes a polarization in the lattice which can be calculated from 
the laws of electrostatics. This polarization gives rise to an attractive 
force on the electron, and to a potential energy of the order —e?/le* if 
the charge p extends over a spherical range whose radius is of the order / ; 
e* is an effective dielectric constant. The restriction of the electron in 
space requires its de Broglie wave length to be of the order / so that its. 


Z2 
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kinetic energy is of the order 47?/?/2ml?. Minimizing the. total energy 
E, with regard to J leads to 


ra] e2  4Am2h? al em R 1 etm (1.1) 
5(-is — > 7 = Apthie® 2? 8 on a Braue 3) tee e 


A procedure of this kind can be correct only if the time the electron 
requires to travel the distance / is short compared with the period of 
the vibrations of the lattice ; otherwise it is not permissible to treat p as 
a static charge distribution. This condition is not fulfilled im most 
cases of practical interest and it is necessary, therefore, to take account 
of the dynamic properties of the lattice. This is carried out in §4. To 
obtain a qualitative picture it should be noticed that if w/27 denotes the 
frequency of the lattice vibrations giving rise to polarization a point 
charge moving with velocity v polarizes the lattice to its full static value 
only at distances large compared with v/w. Thus for sufficiently large 
distances the polarization potential decreases like a Coulomb field ; for 
distances less than v/w it remains nearly constant. The interaction with 
the point charge is thus of the order —e%w/ve*. For an electron this 
holds only if its velocity is sufficiently large to make its de Broglie wave 
length h/mv shorter than v/w. One may expect, therefore, that the 
lowest value the potential energy can have on the dynamic model is 
of the order —e?/de* where d is the length at which h/mv=v/w, i.e. 
d=(h/mw)/?. In §4 it is shown, in fact, that the self energy on the 
dynamic model is of the order —e?/de*. 

The range of validity of the two methods will be discussed in §6; it 
will also be shown then that the linking up of the two methods has not 
succeeded so far in view of certain difficulties inherent in the Hartree 
method. 

Finally in §7 questions of mobility will be discussed. 


§2. Larrick FreLps 


The field of force acting on an electron brought into an insulating or 
semiconducting crystal can be derived from a periodic field of force. This 
latter will be defined as the field of force in an ideal crystal assuming that 
the atoms or ions of the lattice are undeformable and are rigidly fixed to the 
lattice points. FF. Bloch (1928) has shown that an electron can move freely 
through such a periodic field of force though the connection between its 
energy and its momentum in general differs from that of a free electron in 
vacuum. In many cases, in particular for slow electrons in the conduction 
band, the electronic energy y, is a quadratic function of the momentum 
hk so that 

gre IM yk ry ay eles alee a 


as for free electrons although the effective mass m in general differs from 


the electronic mass m,,. This case (2.1) will always be assumed for the 
following. 
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Such an electron for which the periodic field of force has already been 
taken into account is subject to further forces resulting from deviations of 
the field of force from strict periodicity. These are due to displacements. 
and deformations of atoms or ions (lattice defects will be disregarded here). 
Such displacements can be described by attaching a displacement vector 
to each lattice point. Moreover for electrons with sufficiently low energy 
the electronic wave function changes only very little in one lattice distance. 
Therefore the lattice can be treated as a continuum. The displacement 
vectors then become a continuous function of position, i.e. they form a 
vector field. 

The present article is mainly concerned with isotropic ionic crystals. 
Here displacements leading to electric polarization have a much stronger 
interaction with electrons than other types of displacements. The latter 
will, therefore, be disregarded though their influence on the mobility will 
be mentioned in §7. Let P(r) be the electric polarization at the position r. 
Its interaction with the magnetic field of a slow electron is negligible 
compared with the electrostatic interaction. Therefore the rotational 
part of P is irrelevant so that the condition 


curl P(r)=0, bmi tate Bb cae “ais (Bons xp ey (ao ae) 
can be imposed. Hence if electric displacement D(r) and electric field 
E(r) are introduced in the usual way by 
D(r)=E(r)+47P(r), ee aee ne tae oe (2k) 
the only sources of D are free charges, i.e. in the case of a single free 
electron in the crystal, 
HE eet ciate Fe a? Bs slots (x2) 
Le. curl D=0; div D=47es(r—r,,), 


where r,, is the electronic coordinate and +e the electronic charge. In 
the absence of electrons D—0, therefore,so that D represents the ‘external’ 
field acting on the polarization of the crystal. Since P(r) d?r is the electric 
dipole moment of a volume d’r around r, clearly if a potential P(r) is 
defined by 

de Pirjarad O(6)y4 oe id ee ye 2 (205) 
the interaction energy between electron and lattice displacements is given 
by (using (2.4)) 


— | Dir, ra)P(r) d?r— | ered ieee P(r) d®r=eO@(r,,). . (2.6) 


Since the lattice polarization does not make any contribution to D, it 
- follows that —47P represents the electric field due to the polarization, and 
@(r) is its potential. 
It should be emphasized now that the polarization P is determined by 
the dynamic properties of the lattice. Thus while D is always given by 
(2.4) independent of the state of polarization of the lattice, the electric 
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field E depends on it, and can be found from (2.3) provided P has been 
determined from the dynamic equations. In reasonable approximation 
(cf. Fréhlich 1949, § 18) P can be considered as composed of two contri- 
butions 

P(r)=P,(r)-E P(r) = eo eae 
each of which satisfies a harmonic oscillator equation with frequencies 
i)/27 and w/27 corresponding to optical absorption in the ultra violet 
(P,) and infra red (P,,.) regions respectively. The former is mainly due to 
deformation and the latter due to displacement of ions. In the presence 
of an external field D, 


P..(r)-LwP,(r)=D(r, rj)  - se 2 Sate) 

P,(r)-}a2P,(r)=D(r, r,,)/9, 1g coe aon 

where y and 5 are constants closely connected with the energy due to the 

respective displacements (cf. eqn. (2.20)). The frequency w/27 is corre- 
lated to the frequency w,/27 of residual rays by 

w=lelen) Aas Trde oi Deets ee 


where « is the static dielectric constant and «,,!/? is the refractive index in 
the optical region before onset of ultra violet absorption. 
To find y and 6, consider the static case in which P=9, 1.e. 


4P(r)=(e—1)E(r)=(1—I/e)D(r), . . . . (2.11) 


making use of the static relation D=c«E. Furthermore, the high fre- 
quency dielectric constant ¢,, is defined by D=c.,E on the assumption that 
the frequency w,./27 of the external field satisfies the condition 
Wy>W,>w;, In this case equations (2.8) and (2.9) show immediately 
that the infra red polarization can no longer follow and hence is negligible, 
P,,.~0, but that the optical polarization P, can follow nearly adiabatically 
so that Py can be neglected as compared with w Py. Hence P, has 
nearly the same value as in a static field of equal strength, 


47P)(r)=(1—1/e,,)D(r), (static case) oo See 
so that from (2.11), (2.12) and (2.7), 
1 1 ‘ 
4rP,,,.(r)= (= — :) D(r) (static case). . . . (2.13) 


Comparison of (2.8) and (2.9), for the static case (P,,=0, P,=0). with 
(2.12) and (2.13) respectively leads to 


1 afl 1 

ters a(=- -), (2.14) 
ne wo” 1 ii 

S per eS ( wae ~). . 5 . . * Fy (2.15) 


The equations of motion (2.8) and (2.9) of the lattice field have to be 
supplemented by the equation of motion for the electron which follows 
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immediately with the help of the interaction energy (2.6). Thus classi- 
cally 
Wie pred Dire hs. s,s) (2.16) 


All three equations of motion can be derived from the Lagrangian 


| D(r, r.,)[Po(r)-+P,,(r)] d’r+4mr,2,  . . . (2.17) 
by putting 
SE ONE, weiss Gh PCRS) 
and treating the components of r,,, Po(r) and P,,(r) as variables.t It will 
be noticed that in contrast to most other types of fields the spatial 
derivatives of the functions P,(r) and P,,(r) do not occur in the Lagrangian. 
Therefore in a straightforward way the components of these functions at 
different positions can be taken as different variables and (2.18) then 
yields immediately 


(2.19) 


if g represents any of these variables ; the differentiations with respect to 
them follow by considering the integrals as limits of sums. Equations 
(2.19) then lead to eqns. (2.16), and to (2.8) and (2.9) for all r. 

The moment conjugate to q is defined as 0L/dq so that yP,,(r), mr,, and 
5P,(r) are conjugate to P,,(r), r., and P,(r) respectively. Hence the 
Hamiltonian is given by 


W=25-j-1=h { [P,.2(r) +-w2P,,2(r)] dr 
q 


+43 | [Po%(r)+92P,%(r)] dr 
[Dery rallPnle)+Po(r)] Brtdmig® (2.20) 


In the classical case of an electron at rest, when P, and P;, are given by - 
(2.12) and (2.13), and D by (2.4), H’ becomes infinite because the 
integrands in (2.20) diverge as r—r,,. This arises from the disregard of 
the atomic structure of the substance. In a rough way this can be 
remedied by developing the polarization vectors in a Fourier series, and 
omitting terms with wave lengths shorter than a critical one whose 
wave length is of the order of the lattice distance. 

In the quantum mechanical treatment this divergency does not occur 
as will be seen in §§4 and 5. Nevertheless, for simplicity, a classical 
argument will be used to eliminate the contributions of P, to the 


+ In view of conditions (2.2), Py and P;, actually have only one available 
component each. 
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Hamiltonian. For this purpose it will be noticed that the force acting on 
a classical electron at rest must vanish, so that ®(r) has a minimum at 
r=r,, (provided the above mentioned ‘cut off’ at short waves is used) 
and P(r,,) vanishes. For a moving electron, however, P(r,,) does not 
vanish because of the time dependent terms in the equations of motion for 
P,,(r) and P,(r). Inspection of eqns. (2.8) and (2.9) for a time dependent 
D shows that for slowly moving electrons P,(r,,) remains nearly zero and at 
any rate is much smaller than P,,(r,,) because w»? is very large compared 
with w,,2.~ Physically it means that the deformation of the ions giving 
rise to P, has a much smaller inertia than the displacement of ions (leading 
to P,..). It follows from the above that all terms in H’ containing P,(r) 
depend only very insensitively on the velocity of the electron. For small 
velocity the quantity 


ees | {45[P,2(r) +e 2Po2(r)]—D(r, re,)Po(r)}dr, . . (2.21) 


can then be considered as a constant and the Hamiltonian of a slow 
electron in an ionic lattice field can then be written as 


H=H’—C= | Qy1P,2(r)+o2P,2(r)]—D(r, ra)P,,(r)} drt mig? 


It is useful to introduce here a complex field vector B(r) such that 


B(r)= (55) (Palr)+= Par), 5. Kairie Bats 
B+(r)= (53) (Pal) —=P.(e)), 2. (2.284) 
curl B=curl Bt=0, 5. tok a 
or 
ae h\v2 
wl)= (5) "TBH + BIN, . 2.25) 
p hw\1/2 | e 
ult)= (3) i[B+(r)—B(r)]. . 2... (2.26) 


Then with (2.22), 


Hho | B*(r)B(r) dr— (SA) [DE eay1BH()+B(e)] dr + PE, 


ee 2.27 
where P.=mr,,. The constants in the definition of B have been meee 
in view of the following quantization. For this purpose B(r) will be 
subjected to a periodic boundary condition over a cube of volume V=L?, 
i.e. considering (2.24), 


= exp (wr) Ww, _exp (—iwr 
Bir) = 2 by Be(r) 2b, Pw) P (2.28) 
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where the coefficients w; of the wave vector w satisfy 


We My OM teased tut). (Cae wans.” (2.29) 
Also defining (cf. (2.5)) ®,,(r) by 
ome ey ara. (rh) eh te aye. (200) 
it follows that (cf. (2.6)) the interaction energy is given by 
~ [D¢r, r.:)P,,(r) d’r=e®,,(r,,), ... . . (2.31) 
or using (2.30), (2.25) and (2.28), a 


e®,,(r.,)=4rt (; 


Hence with (2.27), 


a= “oe thw 2b, tb, 


pay Oe ee ; 
pees DyaV arm exp (—7wr,,)—b,, exp (¢wr,,)]. » (2.33) 


It will be noticed that (2.32) holds only if w40, and in fact terms with 
w=0 must be excluded from all the sums. To justify this it should be 
realized that strictly speaking the case of a single electron in an otherwise 
neutral lattice cannot be subjected to a periodic boundary condition. 
For the volume V contains an electron, i.e. the charge e, a situation which 
cannot be continued periodically in view of the long range of repulsive 
Coulomb interaction. To avoid this difficulty an opposite charge, —e, 
should be distributed homogeneously over the volume V. This just 
compensates the component w=0 of the electronic charge, and hence there 
is no w=0 term in the interaction energy. 

For quantization it will be remembered that the components of P,,(r) 
and yP,,(r) are conjugate variables. Making use of (2.25), (2.26), (2.28) 
and (2.2), this means that (h/2yw)'/2(b,++b,) and y(haw/2y)/21(b,,*+—),) 
are conjugate. Thus using the notation [w, v]|=uv—vu, quantization 
according to Bose statistics demands that 


qih[b,,+-+-b,, by+—by|=th, .. » « « « (2.34) 


and that commutators involving different w’s vanish. It then follows 
that 

[Bw Ow t]=Sw, we 3 be ats eee (2.35) 
Algo, in the Hamiltonian, the terms with b,*+b, should by symmetrized 
and hence be replaced by 4(b,*b,+byow +)=b, by +b, +4, using (2.35). 
This means that (2.33) should be augmented by hws which is a constant, 
representing the zero point energy of the oscillating ions. This term can 
be omitted, however, because it would simply change the zero of the total 
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el 


energy. Also, in the usual way, the components pj; of p,, must be re- 
placed by the operator (h/1)0/0r;.;, 1.€. 
[Dyer 7 er]= —?N8;,- EMM Thr MeO ea), 
It is useful at this stage to introduce dimensionless variables, because 
this shows that the problem presented by the Hamiltonian H, eqn. (2.33), 
essentially depends on a single parameter, the interaction parameter 
which is defined by 


2 2 
w= pet =3(=-2) —. i te SS eeeer 
Use has been made of (2.14) ; w is an inverse length defined by 
lak ==N@; 1.6.5 ¢— (=")" (2.38) 
Define now the dimensionless quantities 
X2=Ur a, v=w/u, S=2 Vat, av Seta 


and express the Hamiltonian in terms of iw. Then from (2.33), using 
2.37)-(2.39), 
aed 
ho 
RE Aaro.\ U 
es c6 FT ok tb ti( 3) 2 ou ; lovt exp (—ivx)—b, exp (7vx)]. 
(2.40) 
where x, are the components of x. 


$3. GENERAL PROPERTIES 
The aim of the present section is to discuss some general properties of 
the polarization field P,,(r) or B(r). Consider first the field in the absence 
of an electron. Its Hamiltonian H,j4 follows from (2.33) by omitting 
the electronic and the interaction terms, i.e. 
fe hey oe i H,=hwb,,*b,. . . . . . (3.1) 


Here H,, will be recognized as the Hamiltonian of a harmonic oscillator 
with frequency w/27, coordinate qg, and momentum p,, where 


h 1/2 _{ Mho\ V2 
Iw= (sizz) (bw Bw) 5 Pai ( x") (bw*—by) - (3.2) 


-_ 


and M is its mass. For then 


Pas 
Hy = 7 +4? iat 6. es 88) 


From the theory of the harmonic oscillator it follows that the eigenvalues 
of b,,*b, are positive integers, or zero, i.e. 


OO Fie pat he) MOP ter make! hor) 


‘Sy ua! bee 
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where | 7,,) denotes a normalized eigenfunction. Also it is easily proved 
that (n is a positive integer) 


ies | Oye penne | Vide pend oo ee et (355) 
and 
Ole ate) tO | OL ySs0.. ert. (3.6) 
In view of these properties the operators b,* and b, are denoted as 
creation and annihilation operators, respectively. 

The field is thus represented by a set of oscillators w, all having the 
same frequency w/27. Alternatively the field can be considered as 
represented by polarization quanta of different wave numbers w, but 
all having the same energy fiw, independent of w. The corresponding 
classical polarization waves have thus equal frequency, independent of 
the wave number, a property which follows also directly from eqn. (2.8). 
Hence 

0w/Ow=0, ee ey 8 eS (oy) 


which means that the group velocity of polarization waves vanishes. 

One would conclude from this that thefield cannot transport momentum. 
In field theory the quantity fw is, however, usually denoted as momentum 
of a quantum because in the presence of an electron the operator 


iW hworttetpy es Be 


is a constant of motion, i.e. 
DVVerri == 0 ee SO ae (OOF 


Equation (3.9) follows immediately from (3.8) and (2.33) using the 
commutators (2.35) and (2.36). For clearly p,,2 and b,*b, commute 
both with W whereas 


b +h 9 5 —D r [D. A exp —iwr )J=—hw exp (—iwr,,), etc. 
[ w w w ] w y: el é 
(3.10) 


Hence the commutator of p,, with the interaction term in the Hamiltonian 
just cancels that of Siw, +b,. In spite of (3.9) it will be shown presently 
that a polarization quantum does not carry any momentum. In fact 
through discussions with B. Szigeti and S. Zienau it was realized that 
in any lattice field the quanta cannot carry momentum because the 
individual atoms or ions of the lattice oscillate around their respective 
equilibrium positions, and hence on an average have zero momentum. 
The momentum of the lattice is always carried by its centre of gravity 
as will be demonstrated presently. Hence (3.9) does not represent the 
momentum law but the law of conservation of the total wave vector. 

To show this it will be necessary to introduce the coordinate r, of the 
centre of gravity of the lattice and treat it as a dynamic variable. This 
means that the developments of § 2, in particular the equations of motion, 
the Lagrangian and the Hamiltonian require modifications because in 
§ 2 it was tacitly assumed that the centre of gravity is at rest. To carry 
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this through let r’ be the positional vector in a fixed coordinate system 
relative to which the centre of gravity has the coordinate r,, and the 
electron has the coordinate r,,’, i.e. 


rf. rl ee) MR ee cre RO 
Hence the polarization vector. P,,(r, ¢) has to be replaced by 
P,.(r’—r,, t)=P’(r, 2), 


oP,,,(r, t) oP’ (r’, t) S250 oo 5 
aoe by pares + (tay P (r - t). ie (3.12) 


and 


With these replacements, the new Lagrangian L’ is obtained from the 
field terms of (2.17) (the P,-terms will be omitted) by adding the kinetic 
energies of electron and centre of gravity. Thus introducing the operator 


bee 0 Pees, 
0=d(t, r, r= at (Fog): oEtcete Sy on BN ce (3.13): 
one finds (VM is the mass of the lattice) 
L'= | {dy OP’)? —4yw?P’2+ DP’} d3r’ +imr,,' 244 iNMr, Pee eae 3 


It will be noted that D(r, r,,)=D(r’, r,,’) because this function depends. 
on |r—r,,| only. The variables in L’ are the components of r,, r,,' 
and P’(r’, ¢) the latter for all value of r’. The conjugate momenta are 
given by 


OL’ 
Pa= Fe =NMégty | oP Pid3r’, . . . (3.15) 
Pe! 4 eis 
Pez = Or, =M"r ,, wy i Ss Rael te Ae, eee a ee (3.16). 
and 
pecs oy ee 
= Op’ =yd(t, Le r,)P. Cee ea ok PE 8. 92 (3.17), 


The equations of motion follow with the help of (2.19). Since spatial 
derivatives of P’ now occur under the integral, the Lagrangian derivative 
with regard to P’ has to be taken, i.e. 


Oost 
OJqury OPy = 2 Ory OOP; fon)? 2°) (3.18): 
Hence one finds for g > r.;', with (2.4), 
oD( ; 
mi, | CU ee) awa — (2p ae’ 


ee one he Boy eee IS 
where 


eO'(r.')=— [ Dilr’—ry P(r) dr’, . . (8.20) 
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denotes the interaction energy between field and electron. 
For q~>7,;, (2.19) yields 


NM +y 5 | O(t, r’, tea Pte) or = 0 meee. (0,21) 

and for q > P’, 
0°(¢, r’, r,)P’(r’)+w2P'(r')=D(r’, r,,/\/y. . . . (3.22) 
The latter follows also Seat from (2. ae with the replacement (3.13). 
Clearly if r,,’=r,, then P’=D(|r’—r,,’ |)/yw? is a particular solution, as 


required. The pOlabions of the Pomepentous equation (3.22) (i.e. of 
(3.22) with D=0), can be written as a superposition of plane waves such 
that for one of them 
Poeustexp (wr —tat) .. . . 7 2 (3.23) 
where 
o=eFrw or  dw'/dw=r,, . . . (3.24) 
the latter with the use of (3.7). Thus the group velocity is now equal to 
the velocity of the centre of gravity, as required. 
To discuss (3.21) it will be noticed that the integral occurring in this 
equation can be transformed. For with the use of (3.13), 


0 * 0 
sald , / Lh hp / / 3 
5; | Alt F's Fa)P(r’) PW) de 


Beer i fom \ ion Ya SON OP’ ee EN OE ae 
= || 422” (Fae) (P+ (eae) ’) ar 7+ (P+ (F Car P mm le 


‘ Boles) (325) 
Now through integration in parts it is found that 
/ , c Pp’ 
{P po arr = tite aes ee =0; similarly | Ps oF dr =0 
22’ oP’ 

Te 2 
and or, Or ae Ps 8h sa had ou eres (os §) 
because all surface integrals vanish. For this reason, also 

: a ed 
{ (tay) (P a P 57) dr at Oe erate | 


Thus on the right-hand side of (3.25) all terms vanish except 0?P’. Making 
use of (3.22) and again of a 26), and of curl P’=0, 


, OP’ , , dP’ | coe I g 11 3y/ 
jar = dir =a | Pe ae +2 (Dod =2|(D5z) Piste 


Ly¥oD. , ; 
aoe) ae Pe sede (ee 2) ts a aad aren (8128) 
the latter with the use of (2.4). Thus from (3.21) and the above, 
NMr,=4reP’(r,,'). Peep tral ae: (329) 
Comparison with (3.19) shows that . 
Pager ia ip see EN Mrs)... (8580) 


dt 
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which represents the momentum law on the usual definition of centre of 
gravity momentum by NMr,. This definition differs from that of p, 
by (3.15) ; this is a pure for angie required by the Lagrangian formalism 
in view of the occurrence of r, in the field terms. In eqn. (3.30) all 
properties of polarization waves have been eliminated. This equation 
thus shows that any change of the electronic momentum by interaction 
with polarization waves is compensated by an opposite change in V Mr, 
ie. by the centre of gravity momentum (in the conventional definition) of 
the lattice. 

The new Hamiltonian now follows in the usual way from (3.14) using 
(3.15)-(3.17), 
Ha 255, LN MEE Hi [i 4y[(OP’)2-++w2P’2] DP} dr". 

(3.31) 

Again a complex quantity B’(r’) can be introduced similar to (2.23)—(2.26) 
and (2.28) by 


eh LYONUE (aR ee ee es ee « 


or 
h 1/2 
P(r)= (==) [B’*(r’)+-B(r’)] (3.33) 
“YW 
hw 1/2 : 
QP’ (r’) = (=) iB (r) Bry sent 2 ee 
Introducing this into (3.15), the momentum p, is found to be 
p.=NMr,—Lhwb,*b,’. gl ere see 


Hence from (3.31) using (3.35), (3.16) and (3.32)-(3.34) by a procedure 
similar to that leading to (2.33), 


7 1 ; , , D4 La , 
sje SWI (pat 2Awh,, +b)? + aa rete Os 


oh ViR eR 
+4771 (= Oe sr) Py ap low’ * exp (—iwr,,')—),,’ exp (twr,,’)]. . (3.36) 


Quantization can now be carried out in thé usual way by introduction 
of the commutators 
ee Oe *1=S wy, wii [P' ergs r” ey l= —th, [Pejs Tgl=—m. . (3.37) 
Thus the conservation law for the total wave number holds as in 
(3.8) and (3.9), i.e. 


[W!;-=0 with AW’'=p,,' +Lhwb,,’+b,’. . (3.38) 
In addition, however, p, is also conserved hecatiens 
[Pigg ELT se Ors Le es ae 


This follows because H’ is independent of r,. It must be noted that 


a i po pigs oot See 
Me x 7 lH t)= WM (Pet 2hwby "by, Lage ee (3.40) 
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so that with the use of (3.38), eqn. (3.39) is found to be equivalent with 
[p, A |=0 where p=payNir, o...- (3.41) 


a conservation law already found previously (3.30) for the classical case. 

On the usual definition, p represents the total momentum of the 
system. Its conservation is clearly a law different from the conservation 
of the total wave number W’. 

In all practical cases, VM is so large that the first term in the Hamil- 
tonian can be neglected, i.e. it is sufficient to deal with (2.33) or (2.40). 
The main purpose of the above discussion was to clarify the physical 
meaning of the conservation laws, and to show how momentum is con- 
served even though polarization quanta do not carry momentum. 


§4. Srarionary Soxvutions I 
The present section deals mainly with the derivation of approximate 
wave functions and eigenvalues of the Hamiltonian (2.33) or (2.40). 
The problem is one of finding the self energy of a single electron in the 
polarization field. The simplest method to find solutions is application 
of perturbation theory (cf. Frohlich, Pelzer and Zienau 1950). This 
method holds only for weak interaction, i.e. when «<1. 


4.1. Perturbation Theory 
Consider first the solutions of the Hamiltonian 


3 2 


0 
=a t fp) rant eee Pe eey meat 
Fy Pore slaee By ( 1) 


obtained from (2.40) by omitting the interaction terms proportional to 
the interaction parameter ~. The wave functions are products of a free 
electron wave function, and oscillator wave functions of the type defined 
in (3.4), 


Sg eae 


on ==S- TY exp (Gqx)il | n,))- 252. + (4:2) 


where (cf. (2.38)) 
eee IN TY pag Ke WV ne eo ne, AB) 


is the total wave vector in dimensionless form (cf. (2.39)), and 
DR op =eiie | Wes tre okie eas seta thee (4a) 
similar to (3.4). Clearly q is determined by (4.3) in terms of k and the 


n,’s and need thus not be explicitly included in the parameters defining 
the wave function. The corresponding energy «, (multiples of fa) 


satisfying ) 
pliKane rnp Wat =z ea( Mipatsen Nye ay TK gaa sry 3 + Die 4:6) 
is given by 
plore > My ae \=q?+2Iny=(k—2vn,)?+2n,. . - (4.6) 
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The interaction term (cf. (2.40)) 
Aira 


Demat 
F y= F —F y= (=) 27 [by* exp (—ivx)—b, exp (ivx)] . (4.7) 


will now be treated as perturbation. This can be done by straightforward 
use of perturbation theory provided the zero order energy of the system is 
less than fiw, i.e. the corresponding eigenvalue «, of Fy satisfies <g<1, 
because then for given total wave vector k the system is not degenerate. 
For in this-case no quantum is excited, hence k=q, and k<1 the latter 
to satisfy ¢g<1. Thus 


[ks MOOR SSA" exp (lacing) yaat [OF aia ee 


is the normalized zero order eigenfunction with wave vector k, and 
e(k-22 00 JSR 2 ee ee 


is the corresponding zero order energy. 
Now since according to (3.5) and (3.6) 
by] 0,)=0; Bt]0,)=]1y), . - . - (4.10) 


and in view of the orthogonality of the wave functions 
g-1 | exp (ik—k’ .x) dk=3,,; (ny|my)=Sy, . (4-11) 
S 


the interaction term F’,,, leads to transitions from a state described by 
(4.8) to states in which a single quantum is excited, provided the total 
wave vector remains unaltered. Their wave functions are thus of the 


type 


| key -Oy, ctedjee at Onan )aeO 8 exp (tk-ow e)bg tat ee 

with an energy 
eg(k 02. 1 (kv) ee ae 
‘The non vanishing matrix elements are thus given by ‘ 
Ko. ]y-+- 03 KI Fing| 3 Oy... Ly...) =8(4rra/8)¥21/y, . (4,14) 


t $ 


From second order perturbation theory, therefore, 


; dra 1 
E(k)=hw (12 aa) iy AA 


is the perturbed energy (in ordinary units, hence the factor fiw), and, using 
(4.8), (4.9), (4.12)-(4.14), 


exp (kx) . [Ara\V2 _ byt exp (—ivx) 
¥(k)= —_—_ [ 1—7 | — peel Wiese PET aba 

= (1 ‘G 4 ati) ne 
the corresponding wave function. 


The summation in (4.15) is easily carried out by transforming the sum 
into an integral, 


S 
i a : . av, 2) ae 0 St Ree cate (4.17) 
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which follows from (2.29) and (2.39). It will be observed that no restric- 
tion in the v-integration limits of the type mentioned in § 2 for the classical 


case is required because of the rapid convergence of the integral for v>1. 
Thus 


n-1 i 
E(k)=hw Gat ) Ree ee (4518) 
Also since*for k<1, 
STs eg | 
- SER A ml |e Meta, oth, 19) 
we have 
h? 
E(k)=E(0)+ = (uk)?+..., (4.20) 
where (c.f (2.38)) 
1 l\ , /ma\1/2 
E(0)=—ahw=— (=-2)@ (=) 5 Ct Te ea (4.21) 
is the energy of the lowest state (k=0) and 
m*—m/| (1—a/6)—m(1+«/6), ait ate ta (4.22) 


may be denoted as an effective mass. For wk is the wave number in 
usual units, and hence (Auk)?/2m* would represent the kinetic energy of a 
particle of mass m* and momentum fuk. 

The perturbation method described above is restricted to weak inter- 
action, «<1. It leads to wave functions of type (4.16) in which never 
more than a single quantum can be excited. It will be shown now that the 
range of validity of eqn. (4.21) goes far beyond «=1. 


4.2. Gurari Method 


Gurari has shown (1953) that a variational method can be applied in 
order to derive eqn. (4.21) for H(0). In §5 with a different variational 
method another value for #(0) will be derived which is proportional to «?. 
‘This will permit us to estimate the «-range in which (4.21) holds with 
reasonable accuracy because variational methods never yield energies 
lower than the exact one. Gurari’s results were derived independently 
by Lee and Pines (1952) who use a method by Tomonaga (1947) which is 
closely related to Gurari’s method (cf. also Lee, Low and Pines 1953, and 
Zienau 1953). 

In the general theory of quantum mechanics it is shown that the 
solutions of 


PRPs Coa Se Sens (4.23) 
satisfy the variational principle 
SJ=0 with [-PFro1, age RH ae 
where 
Te | Pere, SU Gee Pane are v2 426) 


P.M, SUPPL.—JULY 1954 2A 
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The integrations are to be extended over all variables. They are under- 
stood in a sufficiently general sense, e.g. if the b,’s are represented by 
matrices then the corresponding sums must be taken. 

The functions Y will now be chosen through a generalization of the wave 
functions in perturbation theory. Instead of using a lear superposition 
of wave functions (all with equal total wave vector k) containing no more 
than one quantum, Y(k) will now be chosen so as to lead to no more than a 
single quantum for each wave vector v ; any total number of quanta may 
be present; however, provided they belong to different wave vectors. 
Again let k be the total wave vector of the system, and define functions 


b(k, v)=M(k, v, €)[1+S-1/2c(k, v, €) exp (—ivx)b,+], . (4.26) 
where 


N(k;v, £)=[1+8-1] c(h, v, €))RV21.  . . . (4.27) 


The c(k, v, €) are numerical constants, and € is the cosine of the angle 


between k and v, 
€=(kv)/kv. a ap a oe ea earn 


The volume S (in reduced units) is supposed to be extended towards 
infinity—after completion of the calculation of H(k) ; its occurrence thus 
defines orders of magnitudes and it will be assumed that the c’s do not 
contain S. Let then 


P(k)=S-\/2 exp (tkx) J7u(k, v)y,.. =~ ~« (4:29) 


and determine the c’s from the variational eqn. (4.24). It will be noticed 
that Y(k) is already normalized. Also by development of the product into 
a sum it is easily seen that all terms occurring are of the form 


[k;...m,...)=S-2 exp [7(k—Lvn,)x] IT|n,), . (4.30) 


thus, according to (4.2) and (4.3), belonging to the same total wave 
vector k. 

In calculating the integral J, (4.25) with the use of (2.40) and (4.29) 
consider first the pure field terms S,b,+b,. It will be noticed that all 
#(k, v) commute, and that b,+b, commutes with all (k, v’)’s (vv’). 
Hence using y from (4.8) 


Jy= | P+ Fb tb P=S-1 J i d3x (0, | +(k, v)bytbyb(k, v) | OY) 
Vv v/S 


=S" y [ No, b &)Lolb,v, €) 28H dxx81 S| c(h, », ) . ah 


Here use has been made of (4.10) and (4.11), and N2~1 has been introduced 
because S— > oo. It must be observed that this is always possible in sums 
of the type occurring in (4.31) and in similar ones occurring below. In 
the wave function Y(k), however, where a product ITN (v, k, €) occurs in 
the denominator replacement of N(v, k, €) by unity is not permissible 
because the number of factors in the product increases with 8. 
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3) 


The interaction term in F reduces in a similar way to a sum when 
eqns. (4.10) and (4.11) are used. For then 


J .=1(417ra/S8)1/2 | + Y [b+ exp (—ivx)—b, exp (ivx)]¥ 


=iS-9/9(4a)¥2 2 | d®x (0, | b+(k, v)[b,+ exp (—ivx) 
—b, exp (ivx) }%(k, v) | 0, > 
=1iS—1(477ra,) 1/2 tha iy VRS) sO v, C)v 3 vk OF So 
Finally the electronic term is slightly more complicated, 
Ow i ov 

—S fees Wa | = 13 i8- ak | x xe IL p*(k, ¥) 
; 3 eee ‘ 0 
x 5 Tk, v)xtec.+8 | x Bx Wk, ¥) )( 5 TT lk, v)) x 

(4.33) 


Here c.c. denotes the conjugate complex of the previous term. The 
second term in (4.33) is obtained by similar methods to the ones used to 
evaluate (4.31) and (4.32). Thus with y from (4.8) and with use of 
(4.11), 


~i§-1k | xt IT b+(k, v) Z Hk, vx 


——i8—k J | oaclltas (en ain 
v/Ss Ox 
as Kk Video ty es... «0 (4,84) 
In a similar way the last term is found, and one obtains 
Jg=h®—28-* B (kev) | c(h, v, £) [P+S-* Dv*|e(k, v, €) 
+8-2 5 J (vv’) | c(k, v, €) |? | of ie AION Ah ee all Cet) 
vv 
Here the factor (vv’) can be simplified because using (4.28) 
vv =vv’' (€é’-+sin (k, v) sin (k, v’) cos (6—¢’)) . . (4.36) 


where ¢ and ¢’ are the polar angles of v and v’ respectively. In the 
following summation (integration) over v, the term with cos (¢—¢’) 
vanishes because c(k, v, €) is independent of ¢, and 


{cos (d—¢’) dd=0. 


Thus vv’ can be replaced by w’fé’. Hence using (4.31), (4.82) and 
(4.35), 


J=J AJ gt Jg=h?+S8-12(1—2kv+0?)c*(h, v, E)e(k, v, €) 
- 45 (dra) 28-15 (c*(k, v, €)—e(k, v, €))/vtg? . . . - (4.37) 


2Z2AZ 


ae 
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where 


g—S—Doke*(ky v, S)eky vs é): 15 See 

The variational principle 6J=9 requires that differentiation of J with 
regard to all c(k, v, €) and c*(k, v, €) should vanish. Hence 

0=(1—2kvé+2gvé+v?)c(k, v, E)+0(4ra)V2/v. . . (4.39) 


It is useful to introduce here a vector g whose magnitude is g and whose 
direction is the same as that of k, i.e. 


g=gkj/k.... 29s ah) 223) 
Then from (4.39) 
cee i(4nra) 2 7 
c(h, v, 2 kee) (4.41) 


and c*(k, v, €) is its conjugate complex. 

The above solution holds only if c(k, v, €) does not become infinite for 
and value of v. To prevent this the condition 

(k—g)?<1 25 wutey pada oe epee 

must be satisfied. 

With the values (4.41) for the coefficients c the integral J, (4.25), 
represents the energy (in terms of fiw). Hence from (4.37) and (4.41), 
using (4.38) and (4.40), 


E(k)=hw {i494 4naSAD 


v 


1—2kv-+v? 
foes - V)+0?]? 


2 
Sete 
ahi 92 dra 3D eh . (4.48) 


v P12 k—g. v) pe] 
and from (4.38) and (4.41), using (4.40), 


es (vk) 
Ie kg OE ee 


With (4.17) after integration one obtains 


Roce: 
Blk) =f fg? EO) ee) 


and 
sin“? (k—g) l 
ly eee ee EE a es ee a eS ee ; 
M ot { (k—g)? te aa . (4.46) 
i.e. two simultaneous equations for H(k) and g. Again condition (4. 42) 


must be observed. It will be noticed that eqn. (4.46) can also be 
written as 


d sin-! (k— 
en (k—q) 


9= 2% 7 kg oe ee 
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Now if |k—g|<l1, by development of [sin-1(k—g)]/(k—g) (ef. (4.19)), 
eqn. (4.47) becomes 
ka/6 
1+«/6° 
Thus as k-+>0, also g>0. With the use of (4.48), and of a further 
development of the sin“, finally 


g=a(k—g)/6, or k=(k—g)(1+a/6), or g= 


(4.48) 


E(k)=F(0)+-h?(ku)2/2m*, be <1) See (4.49) 
where 
E(0)=—ahw=— (= — -) Ee (se) (4.50) 
(Sea G 2h : Fang ; 
and 
m*=m(1+-«/6). Set dary teke, eee ato ky 


Thus £(0) is identical with the result (4.21) of perturbation theory, but 
it is no longer restricted to «<1. The limit of « up to which the result. 
is reliable will be studied in §5 with the use of a variational method 
which gives the energy proportional to «2. Also m* is seen to agree 
with the result which perturbation theory yields for «<1. Again this. 
value now holds for larger «. 

Expression (4.49) had to be restricted to sufficiently small values. 
of k to make use of the development of sin~!(k—g). This would not 
restrict, however, the validity of eqns. (4.45) and (4.46) which have 
solutions for all k. In particular E(k) oc khw as k > o. Such solutions, 
however, have no physical meaning because H(k) has to be restricted 
to values for which 

Ek) —E(O)<fuer (4552) 


which means that the energy of the system should be less than one 
quantum above the energy of the ground state. If this were not fulfilled 
then the trial wave functions (4.29) (with (4.26)) would not be the only 
ones leading to the energy H(k) for given k. Others would exist as well 
in which at least one of the #%(k, v)’s would be replaced by a function 
in which the state with one quantum is principally excited (and not 
the state without quantum as in (4.29)). A degeneracy would thus 
occur. It would require a new approach which has not been found yet. 

Equation (4.49) suggests that the electric current carried by the system 
is given by ej=ehku/m* where j represents the average electronic velocity. 
In fact with the use of (4.29), (4.26), (4.41), (4.48) and of the orthogonality 
relations, = 


= ehu OP eh Ai v 
ehuk chu 

rr — — er . . ° . ° ° ° . . 4,53 
mE =m (4.53) 


as suggested. Use has been made of (4.44) and of the fact that the 
direction of the vector represented by the sum must be k/k. In the last 
step (4.51) has been applied. 
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It is also of interest to calculate the polarization near an electron. 
For this purpose it is best to derive an appropriate expectation value 
of the polarization potential ®,,(r) defined in (2.30). Clearly with (2.25) 
and (2.28), using (2.37)—(2.39) 

e®,,.(y)=i(47ra/S)*/*hiw(by+ exp (—ivy)—b, exp (ivy))/v, . (4.54) 
where . 
y=an i, <a C2 


has been introduced. Let now the integration over all variables of the 
wave function Y be split into an electronic one, denoted by Jd*x, and 
one, denoted by the symbol ¢ ), referring to the field variables. 
Thus the integration over all variables introduced in (4.25) can be written 
as 

[=| a@ex< Wr es 
It follows then that (Y+¥Y) is the probability of finding the electron at a 
position x irrespective of the polarization, and « ¥*@,,(y)¥%) represents 
the expectation value of ©,,(y) at a position y and simultaneously finding 


the electron at x. 
Define now 


B,,(x, y= (PG, ly) P| (PP )=S (PD, fy), (4.57) 
which is the expectation value of the polarization potential at a position y, 
provided the electron is at x. For this purpose it will be noticed that 
from (4.29) and (4.26), with the use of (3.5), (4.10) and (4.11), 
(Pb WY =S-4 (y+ IT p(k, v)by* I bik, v)x) 
=S-1(0,[1+S8-/2c*(k, v’, €’) exp (iv’x)by] 
xX by t[1+S-2e(k, v’, €) exp (—iv’x)b,-*]0, > 
=S-*/2e*(k, v’, &) exp (tv'x). a ay en ae 
Hence from (4.54), (4.57) and (4.41) 
e®,,(x, y)=t(47ra) 1/28 iw XY [c*(k, v, €) exp (iv .x—y) 
—c(k, v, €) exp (—iv .x—y)]/v 
2 cos (v . k—y) 
=—4 ee MAE EE 
ma akg voy 
In the simplest case k=0 (i.e. also g=0) one has with (4.17), 
2 cos (Vv. x—y) 1 cos (v| x—y |») 
gs-1y Y |”) 
tos i.e fi, Saree Orr 
sin (v | x—y |) 
~ 7x alle v(1-+v?) dv 
__ |—exp (—| x—y ]) 
7 | x—y | 


(4.59) - 


sights ieeme nt tO) 
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Returning now to ordinary units with the help of (2.37)-(2.39) (denoting 
then ®,,(x, y) by ®,,(r,,, r)) eqns. (4.59) and (4.60) yield 


Lt a 
e®,,(r.,, r)=— (=- -) Trg] [l1—exp (—u|r—r,,|)]. (4.61) 
el 


cO 


Thus if |[r—r,,| > L/w, the polarization potential of a point charge is 
found. If lr—r,, | < 1/u, however, one finds 


1 1 
US Cg oe mr pe e2u, lr—r,,| < l/u, (4.62) 
i.e. the polarization potential of a charge extended over a sphere whose 
radius is of the order l/u. The case that k is different from zero leads to a 
slight modification of (4.61). 


§5. SraTIONARY Souutions IT 


The eigenvalues of the Hamiltonian (2.40) derived in the previous 
section are linear in the interaction parameter «. In the present section 
a different method, originally due to Landau (1933) (cf, also Pekar 1946, 
1949) will be used which has been discussed qualitatively in §1. It will 
be shown that it can be derived with the help of the variational principle 
(4.24) and will be found to yield eigenvalues which are quadratic in «. 

For this purpose use the trial wave function V as a product of an elec- 
tronic wave function 2(x) and a lattice wave function ®, 


ee OID arte ee wees cee as 50) (5.1) 


where {2(x) depends on the electronic coordinate x only, and ® on the 
coordinates describing the polarization waves. It will be assumed that 
both are normalized, i.e. using the definition (4.56) for the respective 
integrations 


| _d%x2+09).2%)=1 ; (G+P)al 2. (8:2) 
S 


Clearly with (4.1) and (4.56), 


3 02 
[verw=— 2 | D+x) a5 
1/8 v; 


v 


Q(x) d’x-+_ F (O+b,4b,8), (5.8) 
and with (4.7) 
[PFige?=i(4era]S)12 (D+ { d9xQ+(x) 5 [b,+ exp (—ivx) 
S v 


—b, exp. (ivx)]Q(x)@). 5 oe ws wes (5.4) 


The conditions (5.2) can be satisfied with the use of the Lagrangian 
factors 4, and A, so that the variational principle (4.24) requires 


{| v+FYA, | 2+2 dr, (O+0)} =o, ee es) 
Oe 
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or using (5.3) and (5.4), 
- ay 
5 —| dxa+(x) (3 42,) Q[x)-+3¢G+(Sb,tby—Aq)® ) 
s 1 0%,” v v 
+(or| dx Q+(x) F ing Q(X) x)®)b=0, eee ee ee 
S 


where the functions 2, ® and their conjugates have to be varied. This is 

simplified because only the last term depends on both. Thus variation 

with respect to +(x) yields 

pe! 
0x, 


Similarly from variation of @*, 


—2 ) 4 (G+F DB) Ax)=A, QO")... . (BT) 


Zb,+b,®+ | d®xQ+(x)Fy,QG=r,P. . . . (5.8) 
v S 


Variation of 2 and @ leads to the conjugate equaticns. Equation (5.7) 
is the Schrédinger equation for an electron moving in a ‘ self consistent ” 
field «®+F,,,®). Similarly (5.8) represents the motion of the field of 
polarization waves in a ‘self consistent’ source. Clearly in view of 
(5.4), ® can be written as a product, 


®=]1D,; also Ay=ZAy . de is ee 


where ®, depends on the coordinates of the v-th field component only. 
It satisfies 


b, +b, D, +4 (4era,/S)¥/2(1/v)(b,+ | Q+(x) exp (—ivx) Q(x) d3x 
S 


v 


—by | 200 exp (tvx) Q(x) d§x)®,=A,D, . . (5.10) 


which is the equation of a displaced harmonic oscillator. Hence 
( B,+b,+®, ) =i (4rra/S)¥/%(1/v) | Q+(x) exp (ivx) Q(x) d8x. . (5.11) 
Ss 


With the use of (5.4) and (5.11), therefore, the self consistent potential 
energy of the electron becomes (use (4.17)) 


« D*Fing® ) 


=— ale dx’ Q+(x') Q(x’) [exp (iv . x’ —x)-+exp (—iv . x’—x)] 


dara 2 cos (v . k’— 
eos se OE a BHx) 20°) 


=-2 J = Sa 
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It represents the energy of the charge e at x placed into the polarization. 
potential due to the charge distribution Q*(x’) Q(x’). This is easily 
seen by going back to conventional units with the help of (2.37)— (2.39), i.e. 
1h oe yl Q+(r)Q(r 
hw <{(O+F, 4.2 >= — (= —2)e ay = ate. see (Bard) 
Eco : | Ma ter | 
Introduction of (5.13) into (5.7) leads to an integral equation for: 
Q(x). Qualitatively, it is clear that Q(x) must be a function which is 
large in a limited space only. For only in this way can the ‘ self consis- 
tent ’ potential energy have relatively large negative values. To find 
an approximative solution use will be made again of the variational method 
in its original form (4.24). For this purpose let 
(x) =(B3/87)1/2 exp (—B|x—x,|/2), . . . (5.14): 
where xX, is an arbitrary radius vector. Then with 
ys (| d®xD+Q+(x)FQ(x)®), . . . . (5.15) 
Ss 
the variational method requires that 
et ae Re Weg gre tree yen iy, (0-1 OY 
provided ® and 2 are normalized, and @ is a solution of (5.8). As 
mentioned it is composed of displaced harmonic oscillator solutions. 


Now using the expression for F according to (5.3) and (5.4), together 
with (5.8)-(5.10), the integral J, ae in (5.15) becomes 


ee I. Qt was 5 Ox Silex tes ga ssh (B,17} 

Also from (5.10) for the lowest states of ie displaced harmonic oscillators, 
| 

y2 


=—(4ra)/8) { dx} 24x) 2x )aZ(1/0*) | ax 24x) 2 (x’) 


DAV= — (4774/8) +(x) exp (—tvx) 2(x) dx 


[exp (iv . x’/—x)+exp Binet 
abi | Qt) Q(x) (OFF, 0) dix, ww es (6.18) 
2) 5 


where eqn. (5.12) has been used in the last step. Thus the potential 
energy SA, appearing in (5.17) is obtained by averaging the potential 
energy (5. 12) over the density distribution 2*(x)2(x) of the electron, 
and multiplying it by 4. This reduction by a factor } arises from the 
positive contributions in displacing the oscillators. 

Introducing now (x) from (5.14), one finds with (5.18) and (5.12), 


— a= exp (—B | x—X, |) 
Zdy=—a(B9/87)? Jae f ane DBL He =a 0 


Pan ee ee ee (519) 
— 7 xP: ° e o e « e ° 
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Furthermore, 


ee 
ssi QM) Ey Ox) dx=38% (6.20) 
1 

Hence with (5.17) 


5 i z 
FB eta pe aA ae = 521 
j= 16 %B + GP. eee (5:21) 


It will be noticed that the wave functions used here are all normalized 
so that from (5.16), 

e dh j= Psat: (5.22) 
p= 3%, and hence od a eg sae ee 
In ordinary units, therefore, with the help of (2.37)-(2.39), 


25 25 Al  1\%e*m 
ze ee peat ce | 5.26 
E(0)=hod = — se hw=— =, ;(= *) yz et) (8-28) 


The polarization potential ®,,(y) at a position y defined by (2.30) and 
(4.54) is an operator whose essential part is now, however, a c-number. 
For from (5.10) and (5.11), 


by += i(4ma./S)¥/2(1/v) | exp (ivx) Qt(x) Q(x) dk-+d,+, . (5.24) 
S 


where d, and d,+ are operators satisfying a harmonic oscillator equation 
a, tay | ven, | ayy Whee in, Pe ee ee ee 
Hence from (4.54) and (5.24) 
0 y(y) = —(4ra] Shen (1]0) | d®x exp (iv . x—y) 2*(x) Q(x) 
+¢.c.-e@,,’ c : 3 eee rar ae ma PS 


where @,,’ depends on the d,’s in such a way that its expectation value 
vanishes. The main part is a c-number and is equal to expression (5.13) 
(cf. (5.12), replacing there x by y and x’ by x). Hence @,,(y) is the 
polarization potential of the static charge distribution +(x) Q(x), and 
the polarization P,;, can similarly be derived from the electric dis- 
placement D due to the charge distribution Q+(x)Q(x). Thus from 
(5.26) is is found that (c.c. leads to a factor two) 

3 o2 O 

23 Pil) = — 4rrahw Q*(x) Q(x), ee Ks) 

1 OY; e 


or in ordinary units using (2.37)-(2.39), 


1 1 
V2®,(e)—tnp(lr—ro)) (==), ete see Poca 


e 
where 


pll ry eu? (S — 2) 24] xx) 2 xo) . . (6.28) 


is the electric charge density. Spherical symmetry of both ®,, and p 
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around the point ry=x,/w has been indicated. Clearly using (2.30) and 
div D=47p(|r—ry|). . . . . . . (5.30) 
one finds with the appropriate boundary condition 
P,.=D/yw? SN ce ar a 
using the definition (2.14). 

This shows that the treatment of the polarization field in the present 
section is essentially a classical one based on the static charge dis- 
tribution p. In fact in this case (5.31) follows immediately from (2.8) 
when P,,=0. 

This classical approach can be used to estimate the effective mass 
of the system in the present model (Landau and Pekar 1948, Pelzer 1950). 
For this purpose the electron will be considered as a rigid charge 
distribution moving with a small velocity. This means that ry is now a 
point moving through the field, Le. 

er Matra, (5.32) 
where r, is the constant velocity. Since now P,,. no longer vanishes, 
(5.31) has to be replaced by (2.8), or (5.28) by 

some ibe 2 
V2(®;,/@?+ ®,,.) =47p(|r—rot|) (=-*). ae ey D0) 
Now if r, is sufficiently small then 
P= Od heorgiy Ds oo ye a. oe © 2 (5-84) 
where ©, satisfies (5.28) with ry from (5.32), and |6,|<|®,|. With 
this latter condition, and with 


Girt i)—2, "55.2, =e, for, .-., . .(5.35) 
eqn. (5.33) leads immediately to 
BP (ro 2fen2) Dy’. es, ee (B88) 


In this perturbation treatment the total energy of the system is equal 
to the energy in the static case, (5.23), augmented by terms proportional 
to ro2. These terms can be obtained from the classical Hamiltonian 
(2.22). Clearly from (2.30) and (5.34), ; 


4rP,,=grad Dy—(rp2/w?) grad Dy”, Bee MEST) 
and using (5.35) : 
47P,,—=F, grad Dy —(ro3/w?) grad @y"””. +. =. ~«~. (5.38) 
‘Thus since with (5.30) and (2.30) 
a | DP,, d*r= | (OE gh (dy deeb gee CRC) 


the terms of (2.22) which are proportional to r,? are given by (definition 
-of m*) 
dm* 2 = Yni,?+Hy/47?)to?| {(grad ®,’)?—2 grad ®, grad ®,’’} 

—(Fp2/w2) fp By”. . . (5,40) 
A term 4mr,? for the additional kinetic energy of the electron has been 
added. 
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Both p and ®, depend on | ro— rot | only, and the integration is straight- 
forward. Integration by parts and use of (2.14) and (5.28) yields 
Am*=4m-+4(47/?)(1/eco— Ve) | p2(r) dr! (ee GA 
From (5.29), (5.14) 
i p7(r) dv =4re®u9(B3/87)?| “exp (—2Bx) x? da=e7u383/647, (5.42) 


* 27,3Q3 
was “3 (=-:): Me eek: 


Smw" \E, € 


so that 


m 
Introducing f from (5.22), and using (2.37)—(2.39), finally 
m*)}m=1-.4(5/8)8a*. 2 ee 


$6. Discussion oF SOLUTIONS AND OF METHODS 


6.1. Results 


The solutions for the lowest state of an electron in a polar crystal having 
a self energy H(0) were derived in §§ 4 and 5 with the help of the varia- 
tional method. Since the energy values obtained in this method are 
above the exact value (or equal to it) it follows that the lower of the two. 
solutions (4.50) and (5.23) is the better one. Hence 


il 1 Mw \ 1/2 256 
have aes [aha : a a AY 
A afew (- “ye CH) if a< pA (6.1) 
and 
25 25 f1 1\e%m 256 
ee ae ela Py, ee ek ce ae: 2), 
LD HTT NET (= =) Fie DS aee aes 


In order to see which of the two solutions holds in practical cases. 
consider sodium chloride (NaCl). Here approximately «=5-6 and 
€,.—=2:3. The wavelength of the residual rays is 61 « and the correspond- 
ing frequency is w,/27=4X1018 sect. Hence with eqn. (2.10), 
w>4-9x 1013 sect, Using for m the electronic mass, m=m,,, eqns. (2.37) 
and (2.38) lead to «26, w10’ em, hw0-03 ev. It follows then that 
(6.1) applies to this case leading to a self energy #(0)0-18 ev and to an 
effective mass (cf. (4.51)) of m*=m (1+«/6)=2m. It will be noticed 
from (2.37) and (2.38) that «ocl/\/w. Assuming m=m,,, this is the only 
factor which may vary substantially from substance to substance. In 
NaCl, the frequency w is rather small so that for most substances smaller: 
g-values than «=6 should be expected, i.e. in general (6.1) should hold 
and m*/m should be nearly of order unity. The properties of an electron 
in such a substance are then not changed substantially through its inter- 
action with the polarization field, at any rate as long as its energy is less. 
than hw above the ground state (cf. (4.52)). It is also seen that 
u-'~10~ em is a length sufficiently larger than the lattice distance which. 
justifies the approximation of the lattice as a continuum. 
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It is not easy to find actual substances for which « is so large as to make 
(6.2) the better solution. The possibility of large « might occur in con- 
duction in impurity bands where m/m,, may be rather large, since 
% oc(m/m,,)°. Values of the effective mass of the order 100m,, are 
indicated by the work of Breckenridge and Hosler (1953) on titanium 
dioxide semiconductors, but even there it will be seen that the method of 
§4 gives better results. In this substance roughly «~110, «,~8, wr 
10% sec~* and hence w~1-7 X10? (m/m,,)/% em—! and «a~2(m/m,,)¥2; 
the estimate is very crude in particular since TiO, is not isotropic. To 
make here «=256/25~10 requires mjm,,~25. Hence from (4.51), 
m*/m=1-+-«/6=3, or m*/m,,~75 which is of the correct order of magnitude. 
The method of § 4 might, therefore, in this case just have reached its limit. 
To apply the method of §5 would require «> 10, i.e. m/m,,— 25. In this 
case according to (5.44), m*/m=1+4(5/8)3«4=1200, and hence 
m*/m.,<30 000, a value which is much too large. 


6.2. Method 

Clearly it seems desirable to develop a method which leads to a contin- 
uous transition between the two methods, i.e. in which the energy H(0) 
is of the form —hw(c,x+c,«?). For this purpose a variational method 
which employs a combination of the wave functions used in §§4 and 5 
might be successful. Buckingham (1954) has in fact carried through such 
a calculation. The transition has been found, however, not to be contin- 
uous ; up to a value «, which is slightly smaller than 256/25 it was found 
that the results of § 4 are exactly reproduced. At «=a, there is an abrupt 
change in wave function and a discontinuity of 0H(0)/d«. At a>a 
the results obtained in §5 are slightly modified. These unsatisfactory 
features are I think largely due to some unsatisfactory properties of the 
method used in § 5 connected with the concentration of the wave functions 
(5.1), (5.14) around an arbitrary point x, in space. The wave functions 
(4.29), (4.26) of §4 do not have this disadvantage. For on replacing 
there x in the wave function by x—x,, where x, is a constant vector, the 
wave function (4.29) is simply multiplied by a phase factor exp (7kx,). 
To show this it must be proved that the functions (k, v) given by (4.26) 
are not changed when replacing there x by x—x, (without altering x in 
the Hamiltonian F). This means that in (4.26) c(k, v, €) exp (—ivx) 
must be replaced by c,(k, v, €) exp (—iv .x—x,) where the new coeffi- 
cients c,(k, v, €) must be determined by the variational procedure. It is 
then found immediately that 


o(k, v, €)=ey(k, v, €) exp (ivx,), ©... (6.3) 


so that #(k, v) remains unaltered, as required. 
The property of the total wave function ¥, discussed above, namely 


W(x—x,)=exp (—?kx,)P(x), . . . - - (6.4) 


where k is the total wave vector can also be derived directly from the 
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Hamiltonian F as has been shown by Zienau (1954). It is closely con- 
vectéd with the conservation of the total wave vector. The wave functions 
used in §5 do not satisfy this requirement. Nevertheless they can of 
course be used as trial wave functions in a variational procedure, and in 
fact as has been shown, for very strong interaction («> 256/25) they lead 
to better energies than the wave functions used in §4. One might think 
of modifying the trial wave functions used in §5 so as to satisfy the 
requirement (6.4). This is not possible, however, without completely 
changing the method of the self consistent field used in §5. The difficulty 
is closely connected with the fact that in this method the total wave 
vector is not on principal axes. Thus with its definition by (3.8), using 
the wave function (5.1), 


(-i5 aa vb,*b,) 2(x) PARQ) PD.) See) 
where k is a constant vector ; with the wave function (4.29), however, 
(-15, + Svby"by) exp (tkx)JTb(k, v)y=k exp (tkx)ITb(k, v)xy, . (6.6) 


as has been demonstrated with the help of (4.30). 

The difficulty (6.5) encountered in the self consistent method in $5 
arises always in the use of the Hartree method, e.g. when solving the 
problem of n interacting particles in terms of a product of x wave functions 
each depending on the coordinates of a single particle (or an appropriately 
symmetrized combination of such wave functions). This can be demon- 
strated very simply in the case of two particles, e.g. a proton (coordinates 
x,) and an electron (coordinates x,) forming a hydrogen atom. The 
exact solution is of the form 


p=exp (tkx,)o(x,—x,)> < : Sea 
where x, represents the centre of gravity coordinates, and (x,—x,) the 
relative coordinates. The total momentum operator is given by —ihd/éx, 
(=h x total wave vector) and is on principal axes. In the ground state 
$(x,—x,) is a function which is large only near x,—x,=—0. It means 
that the electron is found always near the proton, but this may be any- 
where in the total space available. 
In the Hartree approximation, on the other hand, one uses a trial 


functi 
unctioD Peay (Xa )y Y) «ES Ee eee es) 


where the functions %,(x,) and %,(x,) are to be determined from the 
variational method. This leads to a Schrédinger equation each for the 
electron and for the proton. The respective potential energies are obtained 
from the potentials of the average charge distributions of the proton and 
the electron, respectively. In order to obtain a satisfactory result for the 
ground state, both 4, and %, must be spherically symmetrical functions 
centred at an arbitrary point x, in space and being large in the neighbour- 
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hood of xy only (%, in a radius of the order of the Bohr radius, ~, Ma 
correspondingly smaller radius). The energy of the ground state obtained 
in this way is a very good approximation to the correct value. The wave 
function of the whole system is, however, concentrated around an arbitrary 
point x). In the exact solution the function $(x,—x,,) (6.7) is of course 
also concentrated around the point x,—x,=0; the function exp (tkx,,) 
provides, however, for the atom to be with equal probability at any point 
in space. In the Hartree wave function (6.8) it would not be permissible 
to multiply yy by exp (ikx,) and to use x, as a further available variable 
coordinate vector because all six degrees of freedom have already been used 
in %,(x,) and w,(x,). It is seen clearly now how this is connected with the 
fact that 0/éx, is now not on principal axes. This condition would restrict 
the relative motion of the electron and the proton just in such a way as to 
make three degrees of freedom available for the motion of the centre of 
gravity. In the Hartree approximation, however, there is no correlation 
between the motion of the proton and the electron, except on an average, 
binding both to the same point Xp. 

The wave functions used in §5 correspond to the case discussed above. 
Thus both, electron wave functions (x) and polarization potential @,,. are 
centred around the same point x, as follows from eqns. (5.14) and (5.28). 

A further difference in the use of the methods of §§ 4 and 5 can be seen 
in the behaviour of the polarization potential ®;... In the Hartree method 
(cf. (5.28)) ®,, is essentially the classical potential of a static charge p. 
The value obtained in § 4, however, depends on the dynamic properties of 
the lattice (cf. (4.61), (4.62)). As a consequence the value of H(0) accord- 
ing to the method of § 4, i.e. (6.1), depends on the dynamic properties of the 
field, i.e. on w, whereas (6.2) found from the method of § 5 is independent 
of w. 

A method of assessing the range of validity of the method used in § 4 
has been developed by Lee, Low and Pines (1953), and independently by 
Zienau (1953). Essentially these authors find a Hamiltonian which is 
satisfied exactly by the wave functions (4.29) so that one can investigate 
the influence of the terms which have been neglected. Gurari (1953) had 
already indicated this possibility. The method applied is to transform 
the variables (operators) O occurring in the Hamiltonian F' (2.40) in such 
a way, 


AE OS ie SF en se (6,9) 


that in the new variables O’ the Hamiltonian F does no longer contain 
terms which are linear in the field variables. In (6.9) the operator S, 


must be unitary. 45 
It is of advantage to carry through the transformation in two stages, 


S,=8S,S8,, and choose 


S,=exp (—idvxbytbhy). . . - - + + (6.10) 
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This transforms the operator of the electronic wave vector, —7id/dx, into 
the total wave vector for 


0 F 
~i8, 2 Sypt=—i exp Ee yy x CxP ESE, 


Oa ; 
se aes LM a) er rer ee ANN 
By equating it to the constant vector k it has been indicated that it is 


a constant of motion as has been shown in § 3. Also 


b, =S,b,S,-1=exp (—ivxb,+b,)b, exp (ivxb,+b,)=b, exp (vx), . (6.12) 
because using (2.35) (n is a positive integer), 
bb, tb,= (byt = 1)0y; b,(b,*6,)"*= Ofte 1)"b,, . (6.13) 


and similarly for b,+. Hence b,+b,=0,'+b,’. Thus inserting from (6.11) 
and (6.12) into expression (2.40) for F, 
F=(k—Lvb,'tb,’)? +20,’ +b, +74(47ra/S)1/22(b,'+—b,’)/v. . (6.14) 


In this Hamiltonian the total wave vector of the system is k, but k may 
have any value. In this way the number of available variables has been 
reduced by three in a similar way as the number of available variables of 
a system of interacting particles can be reduced by three by eliminating 
the centre of gravity motion. 
Now choose ' 
S,=exp [—Ze(k, v, €)(6t-+- By]. os se ee AGE) 


and let the (imaginary) constants c(k, v, €) be given by (4.41). Hence 
By=S5b,'Sp-*=exp [e(k, v, €)(by'*+by')]by’ exp [o*(h, v. €)(by’ +6,'*)] 


ale 68h, 0, S)) ae 5! ay te ee 
since 


c(k, v, €)=—e*(k, v, €) and [b',, 6,’ +*)=nb,+*-1, . (6.17) 
It will be seen that with the choice (4.41) for c(k, v, €) the linear terms 
in f, just cancel when (6.16) is inserted into (6.14). Making, therefore, 
use of the definition (4.43), (4.44) of H(k), it is found that 


F=(E(k)/hw)—k?+ (k—LvB,*B,)?+F,+Fo, «satel Oakey 
where 


geek aight - Volk, v, E)e*(k, 0", €')(2By*By—By*By *—ByBy). (6.19) 
and 
Peale reer Ae v, £)By By *By +e*(k, v, E)By *ByBy)- (6.20) 


‘The wave function now depends on the variables 8,, By+. In particular, 


it is seen that a product of oscillator wave functions | n, ) for the oscillators 
8ytB, satisfies 


{H(k)/fico—k?+(k—EvB,*B,)2} IT | ny) 
= {E(k)/heo=k?+ (k—Zvny)*}ITjn,), . . . . (6.21) 
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and is, therefore, a solution for the Hamiltonian F—F,—F,. Hence 
the wave function x, cf. (4.8), satisfies 


(F—F,—F,)y=(E(k)/iw)y . 2... (6.22) 


i.e. the solution of § 4 for a given k is identical with the lowest solution 
of F—F,—F,. 

To estimate the error committed in this approximation one can treat 
P,+F, as a perturbation. This will not be carried out here in detail. 
By second order perturbation including some numerical work, Lee, Low 
and Pines (1953) find a correction 


AB = —(0-00702+-0-01k2x2/(1+-c/6)2)\iow. =. . . (6.23) 


Thus for k=0, one finds «? terms as in § 5, though with a much smaller 
factor. This means that when «256/25, perturbation theory for 
F,+F, is no longer valid. The range of validity of perturbation theory 
for F,+-F, should in fact be restricted to «-values for which 


iitje<tusy (1-8). a®<160, = .: 6 4. x (6.24) 


At 2=256/25=10, «?<100 which is not very small compared with 160. 
For «=6, however the value for NaCl, «?=36, the correction (6.23) is 
sufficiently smaller than iw to make the results of § 4 reliable. 

Equation (6.18) also shows that the case H(k)—H(0)>hw leads to 
difficulties because then the solutions of F—F,—F, are then degenerate. 

In conclusion, therefore, it has been found that the method of § 4 gives 
reliable solutions for energy and effective mass if the interaction para- 
meter « satisfies «?<100 and provided H(k)—H(0)<hw. There is no 
obvious way in which this method can be extended beyond these limits. 
The method of § 5 should give reasonable energies H(0) for «?>100 though 
in this case it may no longer be justifiable to treat the lattice as a con- 
tinuum. This method leads to difficulties inherent in the Hartree method 
and connected with the motion of the centre of gravity of the system. 
The validity of the calculation of the effective mass in this methed is, 
therefore, unclear. 

It should finally be remarked that application of the method of § 5 
to metals has led to an understanding of the superconductive properties 
of a one dimensional superconductor (Frohlich 1954) but does not yield 
the isotope effect which is obtained, however, by methods corresponding 
to those used in §4. This agrees with the observation that the method 
of §5 does not introduce the dynamic properties of the lattice. Thus 
in this case too development of a method linking the results of §4 and 
§5 would be very desirable. 


$7. Mopinity 


If a weak external electric field is applied to the system consisting of an 
electron and the polarization field then an electric current is established. 
Its magnitude can be described in terms of the microscopic mobility p 


P.M. SUPPL.—JULY 1954 2B 
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which is the average electronic velocity per unit field strength. By 
definition j. is connected with the average relaxation time 7 by 


po=er/m. ok Ss ate ED 


The methods for calculating . or 7 are well established in the case of 
sufficiently weak interaction between electron and lattice vibrations, 
ie. when perturbation theory holds («<1). In this case the electron 
(or electrons) is considered as nearly independent of the polarization 
quanta. The electron is then accelerated by the field, and is scattered 
by occasional absorption cr emission of a quantum. Assuming the 
quanta to be in thermal equilibrium, the electron carries an average 
electric current which is proportional to the strength of the field. 

This procedure cannot be generalized in a straightforward manner 
to the case of stronger coupling («>1). For this would require the 
possibility of separating the system into free quanta on the one hand, and 
the electron plus ‘ bound ’ quanta on the other. The methods developed 
in previous sections, however, involved the assumption that the total 
energy (above the ground state) of the whole system is less than iw which 
excludes the presence of free quanta. The difficulties which one meets 
in attempts to go to higher total energies can be seen from the transformed 
Hamiltonian in the form (6.18)-(6.20). There if «?<100, the terms 
F,+ fF, can be considered as negligibly small provided no ‘ free’ quanta 
are present so that the wave function is given by x and the energy by 
E(k), as shown in (6.22). If, however, n, quanta of wave vector v are 
excited, as considered in (6.21), then the previously neglected terms 
F,+F, may become extremely large so that the separation of free 
quanta achieved in the Hamiltonian F—F,—F, has no longer any 
physical significance. 

In view of these difficulties application of perturbation theory only 
will be discussed in the present section. It may be hoped that some of 
the results obtained in this way have a validity beyond «=1, although 
no proof for this can be given at present. This hope arises from the 
calculation of the energy of the ground state in § 4 according to which 
perturbation theory happened to give good results up to «210. This 
means, of course, that electronic energies larger than fiw have to be 
excluded because in this energy region none of our methods could be 
applied. Perturbation theory does yield satisfactory results, however, 
for electronic energies which are very large compared with tw. This 
case was treated by the author (Frohlich 1937) in a paper in which the 
idea of polarization waves was first introduced. The probability per 
second for absorption or emission of a quantum by the electron was found 
very small compared with w so that no objections against the use of 
perturbation theory arise in this case. This no longer holds, however 
when the results are extrapolated to electronic energies near iw i 
agreement with the conclusions of the present discussion. If the electronic 
energy is less than iw, however, then quanta can no longer be emitted 
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but only be absorbed. If furthermore the temperature is so low that 
very few quanta only are present then the probability for this process 
is very small. It may be hoped then that perturbation theory may be 
applied to this case as also suggested above. The calculation of the 
mobility was carried out by Fréhlich and Mott (1937) and corrected for 
a more correct electron lattice interaction by Frohlich, Pelzer and 
Zienau (1950). In this latter paper it was indicated that after absorption 
of a quantum the electron has an energy larger than fiw above the ground 
state and hence can emit a quantum. The probability for this process in 
general is so large that perturbation theory cannot be applied, and this 
is the main reason for the uncertainty of the validity of perturbation 
theory even in the case of low temperatures and average electronic 
energies below fw. The fact that the probability of emission is very large 
suggests that the whole process of scattering should be treated similarly 
to that of resonance scattering of light. Thus a quantum fiw, wave 
vector v is absorbed by the electron and re-emitted with wave vector v’. 
In this process the electron is scattered elastically, changing its wave 
vector by v—v’. This resonance scattering has been treated by Zienau 
(1950).¢ This process can also be treated by replacing the free electron 
wave function by a function of the type treated in § 4, but treating the 
quantum as free, as carried out by Low and Pines (1953). Their result 
suggests that perturbation theory does give satisfactory results in the 
case of low temperatures and electronic energies below iw. Finally it 
should be mentioned that use of perturbation theory for both absorption 
and emission together with an appropriate treatment of the kinetic 
equation leads for low temperatures to the same result for the mobility 
as the method of resonance scattering. This follows from the paper by 
Howarth and Sondheimer (1953). 

The above survey shows that in actual crystals where usually 1<a«<10, 
there is some hope that application of perturbation theory for the calcula- 
tion of the mobility should give reasonable results for temperatures 
T <hw/k and electronic energies sufficiently smaller than hw though as 
mentioned it has not been possible, so far, to prove this. In view of the 
theoretical uncertainties, it seems desirable not to involve the complicated 
theory of resonance scattering but to derive the result very simply by 
assuming right from the beginning that the probability per second for 
resonance scattering, i.e. for absorption and re-emission of a quantum is 
equal to that for absorption alone as would in fact follow from perturbation 
theory where the ratio of absorption to emission probabilities is very small, 
of the order n/(1+n)=n<1. Here 

= 1 
le exp (hw/kT')—1 


is the average excitation probability of a quantum. 


Zire hot ae 2) 


+A mistake by a factor 2/3 in the final integration has been corrected by 
C. Herring (1952). 
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To derive this probability let the electron be in a state with wave 
vector q, so that its energy is iwq? and let 


elim oe ae ee 


Also let », describe the number of excited quanta with wave vector v 
so that the zero order wave functions are given by (4.2). Using hwF ing 
(4.7) as perturbation, transitions into electronic states q’ are possible 
in which a single quantum is absorbed, conserving energy and the total 
wave vector, i.e. 
q=q-++y, eee ee eS 
and 
qd 2=92+v2+ 2qv cosO=q?+1 or v?+2vq cos 6=1. (7.5) 


Here @ is the angle between q and v. Since | cos 6|<1, it follows imme- 
diately that v is restricted to the range 


1<~ged yh o> Roe eee 
Now since 
(ny—1]| 6, | n, =n, vt EO ee MD 


it follows from (4.7) and (4.2) that the absolute square of the matrix 
element for the absorption process is given by 


4a Ny 


| te al 


dics)? hmcban)4 eat ee ee 
The transition probability per second is obtained in the usual way by 
multiplying (7.8) by 27/h and by the number of levels in the final state 
per unit final energy q'*hw. Since the final electronic state is exactly 
given by (7.4) the density of levels is equal to the density of levels for the 
quanta, i.e. per unit energy using (7.5) it is given by 


S vdvdcosddd SS vdvdd 


Came OY ane ee ot! 


Hence with (7.8) averaging n, over a small v range, 


& ; 
BL a= 7 aeih a Ra ee Se a Ki ee (7,10) 
represents the probability of absorption of a quantum in a certain wave 
vector range by the electron in the state q. The wave vector range is 
such that the angle @ of q and vis exactly given by (7.5), the magnitude 
v has a range dv, and the azimuth ¢ with q as axis has a range d¢. 

It must be noticed now that conservation of energy requires that 
re-emission of a quantum by the electron in the state q’ leads to a final 
state q’’ such that Pee ; 

Y fie Sle a AM oot Serene eR PIE 
The direction of the final state is arbitrary and the probability for re- 
emission follows from (7.10) by replacing there x by 1-+-n. Now with 
(7.3) it follows from (7.6) that v has nearly the same value v1 whatever 
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the direction of q’’ so that the probability for scattering into any of these 
directions is equal. Thus using (7.10), (7.6) and (7.3) 


: aie . 

P.= | dP,=<— | “d eee op aed Omens (7:12) 
. amd Jo Ji-g V q l=q 

is the probability per second for any absorption process, and hence since 

the probability for re-emission is very much larger (provided 7 <1) P,, is 

also the probability per second for resonance scattering. This latter 

involves elastic scattering of the electron with equal probability for any 

scattering angle. Hence 


Prd Qi4ar=2eund Ql4r 6 >. 3) 46 (7.18) 


is the probability per second for the elastic scattering of the electron in 
a solid angle dQ. 

With the help of this transition probability the rate of change of an 
electronic distribution function f(q) can be calculated. For clearly 


of(q) 


=a = Pa | f(a)—F(9")) dQ" 4m, g?=q"?, . . (7.14) 


since P, is independent of the direction of q’’; dQ’ represents the 
differential of the solid angle of q’’. Inthermal equilibrium, f(q) depends on 
the magnitude g only and hence df/dt=0 in this case. It will be noticed, 
of course. that elastic scattering cannot alter the magnitude of q¢ and 
hence can establish equilibrium only if [f(q)d2 has the equilibrium 
value. This is assumed to hold in the theory of conductivity in weak 
fields where a distribution function f(q) of the form 


A(M=flg+@r/afilg) - - - - - + (4.15) 


is considered. Here q, is the projection of q into the field direction, and 
both f,(q) and j',(q) depend on the magnitude of q only. Clearly then 


| (fla)—f(4") 42” =(f@/a) | (ae—ae"") 42” 
=40(q0/D) filM=47(F(4)—folg)), -  (7-16) 


because fae” dQ’’=0. From the definition of the relaxation time 7, 


0 1 
EL SOS AT ee eee PC 
it follows with (7.16), (7.14) and (7.13) that 
) returns A SP OR AT) 


Also since 7 is independent of q, t=7. 

It should be remarked that second order field effects change the average 
electronic energy. Equilibrium can then not be established through 
elastic scattering processes only. 


360 H. Frohlich on 


The mobility » now follows by introducing (7.18) into (7.1). Use of 
(7.2) and (2.37) leads to 


/2 
ee a eto (Fe) fexp (heofkT)—1], kT <hew .- (7.19) 
(2ham,))'*\en  € m 
where . 
dg==h?/me2=0-5X10-Sem wwe. (7.20) 
is the Bohr radius. Hence if 
O=ha/k sf et eee ON ee See) 
is expressed in degrees absolute, 
=i 3/2 ra) 
(Fla) Be il PANES oases volt sec, T<O. .. (7.22) 
ene ™ /@ 


It will be recalled that w/27 is the frequency of longitudinal vibrations 
which is connected with the frequency of residual rays by eqn. (2.10). 
Hence if c is the velocity of light and A is the residual ray wave length 


then from (7.21) and (2.10), 
h 2nc € 
se fs ae, ee 
@ fev Jz ( + 


Before comparing the theoretical value (7.22) for the microscopic 
mobility , with experimental results the above mentioned difficulties _ 
in the theory must be emphasized again. Furthermore it should be 
remembered that there will always be contributions to the scattering of 
electrons due to mechanisms of a different nature from the one discussed 
here. Thus scattering due to the so-called acoustic vibrations will make 
a contribution to ~ proportional to 7'~*/? and should thus be decisive at 
very low temperatures. The magnitude of this type of scattering is in 
general very small, however (cf. Mayer and Polder 1953). 

Experiments on NaCl have been carried out by Redfield (1953) who finds 
p=250+50 cm?/volt sec at T=82+2°K and w~ 40 at T=200°K 
Theoretically with «=5-6, «,,.=2-3, A=61 x 10-4 cm, one has 9=340° ae 
hence w= 280(m,,/m)?/? and 4~=15(m,,/m)>/? respectively. The agreement 
is very sera (with m~ m,,). Measurements on AgCl by Haynes 
and Shockley (1953) can be represented by p=30[exp (290/7')—1]. 
- Theoretically with «=12-3, «,,=4-0 and @=290°, the factor 30 is replaced 

by 4(m,,/m)?? which has also the correct magnitude. The measured 
temperature range is not wide enough, however, to exclude a 7-3/2 law. 
Agreement with the theoretical temperature dependence over a wide range 
was obtained by Breckenridge and Hosler (1953) on TiO,. Their 
absolute values indicate very high electronic effective mass probably 
due to conduction in impurity bands. 

Perturbation theory which was used in the derivation of expression 
(7.18) for 7 is also valid for electronic energies which are very high 
compared with hw. The actual expression in this case differs, however, 
from (7.18) for various reasons. This case will not be Sonaidered here 


ee le ie 
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because interaction with very short polarization waves is predominant 
in this case so that treatment of the lattice as a continuum is no longer 
permissible. This leads to some uncertainties concerned with the value 
of w and with the correct form of the interaction, but on the other hand 
use of perturbation theory is well justified. In contrast to this case, 
however, use of perturbation theory for electronic energies of the order 
hw is quite unjustifiable in particular for alkali halides where «6. 

Further experiments might help to overcome the theoretical difficulties 
encountered in this respect. This would require the measurement of 
electronic mobilities for electronic energies of the order of fiw i.e. in the 
case of thermal equilibrium for temperatures of the order 770. 
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Bent and sectioned to show slip lines. 


Subgrain growth in bent silicon—iron single crystal as a result of va 
annealing treatments. 
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PAUL A. BECK Phil. Mag. Suppl. Vol. 3, No. II, Pl. 27. 


Subgrains in pure Al after creep, as shown with oxide film-polarized light 
technique. [Ref. (11).] 
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Fig. 19 


Macrophotograph of cleavage face of indented Zn crystal, as in fig. 18, showing 
substructure formed on indentation. [Ref. (16).] 


Migration of high angle grain boundary during recrystallization in pure Al. Two 
successive boundary positions in specimen rolled 7-5° and annealed at 500°c for 
25 sec, resp, 32 sec. Oxide film and polarized light. 75. [Ref. (11).] 
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Fig. 51 


Photomicrograph showing partial coarsening in Al+1-1% Mn alloy after final 
anneal of 1665 minutes at 625°c. [Ref. (118).] 


d boundary migration in pure Al on 


Photomicrograph showing strain-induce 
© after 14% deformation by rolling. 


annealing for 5 minutes at 400° 
x75. [Ref. (146).] 
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Fig. 58 


Photomicrograph showing five successive stages of strain-induced boundary 
migration in pure Al on annealing at 350°o for 5, 10, 20, 30 and 150 
minutes, after 12° deformation by rolling. 75. [Ref. (11).] 
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Photomicrograph of pure Al specimen showing many small recrystallized 
grains of various orientations formed in the immediate vicinity of the 
two surface scratches. Those recrystallized grains with favourable 
orientation with respect to the 12% rolled matrix grain grew selectively 
on annealing for 1 hour at 350°c. Electrolytic polishing, anodic etching, 
polarized light. 25. [Ref. (55).] 


Fig. 62 


Photomicrograph showing several small recrystallized grains of varying 
orientations nucleated at the sheared edge. Those of the recrystallized 
grains with favourable orientation grew large at the expense of the 
annealed fine grained matrix with single orientation texture. Matrix 
was produced by 3 minutes anneal at 525°c of pure Al single crystal 
previously rolled 90%. Deep etching. x8. [Discussion of Ref. (157).] 


Wy 


EDITORS 


: Ss SS D. McKIE, D.Sc., Ph.D., 


University College, London. 

HARCOURT BROWN, 
M.A., Ph.D., 

Brown University, Providence, R.I., 
& S.A, 
yuna H. W. ROBINSON, 
Former Librarian, 
Royal Society of London. 
N. H. de V. HEATHCOTE, 

of B.Sc., Ph.D., 


University Coliege, London. 


Science ANNUAL SUBSCRIPTION 7 
£3 3s. Ou. 


Ww OF OR 
QUARTERLY Kee 18s. Od. 
ORY OF SCIENCE PER PART 
THE HIST CE POST FREE 
AISSAN 


I ALERE FLAMMAM 


THE DREATHEMATICAL WORKS 


OF JOHN WALLIS, D.D., F.R.S. 


by 
J. F. SCOTT, Ph.D., U.A. 


‘His work will be indispensable to those interested in the early history of 
The Royal Society. | commend to all students of the Seventeenth Century, 
whether scientific or humane, this learned and lucid book.’’—Extract from 
foreword by Prof. E. N. da C. Andrade, D.Sc., Ph.D., F.R.S. 


Recommended for publication by University of London 


12/6 ve 


Printed and Published by 


VEAYLOR Sac, LTD. 
ee EE 


